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SYNOPSIS 


This. investigation is intended to show that w sama noise is a Pel a repre- 
os ot of earthquake motion. It appears that a white-noise source used in 
conjunction with an analog computer is a convenient method of analyzing 
_ structures ‘subjected to complex ground ‘motions and of assigning probabilities 
to the deformations arising. _ Until more statistical information is available a 
about earthquakes, ‘it suggested that a “standard large 


0.75 ft2/sec4 /cps anda duration of 1/2: min. 


I ION te 


in os lack of hevebetes of destructive earthquakes. Accelograms of a few | 
“large earthquakes are available and these records show very couplex ground 
motions emtirely different from each other. _ They ar are different in energy con- Pe 
tent or size” and different in duration and shape. Earthquakes of the 
Bin —Discussion open 1, 1960. To extend the closing date one 
- month, a written request must be filed with the Executive Secretary, ASCE. This paper — 
_is part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings _ 
of the American Society of Civil Engineers, Vol. 86, No. EM 2, April, aR 
Research Officer, Physics Sect., Div. of Bldg. Research, 
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1 - will not necessarily bear any relation to those in the past and no prediction of 
a the severity and character of future earthquakes can be made. Two earth- : 
- quakes having the same energy content will have different effects on any par-_ 
ticular structure. Each individual earthquake, of course, has a specific 
: character and a specific ons a particular structure, but this character is 
p Bias Many investigations of structural vibration due to ground motion have been 
“made on a deterministic basis, using actual records of earthquakes a and de- > 


termining the deformations that would occur in various types of structures. 
This approach would be satisfactory if a large enough number of records were © 
vailable and the results displayed on a probability basis. Ut Until = is pos- - 
sible recourse to another mode of attack must be made. 
_ Earthquakes originate as ill- defined slip-faults followed innumerable 


; it travels. “Intuitively, one feels that the essential characteristic of the earth-_ 

quake is its irregularity. Many writers2; 3,4 have suggested simulating earth-_ 
quakes by a random assemblage of impulses. In order that white noise sive 
an earthquake it must have the ‘same effect ¢ on the kinds of 


LINEAR ELASTIC STRUCTURE 
he The response of a on structure of many degrees of freedom to ) groun 
motion may be written in of of that ground motion. om Lat 


period of the e nth harmonic; 


“weighting” factor a and is a structure 


= = fraction of critical damping of the nth harmonic ay ee 


Eq - 1 states that the response of : seaeldnasiiaae -of- freedom structure is a 


excitation occurs only in the integral factor and | consequently | this is he only 


wee “Response Spectrum Techniques in Engineering Seismology,” by D. E. Hudson, Pro- 
World Conf. on Earthquake Engrg., San Francisco, Calif., June, 1956. 


by E. Rosenblueth, 


«Response of Tall Buildings to Random Shakes, sf by A. _Cemal Eringen, Tech. Re- 
Office Div. of Purdue Univ., December, 
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‘The velocity spectrum of the ground motion is defined as the maximum — 


value of this integral, is, 

is a function of period and damping only. 

Spectra of several large earthquakes have been calculated and are well Py 

— 5 The spectra are all different in detail but have the feature in common = 4 

fluctuate ‘rapidly, with respect to period T, about a mean line which 
starts at the origin and rises asymptotically to a constant value. G. W. S 


show the shape of these common average spectra, and the scale of 
_. shown is that of the largest earthquake recorded, viz., El Centro, Calif., 
i The writer used an analog computer to construct mean spectra of bursts” 
of “white noise” for comparison with the mean spectra of Fig. % % . “White — 
: noise” is defined as a random signal having a gaussian probability of ampli- 
i. tude distribution and a spectral density which is constant for all frequencies. 
The ‘Spectral density at at any particular frequency may be regarded as 
average power | passing , when the signal is filtered by a narrow band. pass fil- 
_ ter centered at the particular frequency. In practice, , of course, one en- 
counters” “band - “limited white noise,” that is, the spectral ‘density is only 
constant up to a certain frequency and is then rapidly cut off. As long as this q 
“id frequency is large compared with the frequencies of the system being con- : 
- sidered, the response of the system is unaffected by the value of the cut -off 
~ - point. The computer used provides a a band-limited white noise having a ‘spec- - 
tral density flat from 0 to 35 cycles. 


The equation of men of the mechanical oscillator subjected to an ac- 


in which x is the displacement of the mass relative to its zero position a 
WwW is the angular frequency of the oscillator. _ Variables in the computer are — 
_ represented by voltages having an allowable maximum of 100 volts. If xis to r 
have an expected maximum of V ft amis sec the mechanical al equation is scaled — 


here the carens in square brackets are the respective computer voltages. “—_ 


5 “Spectra Analyses of Strong Motion Earthquakes,” by J. Cc. ‘Alford, G.W. 
and R.R. Martel, Earthquake Research Lab. b Tech, Report to Office of Naval Research, _ 
Inst. of Tech., Pasadena, August, 1951. 


° “Limit Design of Structures to Resist Earthquakes,” by G. W. Housne 


Proceed- 
ings, First World Conf. on Earthquake Engrg., San Francisco, Calif., June, 1956. 
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EARTHQUAKE REPRESENTATION 


wit) is actual of the noise source and the computer 


to ere the factor will We (t) per sec2 
a> The actual velocity x in feet per second is given by the : computer volts 
multiplied by v/100, the acceleration by X x Vw /100 and the on i) 
Ideally, a probability distribution should be given to the size ‘and duration 
_ of the earthquake as well as to its configuration. However, with the available 
information, this Ss cannot be done ‘satisfactorily. The available records of large 
earthquakes show that their duration was generally of the order of 20 sec 


. ‘For various values of period T and damping A, twenty separate bur bursts of re 


white noise of 25 sec duration were fed into the circuit. The system was re- 


set to zero after every burst. Each burst of noise e gives a different value ~~ 
the maximum velocity acquired by the oscillator. pie 


lle: The maximum value of x during each burst of noise is the required spec- = 


trum for that particular burst of noise. The mean of the maximum velocity of 
_ the twenty runs was taken in each case. . This gives the variation of the mean d ia 
maximum velocity for bursts of white noise as a function of period and damp- 

However, this mean still contains the scaling factor k, and the 
given to k will determine the size of the earthquake being represented. — This 

Variation is shown by the circles in Fig. 2c The factor k was chosen so that 2 


the: scales matched at an arbitrary point viz., T = 3.0 and A = 0.20. With this 


__ variation of the ke mean spectra of t the noise with both period and damping is 
very similar to. Housner’s common average earthquake spectra. The scales 


_ value of k the position of the other points was determined. It is seen that E | 
cond have been matched at any other point but the value of k would not be oa 


From this it appears that “white” ground accelerations ‘suitable 
representation of earthquakes. 2s Further, by choosing the scales so that the 


-mean-noise results fit the averaged spectra of the largest earthquake | 


corded, we arrive at some kind of a a definition of a a standard large earthquake. © 
_ Individual bursts of noise will give . results above or below the average just os 7 
7 the actual spectrum of El Centro fluctuates above and below 1 the « smooth curve 
ae bye: The spectral density of the noise source used was approximately 2.90 vt?/ 
Se cps and the factor k came to 0. 51. 


tentiometers, triangles are summers, and double - based 
It is to compare the ‘effect of white noise and an actual earth- af 
on a non- -linear most type non- linearity that 
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EARTHQUAKE REPRESENTAT rION 
7 2 will occur in a structure during : severe > ground 1 motion will be plastic y yielding. 
_ ‘Fig. 3 shows the assumed relation between restoring force F and displace- 
=~ .2 x. Starting from zero with positive velocity, the restoring force is pres 
portional to displacement ‘until the yield point A corresponding to a displace- 
_ ment Xp is reached. The restoring force remains constant at this value until © : 
a change in sign of the velocity occurs at some point B. The hos geese . 
A B will be called a yield displacement. ‘The force F now decreases along a a 
line parallel to the elastic part of the characteristic and stays on this line un- 
til C or perhaps B again is reached. If the velocity stays negative the reverse 7 
yield point C is reached and the force stays constant at this opposite value : 
until the ° e velocity becomes positive again and the characteristic proceeds from 
D to E. Areas enclosed by loops represent hysteretic energy loss. 
- Ina recent paper” the writer considers the yield displacements that occur _ 
when an oscillator with 1 this | ‘yielding : stiffness characteristic is subject to the 
El Centro earthquake. These yield displacements determined for this 
specific earthquake are now “compared with the yields arising when the same 
system is subjected to bursts of “white” ground acceleration of spectral 
density 0. 75 ft2/sec4 /cps simulating the standard large earthquake. aes ae 4 
 ‘The- equations of motion may be derived from a consideration of the 
-elasto- -plastic characteristic shown in Fig. 3. There are separate equations _ 
depending | on whether the e motion is taking place in the elastic or the plastic ia 
part of the characteristic. On the plastic part of the characteristic the re-— 
tarding force is constant and has the sign of the velocity shown by sgn x. On 
elastic part of the characteristic the 1 restoring force must be referred 
the position caused by the yielding. of motion 


"where wis the algebraic sum of all previous yield displacements up to the 


considered and the sign of the yield displacement taken as that of k at 


_ the time. Thus, in Fig. 3, when working on the linear ' portion D E, the value a 


. It will be more convenient now if the yield point is specified in | terms ot ; 
equivalent acceleration rather than the displacement Let G be the: equiva- 
’ is the expected maximum value of x, and X is the scaled computer varia- > | F 
corresponding x, then the scaled equations become: 


= 


+2 + 


| 

= 
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eis ‘These equations were simulated on the computer by the circuit diagram 
shown in Fig. 4. The stiffness characteristic is provided by the | diode - limited 


integrator. As =: increases initially the output of the integrator is --X until: aa 
limited by the diode. ‘If K 4 continues positive the voltage output stays at this at 


constant value. However, ‘if x changes sign, the integrator starts subtracting 


= a linearly- increasing quantity and thus provides the term (X-U) until — 8 


its initial zero at any time. This displacement consists of the elastic rom 

_ placement xg added to the algebraic sum of all previous yield displacements. aa 
;, The maximum value of the algebraic sum of the yield displacements will be 
called the maximum displacement and is the maximum 


RECORDER” 


FIG. _4.—COMPUTER DIAGRAM FOR ELASTO-PLASTIC OSCILLATOR. 


ie value G, , the mean of the maximum yield displacements a arising was deter- = 


% ‘mined when the system was subjected to ten runs of 25 sec duration. Fig. 5 to e 3 
&§ , 13 show this mean maximum yield displacement together with the elastic dis- . 


placement : X, . An addition of the two gives the mean maximum total displace- 
ment. The | ull lines are from the paper mentioned earlier and refer to the 


1940 El Centro earthquake. The results from the white noise are shownas 


‘the 

_ andusing this value, it is "le that the noise gives similar results as an actual 

The double- peaked “shape ai appearing in the cases shown in Figs. 8 and 9 is 
¥ er: reproduced by the noise and thus must be an essential function of the sys 
tema and not merely of the particular character of ‘the 
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aS ke burst of noise is of the same energy level but gives a different value of a 


EARTHQUAKE REPRESENTATION 

These of course, compare the mean of the maxima o 

noise with the maximum of one specific earthquake. Ideally, the ‘comparison 


should be with the average of several earthquakes of equal spectral density. _ 


{ ten 


maximum yield lying above or below the averaged values shown . Presumably — 
one individual burst of noise could be found which would correspond to this one 
‘The probability distribution of the maximum yield w was determined in three 
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four times” the average. This: easily accommodates the shown 


between the average noise values and the one earthquake considered. 
“White noise” appears to be a reasonable of earthquake 
for engineering design purposes. ‘The technique of simulating a complex 
structure on an analog computer and feeding it with bursts of noise having an % 
energy level representative of a large earthquake is a straight forward method — : 
_ of analyzing a structure and of assigning probabilities to the various displace- 


a. _ Using this technique it will be possible to simulate multi- -story OUP 7 


_ The probability values given in Figs. 14 to 16 illustrate an important ce. is 
es burst of noise is of the Same energy level but the resulting displace- 
> - for individual shots can differ by factors ‘of 10 and more. Thus two 
of the same content energy or magnitude but of different 
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any Boe) Cases Irom one hundred runs Of noise. ese are shown 1n Figs. 14 to 16 ee 
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one may seriously affect a structure with a a certain set 
: of parameters but have little effect on another one. This may well 4 path of ee 


‘the explanation of the inconsistencies noted in seismic dam: ae q 
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_ By George Gerard 


structures are presented. The first part of the paper is with 
compressive - strength analysis of long flat plates after elastic buckling 
curs. The plates have various boundary conditions along the unloaded edges. 

‘The “results of a semi- empirical analysis of a large mass of test data are 
Ps shown to apply to various forms of composite box-beam structures. 
ae second part of the paper considers the design rather than the analysis 7 


= of composite stiffened ‘plate structures subject to instability. 

principles of minimum -weight design are presented for transversely support- 
eg ed, longitudinally stiffened plating. it is indicated that designs based on these 
constitute efficient forms of limit - structures. 


fs ‘The field of structural concerned with elastic and slestic sta- 
_ bility has reached a high level of development asa result of important contri- 
butions from various engineering disciplines. literature of this 


— the years provides ample evidence that civil, aeronautical, and naval 
Ls 
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es of | the American Society of Civil Engineers, Vol. 86, No. ‘EM 2, Apri, 1960. ee. 
Presented at the 19 
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has focussed nee of minimum-weight on the thin- -gage metal 

tures. Consider that approximately two-thirds of the structural weight of an 
airplane is designed by stability criteria and its importance becomes readily 
apparen This has naturally provided a powerful impetus for a concentrated fs 
research effort in this area. 
In the spirit of continuing to share knowledge in this field, some recent de- 
velopments are presented which are taken mostly from research in the aero- ; 
a field ee the past decade (the 1950’s). It is hoped that these ——* 


hin-gage metal structure 


‘COMPRESSIVE STRENGTH OF PLATE 

ather well known that the compressive strength of att failure 

P, can significantly exceed the buckling stress when the latter occurs elastically 

. e In the plastic region, however, the buckling stress and failure strength are es-_ 

4 sentially coincident. The theory of plastic stability is now ina 

state of development and the ; 


The compressive or crippling strength of flat plates is of particular im- 
portance in the economical | design thin- “gage metal structures. Some r rather 
mass of test data,¥, 
Approximate Strength Exact of the crippling strengthof 
flat plates presents a complex mathematical problem . This is a result of ail 
: non-linear, large deflection behavior of the plate after buckling, in | addition t 
the nonlinear plastic behavior which ultimately results infailure. __ a one 
itis is 3 noteworthy that E. Zz. Stowell? and J. Mayers | and B. Budiansky® have 


In both analyses, the plates were assumed to have supported unloaded eters a 
that do not warp but remain straight in the plane of the after as 


esa ‘value approximately equal to the compressive- yield strength of the m 


Because of the crippling-strength problem, it is con- | 
venient to resort toa a “strength of materials” is ‘semi- 


« Handbook of Structural Stability, Part - Failure « Plates ai ond con El 


ments,” by G. Gerard, NACA TN 3784, August, 1957. 


, 4 «Handbook of Structural Stability, Part V - Compressive Strength of Flat Stiffen 

Panels,” by G. Gerard, NACA TN 3785, August, 1957, 
a -° “Compressive Strength of Flanges,” by E.Z. Stowell, NACA TN 2020, January, 1950. fax’ ee. _ 
“Analysis of Behavior of Stingly Supported Flat Plates Compressed Beyond the 
Load the Plastic Range,” by J. Mayers and B  Budiansky, MACA TN 3368, 
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STRUCTURES 

a the following relation between the ¢ average compressive stress, 6 
and the ee unloaded edge stress, Op, is valid: = =~ 


it” is now w assumed that failure occurs the stress. reaches the 


3) 


5 contains the non-dimensional parameters and (t/b)(E/oey /2, 
which canbe used to synthesize test dataon the cripplingof flat plates of vari- 


= types of materials. The constants m and 8 are determined by a best - fit of 


-Stothe test data, 
"Flat Plates.— There is a considerable mass of test data from many } investi- 


gators on | the crippling strength of long, simply s supported flat plates. These 
data have been summarized3,4 | and include plates tested in V- roove support-_ 


a ing fixtures at = unloaded edges as well as square tubes which act as an as- 


“a sembly of four, ‘simply supported plates. 


-: In both of these cases, the unloaded edges are free to warp after a Se ee rot. 
%s and thus the stress distribution resembles that shown in Fig. 1(b). This 


¢ K - difference in edge conditions between the tested plates and those assumed in 
Big Pe theoretical analyses of Stowell5 and Mayers and Budiansky,® and shown in 
y 
voll Fig. 1(a), results in a significant difference in crippling —" The — 
with warped edges generally have a lower strength. 
. om Available data on V-groove plates and square tubes, of various sinilieneats 

various materials which differ significantly in E and Also fits the 


data within + - 10% limits the 1e vente 


— 
4 iii 
_inthe plate 
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4 
In addition to the data on simply supported plates, there are some test data 
on long, clamped flat plates that can be used to check Eq. 6. These data were — a 
obtained? by K. B. Jackson and A. H. Hall who used a special test fixture — - 


_ provided clamping along the ieaiedae unloaded edges and also permitted i in- 


FIG, ~—8T RESSES IN FLAT PLA’ TES AFTER BUCKLING UNDER 
eit CONDITIONS OF UNIFORM END SHORTENING, (a) STRAIGHT © ar 
UNLOADED (b) STRESS-FREE UNLOADED EDGES, 


al The data in eernd of the e parameters ‘s of Eq. 5 are shown in Fig. 3. It can i 
_ observed that the value of m = 0.85 fits the data well Bs value of 8 was pre- a 


dicted using Eq. 6 in the form 


Since for the simply s supported case k = ‘ and for the clamped 


7 kg = 6.98, it dicted by E ‘8 that for the lamped plate “4 
case kp it is predit cte 


It can be idoied ) in n Fig. 3 that this value of 82 does, i i 


ft to the test data. 
Actually, the data pre: in Figs. 2 and3 are a sample of 


_ mass of test data.? 38 The excellent correlation obtained therein, for many 


a pe “Curved Plates in Compression,” by K. B. Jackson and A. Hall, National — Re- 
of tage by G. Ger: ard, Journal Aeron- 
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‘FIG, 3. DATA ares PLATE ELEMENTS IN COMPRESSION. 
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April, 1960 
materials, of plate elements, and temperature conditions, lends 


considerable confidence to this semi-empirical approach. 


ong tat pat results can now be summarized as follows: ane cr rippling strength 


limits up to the 

The values of B gi given in Table 1 are based on in? and 8 


_ Plate Structures. —Iti is of interest now to consider the application of the cata 
i concerning compressive crippling strength on plate elements to stiffened plate 
structures such as box beams in bending. It can be anticipated that under cer-_ 
conditions, the crippling formulas presented previously may be useful 
i predicting the strength of relatively complex plate structures in a rather simple 
mr. V-groove supported flat plates, and tubular assemblies of plates achieve 35 4 
strengths inexcess in excess of the elastic elastic buckling stre stress byv virtue | of 


TABLE 1.—CRIPPLING COEFFICIENTS FOR 


ae to reach the compressive-yield strength. - In a plate structure, « on the © 


4 
E other hand, various types of stiffening elements are used to support the unload- 
ed edges. If the unloaded edge-stiffening elements are strongly attached to the 
= so as to form a monolithic structure, and the deflectional restraint of 
the stiffener is sufficiently highso as to remainstraight up to the compressive ¥ 
_— of the plate, then the crippling strength of the plate can be achieved. In 
: such cases, the stiffening elements of the plate structure effectively perform | 
the same role as the ii -groove supporting fixture for the plate | elements or the 


__ The bases for these remarks are contained in the work8;9 which demonstrates 7 
the application of the crippling- strength formulas to stiffened panels u under 

_ compression and multicellular box beams under bending. Test data on the 

2 latter may be of particular interest and, thus, some of the test results taken — 


from various sources and presented by writer® are summarized in the 


ge 


= G. Gerard a rd and “Becker, NACA T TN | be 
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The tests under consideration were conducted on aluminum alloy box notin 


in bending. | Some of the b beams were ‘single- -cell tubes while most were fabri- - : 
cated boxes of 1, 2, or 3 cells” s with various angle attachments at the web-cover | 
joint. _ Depending on the web geometry relative to that of the plate elements of _ 
_ the compression cover, as well as the joint details, the cover -plate elements | 
_ received elastic restraint along the unloaded edges in some cases. Thus, the - 
- test data presented in Fig. 4 are divided according to values of buckling coef- a 
- ficient, k = 4, 5.35, and 6.35 pertinent to the test. In Fig. 4, b is the width of 
_ the cover plate between web stiffeners, E and ene are the mechanical proper- 
“ties: of the plate material, and the c crippling strength is computed from G = M c/I, ; 
where the cover material is assumed to be fully effective. 
It can be observed from Fig. 4 that for the k = 4 case which corresponds to 7 
simple support, Eq. 10 fits the data 1 very well when using the value of 8 = 1. 42 ae 
a in Table 1. The lines for the other two values of k shown in Fig. 4y were — 
computed by use of Eq. 8: for k = 5.35, 3 = 1.60 and for k = 6.35, 8 = 1.75. In| 
all cases, the use of Eqs. 10 and 8, basedon plate elements 5 under compression, 
‘predicted the behavior of the multicell box beams under acceptable 


 —- test of the predictive value of Eq. 10 is afforded by the sins summa- 4 
-rizedl0 ’ by J. Vasta on tests of steel ship hulls to failure. One hull was tested | _ 


in the sagging condition and experienced a compressive failure of the deck. — 

This is shown as point 1 in Fig. 5 where the lines correspond to Eq. 10 for 
various values of k given in Table : The relatively light framing used in the br { 
region c of the deck failure did not provide any significant elastic restraint bee : 


sulted in compressive failure of the bottom plating adjacent to the keel. .In both — 


a cases, the plating was bounded by relatively deep longitudinals which apparent 
: ly provided significant elastic restraint along the unloaded edges of the plate. — 
Thus, points 2 and 3 in Fig. 5 correspond to k of approximately | 6, which ap- | 
 m From Figs.4 and 5, there appear to be some positive indications that ‘Eq. 1-10 
can be applied to plate structures. It must be remembered, however, that the — 
_ method is semi-empirical in nature and requires verification in cases whit which | 4 


on may depart eer from the tests on which it is based. an to 


4 _ The other two hull tests were conducted for the heating condition and re-_ 
4 


principles, which require simultaneous yielding of the “gtructural 
elements. Alternately, minimum-weight designs, according to Stability con-— 
cepts, the most economical form of limit design, ~~ 
i: _ There is, of course, a long history of civil engineering interest and develop- 
-ments in economy of design. Particularly in recent years, modern developments" 
in plasticity theory have resulted in renewed interest in the economical design f 


« 


_“*“ “Lessons Learned from Full-Scale Ship Structural Tests,” by J. Vasta, Soc. of 
Naval Archts. and Marine Engrs., Preprint No. 1, November, = 
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In the design of thin-gage metal structures under compressive loading, sta- 
— bility rather than yield considerations generally constitute the limiting strength _ eee 
2) conditions. erefore, in such cases, economy of design must be based on [i 
— 
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4. _—DATA (ON BOX BEAMS UNDER BENDING WHICH FAIL BY CRIPPLING 


OF THE COMPRESSION SURFACE. 
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PLATE STRUCTURES | 
“of beam and shell structural elements. however, exclude sta- 
bility as a mode of failure and, consequently, apply only to relative thick-walled 
—— Because of the interest in thin- -gage structures, ‘it is pertinent to 


turn nov now y to some of the developments in the aeronautical field concerned with» 
the minimum- weight design of such structures. 
Design Conditions. —In the design of a structure, certain of the overall di- 4 
e- are, generally, fixed by considerations over which little choice can 
. be exercised, such as span or panel width. However, as the structure is divid- 
ed into components by the introduction of longitudinal and transverse stiffening — 
elements, it is then possible to influence the weight of the s structure by the spac- 
In addition to these geometrical considerations, the overall loading features” 
are, also, generally fixed by design conditions over which the structural de- 7 
- signer has little | control. These factors are governed primarily by the per- 2 5 
formance requirements of the structure. 7 Thus, in approaching the structural 
design, the overall geometrical a and loading features are essentially ieee 4 
7 By combining the loading and geometric c parameters into a suitable lading 
in index, it is possible val conduct minimum-weight analyses of various types of 
compression elements. A relatively wide > variety of such cases is available, 11,12, 5 
based on appropriate indices. For the purpose of presenting a review of some 
of the essential aspects of minimum -weight design, the analysis of a wide, 
longitudinally stiffened panel is developed in the following. These results are 
* then combined with those for a lateral-stiffening system in order that the struc- 7 C 
tural efficiency of ‘composite compression surfaces can be considered. cs , 
_ Wide Longitudinally Stiffened Panels .—In view of the large possible number 
of buckling an and failure modes ofa wide longitudinally stiffened panel, it ap- 
pears that any analysis aimed at predicting the strength might be limited in 
- scope. It is remarkable, however, that such analyses are in substantially good 
agreement with test data. In an attempt to arrive at conclusions concerning 
_ the major factors to be considered in the design of efficient, wide, longitudinally 
stiffened panels, a simplified analysis first presented13 by A. Zahorski leads 
The assumed |conditions for minimum- -weight design of a compression s struc- 


ture are that the applied stress andallowable stress be equal and that the po pos- wer 
sible of buckling occur simultaneously. Thus, the applied 


_ It is to be noted that the loading | N may contain appropriate safety factors. _ 


Y For the wide, stiffened panel, it is assumed that the flexural type of insta- 


(11 


“Weight- -Strength Analysis of Aircraft Structures,” by F. R. Shanley, McGraw- 
_ 12 “Minimum-Weight Analysis of Compression Structures,” by G. Gerard, New York 
13 «Effects of Material Distribution on Strength of Panels,” by A. Zahorski, Journal 
Sciences, Vol. - 11, No. 3, 1944, pp. 247 
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April, 1960 
stiffened panel is then considered as ane nee of plate elements sand 
type of | buckling of these elements is given 


ts 


iy 


en tes 


; and denoting the stress under the conditions of 7 16 as. 6, Eq. 15 can be re- 
= pts 1/4 NE 1/2 
(18) 


£ generally f fixed by those considerations which govern the structural loading in- — 


dex N/L. - _ Consequently, ‘it is possible to arrange the material in | the cross © 
B Section to arrive at an optimum distribution that results in the maximum value — = ; 
for the configuration and, therefore, optimum stress according to Eq. 18. ee 
a commonly used longitudinal. stiffening arrangement is shown in Fig. 6. | 


= and te/tw = 1. 0, the are the ‘geometric 


ofa) __ On the condition that the applied stress and various buckling stresses are -- 7 
“| 
Panel Efficiency.—The_panel-efficiency coefficient, ap), contains terms, 
| 


(0. 633 + 0.37 


5. 19 and 20, Tb denotes and rt represents By combining 
Bas. 19 and 20 with 18, the e panel- efficiency the stiffened 


kr? 633 + 0. 37 ry 


An approximate method of analysis, which yields excellent results for the _ 
- Z- atiftened wide panel, is based on the assumption that the various elements © 


substituting ‘into 21 and taking k= 4 for pin- jointed “webs 


[ret (0.033 633 + 0. 
Upon maximizing 0, , the results are caine: 


and the results of an extensive test program conducted by the National Ad- 
_ visory Committee for Aeronautics12 is presented. In observing Fig. 1 it must 


4 


_ The discrepancy in the region of N/L = = 0.2 ksi is virtually eliminated when 
the proper value of the plasticity reduction factor, 7 in Eq. 8, is used in the 
region: between the proportional | limit and the yield strength. 

_ By substituting the values given by Eq. 26 into Eq. 19, we ca can deduce that 
the optimum distribution of material in the longitudinally stiffened panel is 
415 t in the plating, and 0.585 f in the longitudinal stiffeners. 


It is interesting to note that this method of analysis c can be extended toother 
cross-sectional configurations. J. Farrar has given!4 results based on a 


a 14 “The Design of Compression Structures for Minimum Weight,” by D. J. F 
Journal Royal Aeronautical Soc., Vol. 53, November, 19 1949, 1041- -1052. 


stresses for the skin and stiffener web. Hence, 
— 
24) 
2) 
| 
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somewhat more exact method of analysisfor Z, hat, and Y stiffened wide —_ 
In practically all cases, the results for common stiffening elements = 
that ap = 1. Again, as shown by the writer, 12 excellent agreement is obtained 
_ between theory and test data for the other stiffener r shapes. 


Transversely Supported Panels | Optinaem Design. —In the 
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“the loading index, N/L. Since L is the mie of the panel between tr ansverse 
<a 


= which results inminimum weight of the composite compression system, con- 
sisting of the wide panel and the transverse s 


— 
— 
— 
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6. LONGITUDINALLY STIFFENED WIDE PANEL 
i 
FIG. 7.—COMPARISON BETWEEN TEST DATA AND THEORY FORWIDE, 
LONGITUDINALLY Z-STIFFENED PANELS OF MINIMUM WEIGHT 
hne the ramaining nart af thic nrohlem ic tan determine the ontimum cnacing 


‘This type of analysis has been presented in rather complete detail.’ For 
purposes of illustration of the | of the essential results are pre- 


he minimum deflectional r rigidity of the simply Supported transverse stif stiff- - 
ener that is required to nodes in the atiftened and, then, 
4 


general instability, is given!? 


combining ant 28, the e effective thickness of the tran r 


verse stiffene 


Since t = N/o, om Be u, and by substituting Eq. . 18 for ¢ and Eq. 29 for ir, 

=. 

_ Eq. 30 can be transformed into the following non-dimensional form for T= 1 


is the effective thickness of the composite structure. j= 
Since the minimum- -weight design is of interest, Eq. 31 can be minimized — 


terms of the transverse, stiffener spacing L by performing the 


for ‘minimum weight 


The resulting effective, 


It can be observed Eqs. (t and that the optimum spacing 


LY 


and optimum effective thickness (te/w o) are in terms of the geometric pro- 
a portions of the transverse stiffeners (w/h) and the loading index (N/w). ). Thus, di 
= original objective of determining the solution in terms of the prescribed 
loading (N) and geometric quantity (w) has been achieved. 
a Be: In Eq. 33 the terms in the last bracket represent the effective thickness of i 
_ the stiffened panel and that of the lateral stiffeners, respectively. Thus, the — 


— 
a _ For the long, transversely supported stiffened panel shown in Fig. 8, the [am 
 andc=1: 
— 
— 
= 


‘the analysis: for transversely stiffened panels of optimum | phere it is of = 


FIG, 8. _—TRANSV ERSELY SU PPORTED LONGITUDINALLY = IFFEN 
panels xs for a yield limit rather than a stability limit, and, also, 
to illustrate the difference in weight between these two design conditions. 
a wide longitudinally stiffened panel with effective thick- 


the e stress level is fixed at the yield the transverse 


Eqs. . 34, 35, 29 30 and the: 


(= Ww ae 


&§ 
— Longitudinal stiffeners .............. 0.47 
— 
strength. 15 OF lmporrance in two respects. it Provigades dlidiySisS Of 

§ 
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ipl PLATE 
‘The last complete. term in the brackets represents the effective thickness of 
the transverse supports. As compared to Eq. 33, it is not a fixed fraction of 


7 . ; the total effective thickness, but varies varies with the loading index and characteris- 74 


tics of the material used. ee 


Typical results for the compression structures designedaccord- __ 
_ing to stability, Eq. 33, and according to yield strength, Eq. 36, are illustrated a 
a in Fig. 9. . In the computations, an idealized stress-strain curve wasused which | 
_ was composed of an elastic portion and a horizontal yield strength limit of © 
itt can be onnerwed, from Fig. 9 th at the lower values of N/E w, “the opti 
provide a structure of 


YELD STRENGTH 
Ss 


ver” 


a, 
whe 20 


FIG, 9. EFFICIENCY OF TRANSVERSELY SUPPORT TED, 

siderably less the yield- difference is pri- 
_ marilya result of thefact that a relatively large weight of transverse structure _ 
is required to support the stiffened panels atthe yield strength, = 
At values of N/Ew equal to 10- the yield strength and optimum designs 
are equal. Beyond this value of N/E w, , the yield- strength design must be this. 
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= types of design informationthat c | 


ses. “This app approach c can be app applied to; of compression on struc- 
tures.12 Not only can the minimum -weight configurations of a particular form — 
types construction can also be investigated. Finally, , such a analyses 2 are 

= value in studying the comparative — of materials of 4 
different elastic moduli and ‘strength 


a 


stiffener spacing, in ‘inches; 


= height of transverse 'Stiffener, in inches; 


s moment | of inertia, in inches4; 
= buckling coefficien 


= length, in inches; 


= subscript denoting transverse stiffener; 


t 
subscript denoting skir skin; 


effective thickness of stiffened panel, in inches; : 


= effective thickness of panel and transverse stiffener, in inches; 
= volume, in cubic inches; 
panel width, in inches; 


subscript denoting 


v = Poisson’ s ratio; 


— pnent; 

— 

2 
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PLATE STRUCTURES 
= radius of gyration, in inches; — =. 

applied stress, in pounds per per square inch; 
column stress, in siete per square inch; 


= ae stress, in pounds per square | inch; 


= stress, in ‘Pounds square incl inch; 


0 average or ‘effective stress, in pounds: per r square inch; 


= crippling ‘Strength, pin gene per square inch; and 


plasticity reduction factor for col umns. ll 


im 
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DYNAMIC ANALYSIS OF ELASTO-PLASTIC STRUCTURES CTURES 


By Glen V. Berg,! F. ASCE and Donald A. diane A.M. nc 


A numerical methodis presented for determining the response sponse of anelasto- 


"plastic structure of a finite of of freedom to by 
a a speed digital computer. — 


In recent years the problem of analyzing the r response of a non- stlitins sys. oo 
tem to dynamic forces has received increasing attention. To the structural - 
engineer the elasto-plastic response problem is perhaps the most Important 
case, and recent contributions in that area have been numerous. 


The present paper looks into a special case of the elasto- -plastic response ~ 
problem, namely, t the multi- ~story elasto- ~plastic structure responding to lateral 
dynamic forces. A complete | and practicable machine method of ‘Solution 
presented. 


‘Note, —Discussion open until September 1, 1960. To extend the closing date one 7 ’ 
_ month, a written request must be filed with the Executive Secretary, ASCE. This paper 2 
is part of the copyrighted Journal of the Engineering Mechanics Division, 
“of the American Society of Civil Engineers, Vol. 86, ‘No. EM 2, April, _ ; 
“i 1 Assoc. Prof. of Civil Engrg., Univ. of Michigan, , Ann Arbor. we ae 
2 Ass’t. Prof. of Civil il Engrg., Univ. of Arizona, Tucson. 
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elasto-plastic constraints. Uniqueness of the elasto-plastic 
qrRODUCTION 
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‘one SOLVING THE EQUATIONS OF MOTION FOR 


The dynamic system considered in this paper is a lumped mass system of 
_n degrees | of freedom, subjected to external forces and/or motion of the base. 7 7 
“Sucha system of order twois shown symbolically in Fig- 1, in which the springs : 


= displacement of of the ‘ith ieee relative toa , coordinate system 


fixed in the base; 
= of the base, relative to a “fixed” frame of ref-_ 


Fy(t) external force applied to the ‘ith pody; 
q OR ix, x,t) = = resistance force exerted on the i-th body | by t the spring and 


= 


and are “displacements relative to ‘base 
FIG, 1.—DAMPED ) DYNAMIC | 
and dots denote differentiation with respect to to time — The equations of mo 


for me — are then 


ck 


1* 


| 
— 
alt 


damping forces or negligible, and therefore the resistance func- 
tions Rj do not involve the velocities x;and (b) the driving forces fj canbe ap- 
proximated by smooth curves through ‘sets of discrete points equally — 


in time. Blast | problems, for example, often meet these requirements. 


[5x + -t) 
2 x;,(t) x;( + 10 X(t) + (t-n)] 
‘By recurs recursive use of these t two formulas one advances “step: by ste step through the 
nag For problems which involve » very ry irregular driving forces, suchas the earth- 
_ quake response problem, it is advantageous to use a single-step method—that . 
is, one which projects to time t +h from the values of xi, X and Xj at t time t, 
and the functions fj, without using the values of the variables at earlier times. 7 


The Runge-Kutta isa step suited to vai 


computers. 
For the equations 


oma 


The the constants m, n, p, and r are 
- 


Princeton Univ. Press, Princeton, 


s,” by . Levy and A. Baggot, Do 


er Publications, Inc., New York, 1950. ee 


ere are many methods of approximating the solution to Eq. 1 numerically, 
and no one method can be claimed to be the best for all cases. 
: 
— ‘ 4 
— 
the formu — 
i 
—— = 
— 
‘ 
q 
— 
— 
— y W. 


€ 


‘This makes m = (q-r) = 0 and thus us reduces the m number of arithmetic opera 


required. A second due to S. Conte and R. 


= 0.856 85614352807 u 


In 


To adapt Eqs. +1 to the to the -Kutta pre one can them! in the form 


to the restrictions that the r number of degrees | of freedom must “be finite fase 
> the he resistance functions Rj must be single- -valued and piece- wise continuous. 7 


‘THE RESISTANCE FORCES FOR AN ELASTO-PLASTIC BENT 


The specific problem considered is the response of a multi- story 
bent to lateral dynamic forces. . The bent isa rectangular plane framework of 


‘and in ine plane of the bent. Connections may be either pinned or fully re- — 

strained. The mass of the structure is assumed d to be | concentrated lat the joints. a 
Damping is assumed to be viscous. Deformations due to axial forces and shear 
forces in the members are neglected, and the effect of axial forces upon the 


5 Kutta Third-Order Procedure for Solving Differential Equations Requiring 
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- sumptions yield a dynamic system having degrees of freedom equal to the num - 
b= 7 ber of stories in the bent. - The principles used can be extended to more general 
problems, such extension is not considered herein. ste: 


with the assumption of viscous damping, the resistance functions 


structure were passing through the equilibrium with unit velocity 


| 
The forces Qi ..nj)are the forces necessary to hold the frame in static 
equilibrium inthe Prive defined by the lateral deflections x;(j = = 1,2,. .n). 
=4 The jth column of the matrix k represents the set of static lateral forces that : 
_ would be required to hold the frame in equilibrium with ‘unit lateral de deflection — 
: at the jth floor and zero lateral deflection at all other floors. A conventional _ 
= analysis is sufficient to evaluate the stiffness matrix. The effect of ec- _ 
centricity of dead load can be taken into ac account in in evaluating | the matrix ‘if de- 
Let the subscripts k and locations in the such that each 
= of the subscript denotes the end of a specific structural member at a 
specific joint;and let My bethe moment at location k, considered positive when | > 
| tends to rotate the end of the member clockwise. Then the end moments in 
i the members of an elastic frame can 1 be expressed in termsof a matrix of in- 
fluence such that 


jt) column of the ‘matrix is he: set of moments that would exist 
‘the members if the frame had unit lateral deflection at the jth floor and zero 
lateral deflection at all other floors. The matrix » can be evaluated by con- : 
If the lateral deflections become sufficiently large, the elastic limit may b a 
exceeded at one or more locations in the frame. This will, of 
the restoring forces and moments 


Be 4 _ restrained, and that all members in the frame, columns and girders | alike, wad - 
_ ideal elasto-plastic moment-rotation characteristics. The typical moment- _ 
rotation curve follows a linear elastic branch until the moment reaches the 

- plastic hinge moment - Upon further deformation, the moment remains constant — 
a plastic hinge forms, making a in the reversal of 


i. Be this paper it is assumed that the connections a are either pinned or wd 


EM2 b ELASTO-PLASTIC STRUCTURES > 
ae 
— 
_ locities, and restoring forces Qj. The dampingforces are # 
— 
_ The restoring forces for an elastic frame can be expressed in terms of a 
| matrix of stiffness coefficients k suchthat = 
| 
4 
‘ 


the plastic hinge saoment, either positive or negative, is then 
follows a plastic branch as before. reversals cause the moment - 
rotation curve to follow a “nee 
the frame behaves just as 

~ though the member were hinged at that location, with a constant moment applied 
— to the hinge. This condition prevails until the next hinge forms or until a de- 
crease in strain occurs. at an existing plastic hinge location, causing the hinge — 

i _ E. Cohen, L. S. Levy and L. E. Smollen developed§ a a procedure | for adapt- 

ing the method of normal modes to an elasto-plastic frame with infinitely rigid 

a girders, and Leo Schenker, M. ASCE,? indicated a method of extending this to 


ie include the frame with flexible girders. In both procedures, an elastic solution 


‘FIG. 2 IDEALIZED MOMENT - ‘ROTATION DIAGRAM 


-_—- at some location. The frame is then modified by inserting a real = doll 
at the plastic hinge location, the new normal modes and frequencies are com- 

Be cone and the solution is continued as a superposition of the conditions at the | - 

_ instant of transition and the changes found by solving the modified frame. Each > 
time a new hinge forms or anold one disappears, the frame is modified accord- 
ingly. In this way a continuous solution can be 
For a step-by-step numerical solution, an analogous and somewhat simpler — 

mee can be devised. The conditions at the end of a time step can be — 
found as a of the Conditions at the of the step and the 


7 nf 6 “Impulsive Motion on of Elasto- -Plastic | Shear Buildings,” by E. Cohen, L. s. ‘Levy and : 
E. Smollen. Transactions, ASCE, Vol.122,1957.0 
.  7*The Dynamic Analysis of Tall Structures in the Elastic and Inelastic Ranges,” by 

Leo 


‘strain after reaching plastichy, follows a path paral- 
— 
- 
if 
4 — 
the “elastic i ements” -that ic the changes which would he caused 
x 
— 


ELASTO- PLASTIC STRUCTURES 


by the motion the time step if plastic hinge rotations did not occur. 
_ resulting moments may exceed the raatie hinge moments at one or more lo- A 
cations in the frame. To remedy this, ‘ “corrector” solutions are superimposed 
_ in which the frame has hinges at these locations, and these hinges ; are rotated — 
: in such a way that the total moments, obtained from superposition of the elas- 
_ tic and corrector solutions, nowhere exceed the plastic hinge moments. ‘The 
superposition of elastic and corrector solutions is illustrated in Fig. ae 


1 


7 


FIG. 3. _-ELASTO- PLASTIC SOLUTION BY SUPERPOSITION. _ 


on was s perhaps the first to apply the concept to the 
e _ type of problem considered here,9 and H. H. Bleich recently formulated a 


8 “Theory of Perfectly-Plastic Solids, by W Ww. ‘Prager ai and P Hodge, John Wiley 


‘Sons, 
Response of Beams in the Elastic and Plastic "by W 


— 
— a 
or 

— — 
— 
— 
— | 


7 method of analysis on the same —_— in a somewhat broader prob- 
lem.10 However, no method of solving the resulting inequalities was given in 
the Bleich paper, and the existence and uniqueness of the solution were as-_ _ 


yy The effects of plastic hinge rotations upon the forces : at the floors and the 
‘ bending moments in the frame members can be expressed in terms of two ma- 
trices of influence coefficients b and v, evaluated as described ae 
ret ds to clockwise rotation of the ‘end of the member. (The hinge 1 ro = 
tation in Fig. 3 is positive.) Take the bent in its initial (unstrained) position, © 
insert a hinge at location £, and rotate the hinge through an angle @ ¢ = 1, with — 
. zero lateral deflection at all joints. _ Evaluate the lateral forces at all floors 
and the end moments for the members in the frame - The lateral force at the 7 
_- floor is the influence coefficient bj9, and the moment at location k is the © 
influence coefficient | Both matrices b and V can be evaluated by conven-_ 
tional static analysis. It will be seen that only those locations at which plastic _ 
hinges form are significant tothe dynamic response. 
If the lateral deflections x and the plastic hinge rotations @ were known for — 
_ some instant of time, one could evaluate the restoring forces and moments as 


one must also find the plastic hinge rotations which satisfy all the conetreiote. 
f of the elasto-plastic system. Aniterative methodof solution is developed below. 
a ae _ Consider the time rates of change of the forces Qj, the moments Mx, and the | 


plastic hinge rotations ¢,. Let the plastic hinge moment at location k be Px 


The moment M cannot exceed the plastic hinge moment P in magnitude. — 


if a plastic hinge exists, the moment must be either constant a 


_ consistent with the moment. If the moment is decreasing, the hinge ciel 
3. Equivalent conditions | exist for a negative plastic hinge. 

- 4. If the moment is less than the plastic hinge moment in magnitude, a 


‘These. may be expressed mathematically as the following set of c 


eae 10 “Response of Elastoplastic Structures to Transient Loads,” by H.H. Bleich, Trans- Teams 
a 


ctions, New York Academy of Sciences, Ser. Vol. 2, 1955, PP. 
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3. = then > 


and given any x’s whatever, and any M’s not “exceeding the plastic hinge 
_ ments, * there exists one and only one set of ’ s and and ¢’s s which 1 satisfies Eq. 18 


and Constraints 17. The proof follows 
First, observe that if for one or more values k, Mx < Ps then 
hee 0 and there is no constraint on Mk. _ Thus the solution for $, is unique, it~ 
_ does not affect other locations, and tl the effect of conditions at other ‘locations — 


im upon My is immaterial. Hence it i is sufficient to consider only those locations | 


the Kronecker 


* The ae in which all members at a joint become hinges simultaneously is excluded. Da 


eness theorem can be established: 
if 
‘Se 
— 


e vector ais unrestricted . Toe establish \ 


: scot a, there exists one and only one » pair of vectors y and z that satisfies | 
; Eq. 20 and Constraints 21. . This amounts to a partition problem in n-dimensional | 
_ Consider the columns of the matrices -6 an and E as vectors mM, NQ +++ n> 
E2,-++-tn- That is, let {nj} = {-5; jit and = i} Then let 


be a set of vectors such that every ‘a is either ‘ni or There are 29 


. sets. If the n ’s and §’s partition the space, every vector a in the | space can be 


set of a;’s for which 
where every qj is nonnegative. If this is s so, then on those i for which aj = 7j, 
one has yj = qj and zj = 0; and for those i for which aj = &j, one has yj = 0 and an 
Thus the vectors y and z corresponding to any given vector a are unique 
To illustrate this concept, consider the two-dimensional case. Suppose posi- 
tive plastic hinges exist at locations 1 and 2, and all other locations are | elastic. 


7 -16,580 105,560} + | 
i‘ — columns ot of the. ‘matrices in this e equation are plotted as the e vectors C1, oa, 
ag , and 19 in Fig. 4. The four vectors partition the plane so that any — 

in the entire plane can be formed by linear combination of one and only one 

set of aj’s with nonnegative coefficients. 

_ By virtue of a recent theorem due to H. Samelson, , R. M. Thrall and O.— 

-— necessary and sufficient conditions for the columns of & and -6 to 


- Partition the space are that every principal minor of the matrix E be positive. 


coefficients.” In other words, f any given vector a there is is one and one 


5 "equilibrium configuration. This is the frame used f for the “corrector” ee, 


a hr: If any coefficient is zero, say qq, = 0, obviously one can have either 4,, 

a, Se _™y- However, in the case ym = Zm = 0 and uniqueness is preserved, — 
me oy “A Partition Theorem for Euclidean n-Space,” by H. Samelson, R. M. Thrall and 
O. Wesler Proceedings, AMS, Vol. 9, No. 5, October, 1958, pp. 805- 


iniqueness, it must be shown 
— 
“4 
— 
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— 
vw bee 
— 
— 
— 


4 4, 4.—PARTITIONING VECTORS, 
y above. Let the hinged joints in this frame undergo hinge rotations ¢, with later- — 
al deflection 1 of the Goors 1 eee The resulting strain energy in the frame 


(23) 


OK My . 


— 
— ||| | i 
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Because Eq. 25isa strain-e ~energy ‘equation, itis a a positive definite. quadratic 
form. The principal minors of y are therefore positive. Finally, since £ was a 
obtained from v simply by changing the signs of the off-diagonal elements of 
certain rows and the corresponding columns (see Eqs. 19), the principal minors — 
 - ¢ are positive. It follows from the Samelson- Thrall-Wesler theorem thata _ 
ae The solution can be found by the 1 following iterative scheme: Let zP be th 
ae ae! " approximation to the vector z , with all elements of zP nonnegative. (Zero - 


is a convenient and first approximation. ) For in no- 


tation, 


4 
0, 


UR , let 
Repeat this for k = 


The to the vectors y and z is then complete. It satisfies 
Eq. 20 and the constraint z_ = 0, but may | fail to satisfy the other constraints — ; 


in Constraints 21. _ The process is repeated until the amount by which the ap 
~ proximation fails to satisfy all the constraints is insignificant ‘a aoe 


the vector z is determined (and therefore e s) the rates of 


4 
= 


iy 


“iterative method of solving linear Eq. 20 and Constraints 21 reduce 
finally to a system of linear equations, since at least one of Zk, Yk Must be zero” 
for each value of k. The nonzero elements form a system of n (or 
equations in the same number of unknowns. In this application, the Gauss - Seidel — 

Aa procedure is s used to set 1 up the equations as well as to solve them. The pro- 

cedure lends itself well to machine computation, and it turns out that converg- 

ES ‘Using this procedure, at each step one must perform three matrix multip 

_ cations and solve a system of equations equal i in number | to the number of plas- 


4 
— 
— 
| 
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<- In developing the above procedure, rates of hens of deflection, hinge ro- 

tation, and moment have been used. As soon as one introduces finite time 
tervals in place | of differentials, one loses the claim to uniqueness, since it _ 
‘cannot be claimed that reversals in strain do not occur within the time inter- 


val. However, by taking the intervals sufficiently small and considering the 
L strain rates ase everywhere monotonic within } any on one | time interval, the exact 


4 


s+1 


be 


Then the finite - difference equations are vd 


which correspond to Eqs. 29 and jar’, 


‘then either am, 


then either sand A 

< AM, < Px Miss and d =0 

foe 

ADAPTATION FOR THE ‘COMPUTER 


In using the computer, both storage requirements and ‘computation time can 
be reduced by expressing the moments in nondimensional form. — _ To do this, 


. 
— i 
4 Stem can be approximated to any desires 
a) degree of accuracy y tem the subscripts s and s+l denote 
— condition at the end of the sth and Aine, 
the value of the variable xj at the end of the 
— 
— — 
— 
Py - Myg and Ad <0 83) 
or - My, < AM, < Py 4 
> .34) 
— 


vy Py 


do not the of the solution. In this 1 no- 


> wes 


and ‘the elasto ) -plastic. constraints (33) and become 


k,s 


lation of the problem. Let 
of moment locations order a the matrix o) 


> te 

rz = pth 


4 
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BLASTO- PLASTIC STRUCTURES 


‘Then, starting with = p=1,compute 


+ 


- venguene test is satisfied, one accepts the last aporoximation as the correct 


Solution f for the and AQi and ‘Ap, from Eas. 36 and 
a 


NUMERICAL EXAMPLE 


the Mowing mum bent shown in Fig. 5. fy 36, 000 one “one for 
= 
= Pp = Py = 13 1308 in.- in.-kips 
= 2072 in. .-kips 
- Because the bent is ‘symmetrical and the effects of axial forces on the member 
_ stiffnesses are ignored, the stress distribution will also be symmetrical. Hence — 
only half the moments need be recorded. Moreover, the moment at location 5 — 
cannot reach plasticity and therefore does not need to be recorded. ia 


plastic hinge mom 


a 

The matrices of influence coefficients kij; Dig Ukj, and V_ ¢which describe 


elasto- system are the forces and moments showe 


— 
— 
— 
sThen 
— 
— 
which is an indicator Of the amount Dy which the solution falls to 
Constraints 38 and 39. If eP is larger than some predetermined errortest 5. 
. 
— 
— 
FF if 


EXAMPLE, 


dead | load on the moments 


-24.11 20.20 | 


[as 


~1584 1021 ian 


“in, -kipe/in. 


69,700 25,720 -9,140 
25,720 88,300 -31,370 -56 


,930 
| 2 40 -31,370 80,010 -48 
—L-16 580-56 930 -48,640 105, 


Converting the moments to nondimensional form according to Eqs. 35, one ob- i. 


— 8 — 
— 
— i — 
— 
— 
in Figs. 6 and 7. The effect of eccentrici yo 
— 
— 
— 
— iz 


ELAS'.0-PLASTIC STRUCTURES 


29.12 31.97 29. 25 


-.5987  .2012 | 

-.0833 


L 0000 -.3921 
23553. 0000 
L-.1605 4102 
t at time to the system is in the cor contighvetiog 


Py 
= 10.00 in. 


and that plasticity y has not yet been encountered. at | any lo location. 


zy. = fo = = 10.00 kips, 


and that these forces dec decay linearly to zero at time to + 0. 1 sec. The complete 


_ process for ‘determining inelastic behavior and executing | one time ‘step 0 
_ The Runge-Kutta constants for this example are those given in Eq. 9, and — 
time interval of h = = 0. .02 sec is used. These initial conditions and time interval 
: have been chosen for illustrative purposes, and are not necessarily typical or 
_ Because plasticity has n not yet been reached, the restoring forces at time 
to are given by Eq. 13, ak 7 


moments at time to are 


0.5194 -0.4740 

_| 0.7222 -0.5664 

=| -0.5987 0.2012 
-0. 0833. 


— 
— 
— 
— 
— 20.20) L4.900J ~ L 36.29) “PS 


qs. 11 and the initial conditions, 
= 10.00 in./sec 


+ 25. 24) = 4453.4 in./sec? 
10 - 36.29) = 4 /sec? 
= (0.02)(10.0 200 in. 


(to) = (0. .02)(25. 00) = 

ve 


ijand 


14 
— 
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— 
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ELASTO- PLASTIC ST STRUCTURES Le 
The superscripts x and Zz identity the variable associated with x. - Before pro- 


+= K 
to the second of pone must the functions 3 “io” 


which, ‘in require the of the As. The first terms in 40 


0. 7222-0. ‘fo. 
-0. 0833. 0.2464 ‘orn 0. 

The iteration of Eqs. . 40 ‘and 41 yields, for init 1, 


“0.0444 - 


0064 + (-0.1311)(0. 0166) 


uy 0.570" - 0. 0086 = -0. 5793 


= 


— 
fe | 
— 
— 
— q + = -0.9722 --1016 
—_— 
el =) 40, | =0.0402 


1880 


PE a 
— 
the Process | @ 7 


_ ELASTO- PLASTIC STRUCTURES 
One can now return to the second of Eqs. 6 and compute 


10.784 in. 
—— going to the last of Eqs. 6 one must evaluate foci 
Eqs. 36, 37, 40, 41, and a again, just as illustrated earlier, 


Ay 


= = 0.00 
The last of Eqs. 6 is now executed, yielding ee 


= -10. 856 in./sec. 
@ Eq. 1 now yields the vaeeeeeres and — for the end of the time 


— — 
| 7 = 
Ap, = 1.000 
— 
| 4 
— 
— 
1 
— 


The terminal values of Q and ip tor this. time step must now be computed 
, 41, and 42, as before. es 


- 


h) = -23.83 ki 
= 


One now advances to the next time step, taking these terminal values as the 


initial values for the next 


he multi- -degree elasto-plastic response problem has been programmed 
s ie solution on an IBM Type 650 computer at the University of Michigan. Both 

_ the Milne Predictor-Corrector Method and the Runge-Kutta Method have been 

- used. The programs were written for the basic IBM 650, without index — - : 

_ tre, automatic address modification, or floating point hardware. If the pro- 

grams were to be rewritten today, they could be simplified a great deal bytak- 
The 2,000-word memory was able to accommodate a system with as many 
3 as 16 degrees of freedom and 16 plastic hinges using the Milne method, and | 


the Runge-Kutta method the system could go to 12 degrees of freedom and 

_ 16 plastichinges. The more severe limitation for Runge-Kutta occurred partly 

_ because viscous damping and a form of inelastic deformation n not considered in 

this paper were taken into account in that program. 

_ The four matrices, kjj, Big, and were stored in the memory initi- 
ally as part of the data for the main program. A separate program was written 

‘ to evaluate all of t these matrices. Input for the matrix evaluation program | 

ie comprised the number of ‘stories, number of bays, story weights, ain call i. 
stiffnesses, plastic hinge locations, and plastic hinge moments. Acting on this a 

— i information, the machine wrote the slope deflection equations, solved them by é, 

iteration, evaluated the matrix elements, and printed the matrices on loadcards 
ready for use with the main program. 
the matrices A_j and one need consider only those locations 

a at which plastic hinges form during the response. There is no need to il 

% late or record the moments at locations that remain elastic. Thus the limi- 

; ; tation of 16 plastic hinges is not as severe as it might seem. Intuition andex- _ 
perience with smaller systems may enable one to predict which locations are “ie 
most likely to become plastic during the ‘Tesponse. The unlikely locations can 
_ be omitted from the matrices. One can process the output data after the run | 

5 _is completeto check whether or not all of the unrecorded locations did, in fact, 


— 
+h) = 1.0000 — 
— 
— 


The programs were applied to the problems of response to 

and earthquake. . In each case the driving forces were put in punched cards 

. which were read by the ‘computer as the solution progressed. ' For the blast - 

problem, the equations governing the dynamic pressures were given to the ma- , 

chine in a separate program and the machine calculated the driving forces and 

- converted them to punched cards ready for input for the main program. hy For 
the earthquake problem, recorded accelerograms were approximated by piece- — 

wise linear functions, and the time-acceleration coordinates of the yy cool 


chine accepted this input data and interpolated where necessary to tind the de- 


See with both the blast and earthquake problems chews ts that the 


equipment that exists today. Tomorrow’ s equipment, or even a a late version 
and permits many of the 


Michigan Research Institute under contract with Air Force Special Weapons 
Center, Albuquerque, N. M., and in part upon research conducted by the = 
_ author for his doctoral dissertation at the University of Michigan under a doc- 
toral committee of which Bruce G. Johnston was chairman. Computation w was 7 

_ performed at the Statistical and Computation Laboratory of the University of af 
Michigan. The authors wish especially to thank Bruce Arden and Mrs. Sarah 2 


Brando, the for their unlimited cooperation. 


are | here for reference. 


damping coefficient; 
force 


ient; 


= time in numerical 


= subscripts associated with mass or force; = ; 
= subscripts associated with bending ‘moment or ‘plastic 


= stiffness coefficient; = 


— 
im 
— 
— 
4 
Seal ear in the text. Those used fre- “ey ee 
— 
— — 
— 
m,n, p, q, r = constants in Runge-Kutta method of int 


_ = plastic hinge moment; 

force; 


= rotation - moment coefficient; | 


dimensionless bending and 


= displacement - moment coefficient (dimensionless mo- 


= rotation - moment ‘coefficient (dimensionless moments); 
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BEHAVIOR OF BUCKLED RECTANGULAR PLATES? 


A review is presented of recent work done by the writer on the post- 


‘Unlike simple columns, rectangular plates which are supported | on all 
oY edges may |carry considerable load beyond their buckling load. Under some | 

it may be advantageous to utilize this additional load- 
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— of a separate analy 
plates in longitudinal compression. Results ted. Adetailed 
P ttern are presented and interpreted. 
~nomenon of change in buckle in sheay of a plate buckled incompres- 
q 


re: In discussing the post- aiiten behavior of plates the basic equations <a 
_ be presented and the method of solution will be indicated. | _ For the problem of 
simply supported rectangular plate subject to longitudinal 
theoretical load- curves will be presented which indicate possible 
: changes i in buckle pattern. - Results will be presented of a separate study of © 
change’ in buckle pattern. Comparisons will then be drawn between the results 
of the plate problem and other theoretical results and with experiment. Inas- — 
- much as the details of these analyses | are available elsewhere, 2,3 ' they will ll be 
; omitted herein. Finally, a a discussion is presented of the stiffness in shear of 
an infinitely long, simply supported plate already buckled by longitudinal 
compressive — The details | of this er are presented in the Appendi az 
B EQUATIONS AND METHOD OF SOLUTION 
os © von Karman has presented the basic nonlinear differential equations for 
Sp plate undergoing inegee deflections. is. For 2 a plate with no lateral load the aie 


oN 


where Ny, x» Nos N,, are the resultant forces, w is the lateral deflection, Dde- 
is Young’s modulus, h repre-| 
notes the plate  stiftness D= 12(1 - pre- 


tate x and | y y when they appear as a subscript following a comma ‘indicate par- 


_ tial differentiation of the main symbol with respect to x and y. Only — a 
initial eccentricities subject to in- plane loading are considered. ; 


+Wy,Wy, 


= bs 2 “Loads : and Deformations | of Buckled Rectangular Plates,” by Manuel Stein, NASA 7 


‘Stein, NASA Technical Report, R39, 1958. 


= 


2 


| 
| 
| oa 
As 

| 
proper substitution in Hooke’s law for a plate the forces can be expressed 
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BUCKLED RECTANGULAR PLATES 


= where u v are displacements in ‘the xa and y , respectively. 
«Eqs. 1 ahd 2 together with a complete set of boundary equations determine _ ; 
‘Notice that, if the N’s from Eq. 2 are substituted into Eqs. 


a a perturbation method the nonlinear equations are converted into a 
‘Set of linear equations by assuming that u, V, and w may be expanded in a 


‘power series” in terms of an arbitrary parameter « e. For the present purposes 
u, Vv, and w are to be expanded about the point of a eran (at buckling € = fan a" 


if Eqs. 2 2\ and 3 are substituted into Eq. 1, - equating power of € leads to an 
infinite set of linear equations. The ones corresponding to the zeroeth and 
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series expansions of these basic equations, except for one study15 that made = 
— 
> 
at 
a" 


to the post-buckling behavior of plates. Solutions of additional 
equations give a second approximation and then higher approximations. Ks The 


advantages of the present method are greater accuracy and 
tical Solution. post post- buckling behavior of of ‘simply support 
r plates in end compression (the edges are constrained to remain : 
Be Pray is investigated in detail by solving the first few of the equations. 
7 Nondimensional load- shortening curves presenting some of ‘the results of this | 


problem are shown for plates of various length-width ratio in Figs. 1 to i. 
_ For plates of finite length- to-width ratio the curves are also labeled accord-— 


were obtained by using the values of m eteieeientian to the lowest buckling 
load. In addition, -load-shortening curves are given for other values 
intersect with these basic curves for the range plotted. > 
7 Ly, Change in Buckle Pattern.—The intersections of the load- shortening curves 
— _ possible changes in buckle pattern. . In order to discuss the phe- ; 
nomenon of change in buckle pattern on a sound theoretical basis, a simple =e 
model was chosen, as shown in Fig. 6, which has the of 


been analyzed? i 
torsional springs and supported laterally ‘by ‘nonlinear extensional 
springs (the nonlinear springs represent plate membrane stiffness). _ Each of 
three elements are of length l, the spring constants are C for the linear tor- 
_ sional: springs, and K and K, for the nonlinear extensional springs, and this 
column is subject either to a load or to a shortening. Pe See 
a.  * shown inthe Fig. 7, this idealized structure duplicates the important 
ie property of plates in that there is an intersection of the load-shortening curve 
* one mode (say the symmetric mode which corresponds to the lowe 
buckling | load) with the load- shortening curve of another mode (antisym— 
metric). ‘The 
a buckle pattern which is neither ‘symmetric n nor antisymmetric) from the — 
metric to the antisymmetric buckling configuration. Change in buckle pattern 
= to occur at a loading corresponding to the intersection of the load- a . 
ill shortening curve for the symmetric buckle pattern and the transition curve. 
Se The load at which the first change in buckle pattern starts to occur is called 
j % the secondary buckling load. Whether the transition curve is stable is shown 
= =—s-_: to depend on the stiffness of the nonlinear springs and upon the method “a 
loading (controlled load, controlled shortening) which is shown by comparing 
is Figs. 7 and 8. — " However, the load (and shortening) at which the buckle pattern a 
. 2 starts to change (secondary buckling) is independent of the method of loading. 
Secondary buckling always occurred fi for loads and shortenings greater than — 7 
i those given by the intersection of the load-shortenin 1g curves for the sym- 
metric and antisymmetric equilibrium configurations. 
‘Stability occurs as one might expect intuitively. For r example, for con a 
i trolled load, if the transition path from the curve for symmetric buckling to 
that 
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path is stable; otherwise, the transition path is unstable. Ino order to enee- 
the stability of equilibrium position for the column problem in a 
vigorous fashion it was necessary to examine the second variation of the total - 
potential energy. It would be expected that for plate problems, stability could 
% = = } == | 
joe = Changes in buckle patterns are not calculated for plates of finite length- -to- 
_ width ratio for the present analysis because of the extensive calculations re- 
quired. However, preliminary calculations have indicated that for length-to- 
_ width ratios near unity the change in buckle pattern would be rather abrupt; _ 
a whereas, for higher length-width ratios, changes in buckle pattern would be 
continuous or at least less abrupt. 
. For an infinitely long plate (Fig. 5) the number of buckles along the length 7 
be. expected to change continuously as the loading progresses. The 
a buckle length for a given shortening would be such that the load is a minimum. : 
- 3 Note that the results giving the lowest load for a given shortening for length-— 
_ width ratios 2 and 4 do not differ much from the infinite-plate results. Indeed, 
_ the infinite- -plate | curves form an envelope for the finite-plate curves. 
_ Comparisons with other theory. —For the square plate buckling tuto 
‘square buckle both Samuel Levy10_ and S. A. Alexeev15 obtained essentially 
exact solutions. — For a square plate buckling into two buckles (m = 2) only 
Alexeev obtained an exact solution. . The present ‘results for the second ap- 
‘proximation agree with the results of Levy and Alexeev. _ te 
_ For plates of various other finite length-to- width ratios, previous results 
. are available for the initial slope after buckling. As can be seen from Figs. 1 
4+ , for ‘some length- to- width ratios used, straight-line load-shortening 
curves based on these initial slopes would give unduly higher loads for given _ 
shortenings everywhere in the post- -buckling range except immediately after 
__, The best available previous results var the infinitely long plates are those — 
of W. T. Koiter.11 _ As shown in Fig. 
and solution of | Kar] 


range Marguerre’: solution continues to the first approximation where- 
as Koiter’s solution deviates slightly from the first approximation in the di- ” ‘ 
rection of the second approximation as shown by the dotted line in Fig. 9. The 
_ ‘Tesults of the second approximation give lower, more accurate loads than = 
e __ Experimental Results and Comparisons of Theory With Experiment.- to 
obtain simply supported loaded edges is impractical in laboratory experi- 
ments. The experimental results \.ere therefore obtained for panels = 
: to af ae loading which results in almost complete clamping of the loaded 
a edges. However, if the panel tested is long compared with its width (say of | 
a -to-width ratio 4 or greater), the size and shape of the buckles near the 
_- center are almost unaffected by this clamping. The experimental results — 
= j. which are compared with theory are for panels that have a length-to-width | 
a ratio of at least 4. Hence at least certain of the experimental values ined 
_ ghould be directly comparable to the simply supported theoretical results. — 
stated previously, the theoretical results for length-to-width ratios 2 and .. 
= not very different trom the results for the infinite plate. Thus the experi- = 
‘mental results may be compared to the theoretical results” for the ‘infinite 
plate. Such comparisons are shown in some of the following figures. Data are _ 


} 


— 
q 
he 

| 


wl 

Ist APPROXIMATION 

— 2nd APPROXIMATION 


> 
FIG, 9.—COMPARISONS OF THEORETICAL LOAD- wig. 
_ SHORTENING CURVES FOR AN INFINITELY 


PLASTIC, 


THEORY 


10, . COMPARISONS OF LOAD-SHORTENING CU ES AS GIVEN BY (ELASTIC 


— 


from th three different specimens. Inal all of speci- 
Re = the loaded edges were ground flat and perpendicular to the longitudinal 

axis of the | specimen. They were compressed “ flat ended” between the platens a ; 
of the 1,200,000-lb capacity hydraulic testing machine at the Langley struc- _ 


tures research laboratory, which applies load through the use of a — 


Plate Supported by Multiple-Bay —In the first test discussed, an 
attempt was made, by using a fixture having lubricated knife edges free to ro- | 

tate, to provide the edge conditions usually specified by theory along the —— 

loaded edges of each panel: fs simply supported straight edges free of in-plane > 


A comparison with theory is presented in Fig. 10 of the experimental load- 

_ shortening curve for this test. The experimental curve shows abrupt changes 
- corresponding toabrupt changes in buckle pattern from 5 to 6 to 7 to 8 buckles © 
while the heoretical curve, is based on continuous in 


Also shown ii is the laa — Pb 2 Gacy rae gages at the crest of a 
- buckle indicated the plate material had been strained into the plastic range. 
_ The type of changes of buckle pattern obtained here with a hydraulic - type | 
testing machine is similar to that described ina previous section for a con- — 
trolled-shortening type of loading. In consideration of the practical difficul- 


ties of measuring total shortening, such as how to account for the bending of a 


— 


testing machine platens, ‘the present ‘agreement between experiment 


a _ Bending strains at the crest of a buckle for the same test are plotted 
|. Again 
ie abrupt changes | in the experimental results do not appear in n the theoretical 
results. The agreement between theory and experiment is good. 


a -Stiffened Panels.—Data on four Z-stiffened panels similar to those used 


-aircraft-wing construction were obtained in the range from zero § strain up 
to several times the buckling strain. These Z-stiffened panels were part of a _ 
Be In Fig. 12 the experimentally measured ‘strains at the crest of the buckle © 
are plotted against stiffener strain and compared to the present theory. From > 
= the comparison of the results shown, it is evident that, although there is scat- _ 
ter at buckling, the theory for simply supported plates gives strains at the © 


crest buckle that agree with practical experiment in the post- 


— Hat- Section Stiffened Panels. - —In Fig. 13 buckle depths measured ieeen a 
ae series of tests at Langley on panels with hat- -section stiffeners are ered a > 


- against stiffener strain and are compared to theory. There isa considerable — 


scatter in experimental results. However, it may be stated that in the — 
buckling range, theoretical sebiiie for the depth of buckle of simply supported 
with experimental results on such practical stiffened 
SHEAR STIFFNESS OF A BUCKLED PLATE 
‘An extension of the compression problem is presented which treats the 


— of the stiffness in shear of a rectangular plate that has been buckled . 


16 “Compressive Strength of Flat Panels Z and Hat-Section Stiffeners,” by 
N. Kotanchik, Robert A. George | Ww. Lender and John Neff, Jr. 
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FIG, 14,.—COMPARISONS OF NONDIMENSIONAL SHEAR 
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Bruno A. Boley who presented!7 estimated curves based mainly on theoretical 


If a plate is unbuckled its shear can be identified immediately 


shear modulus However, if a plate is itis ina 


_ rather complex state of stress, and { it is necessary to define what is meant by | 


have already buckled locally. Previous work on this soutien has been done by 7 


stiffness. it seems reasonable | that ‘the ‘Shear stiffness of a buckled 


shear stiffness obtained in this manner, , and this analysis is further 


of length-to-width ratio 4 or greater. 
_ Having already obtained a solution for the displacement and eniee for s a 
iz rectangular plate buckled by longitudinal compressive loads the method of so- a 
i lution is to go back to the original von Karman equations (Eqs. 1) and to sup- 
i pose an initial state of stress and deflection given by the compression solu- 7 
tion and to require an additional state of stress and deflection corresponding oo 
to a small uniform in-plane shear loading at the edges. Since the additional 
shear load isconsidered small (infinitesimal) the equations may be linearized. a 
F ; Solutions corresponding to the first and second approximation of the « com- 
~ pression problem for an infinitely long plate are shown in Fig. 14. These so- _ 
} lutions show that at about 3 times the buckling compressive load the = 
stiffness of the infinite panel decreases to about half its shear stiffness (G) at 
_ buckling, but that above 3 times the buckling load the shear stiffness begins to a 
increase. _ The present results differ with the recommended ¢ curve of 7 -. = 
‘shown in Fig. at least for lower loads. 
OTHER PROBLEMS SOLVED | 
In addition to the e eompression | problem solved, ite is s available2 the so- so- 
_ lution by the writer of post-buckling problems of rectangular plates subjected 
ig : to auniform temperature rise with three sets of boundary conditions restrict- a ; 
ing in- plane displacements. results of these will “not dis-— 


a. A linear set of equations has been derived to replace the nonlinear large- 
deflection equations for plates. The advantage of this set of equations is that 
Bre a way to obtain accurate results in simple form. The linear equa- 
tions are subject to the limitation that they cannot be used to solve post-_ 
4 buckling problems for plates with initial eccentricities. erties) 


tae 17 “The Shearing Rigidity « of Buckled Sheet Panels,” by Bruno A. Boley, Journal of | 


longitudinal compressive loads. This problem is important inthe stress 
ee. 
ff — 
4 
sulting average unit shearing deformation except that, since the probiem is 
nonlinear, the shear stiffness would depend on the magnitude of the applied 
_ ‘shear stress. It was thought that a practical and useful measure of the shear 
stiffness might be obtained by restricting consideration to an applied shear 
\ 
ad 
q 
| 
| 
\ 
‘ 
— 
x 
— 


or the: compression the second of the 
theory agrees with exact results for the square plate. Results for plates with 
finite, as well as infinite, length-to-width ratio indicate that the effects of 
change in buckle pattern must be considered. For an infinite plate, results ob-— 


tained in first approximation agree with the best previous results for much of 
the range, but results for the second approximation give | lower and more ac 7 
‘The: comparisons made here indicate that, for extreme fiber strains 
‘ ; deflections at the crest of a buckle, the present theoretical results for simply 
Plates with straight edges" free of shear agree well 


sults indicate that intersections between « curves for the 1 various och 
: configurations lead to change in buckle pattern. The load at which change in © 
buckle pattern occurs is shown to be independent of the type of loading, but 
‘the manner of change does depend on the type of loading. Further ‘computa- 
_ tions similar to those done for the simple model are necessary to give pre- 
; cisely the secondary buckling load and the transition loading paths for finite 
| Seta simply supported plates : subject to the various types of — 
nal compressive loadings (controlled load, controlled shortening, etc.) ‘=o Se 
The shear stiffness of an infinite plate buckled by longitudinal compressive ~ 
load has been investigated in detail indicating that the effective shear stiffness _ 
a _ decreases” to almost half the shear modulus at 3 times the buckling load, but | 
increases again to about three- wea” the an modulus at 5 times the 


5 APPENDIX. _IN- PLANE SHEAR OF AN 
PLATE BUCKLED IN LONGITUDINAL COMPRESSION 
An infinitely long plate is considered to be in nan initial state of oten an 
deformation corresponding toa given compression load. Such a state of 
stress and deformation is given by the solution of the compression problem. 
‘The problem is. to find the shear stiffness from ‘the. additional stresses and 
J deformations corresponding to a small, uniformly applied shearing force a 
Pe the edge. By considering the shearing force n to be a small perturbation, the 
7 equations for the additional stresses and defor mation may be made linear by 7 
‘neglecting squares and higher p powers of the perturbations. The linear equa- 
tions obtained in this manner are solved by the Galerkin method. — 


7 lutions be u 


— 
a # ‘The results of an analysis of a simple model have been presented, indi- _ oe 
a 
oy 
have the initial solution ug, vo, Wo of Eqs. 1 and let the desired additional 
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xy = Nyy + 
where the YN s with zero eanoertgte are given by Eqs. 2 if u, v, w in Eqs. 2 are 
replaced by ug, vo, The Ny’ s are found to be terme n?): 


wo, x + | (3 %0,y 1, 


lowing relations to be satisfied, which were formed 
5 ite 1 (again terms): 


with the initial values of the stresses and displacements 


and the e boundary conditions determine the problem. 


Conditions .—The 2 boundary conditions considered are the 
as those considered for the compression problem except for the provision of : 


“the uniform shearing force present at all The of courciante 


zing 
x,0) = w(x, b) = 


—— 
— 
— 
— 
q 
1 
— 
A 


where a/m is the buckle length. hie, the deformations are considered to be 


periodic inthe x direction, 
Solution. Eas. 6 the and second of 7 be- 


- 


Vv 


e 


z x denis the constant A is to be determined by the Galerkin method. Now substi- 7 
- _ tute in the foregoing equations for Wo and w 4 and solve for u,; and v . (Sub- 
Ti stitute for thezero ‘subscript terms in these squations corresponding functions © 
from the writer’s previous work2 according to the degree of approximation — > 


For the assumption (Eq. 15) the Galerkin method gives the following equa 4 7 


fr from which A is 


[ w, - (Ny 


dy = 


— 
— BUCKLED RECTANGULAR PLATES 
x,0) 
iy 
— 
— 
Wi,xx * + 2 Nxyo¥1,xy +! — 


compression ‘problem and using the assumption of Eq. 15 for the additional 


= 


is the wave m b/a and k / 
The average unit shearing deformation, I, is given 


where @ is the average value u along the length 


The stiffness is defined a as nM } which an 
is equal | to the shear modulus G ==> 
(Eq. me the following is obtained 


4+ | u) B 


The curve labeled first approximation in Fig. 17 was obtained by using in Eq. 
Es 20 corresponding values of k and f given in Table 1 of the writer’ 3 previous 
No details are presented for the second approximation curve of Fig. 


_ which was obtained using the second-approximation results for the compres- 7 


ston problem? and using Eq. 15 nal ‘The quite similar but much 4 
ing An extension of this work the first using more 
cated expressions than Eq. 15 for w, have indicated that the results for shear . 


stiffness presented in Fig. 17 are about 5% two- and- -one- ~half times the 
kling load. 


Such a_ procedure as that just outlined using the first approxi ion2 for the 
, ta 1) : 
— 
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Mario G. Salvadori, ASCE, ‘reer Skalak,? and 


A Paul Weldtings 


‘The propagation of waves and shocks in various inelastic media are 
= studied for a semi-infinite body loaded at its free surface. aes “a v 
ot _ The first type of material postulated is a locking medium, which responds 
elasto- ~plastically at low stress levels behaves asa rigid body after it 


& may generate either a supersonic shock wave, or 3 a subsonic wave with an | 
_ elastic precursor, or a purely elastic wave. All three phenomena may occur — 
_ in the foregoing sequence. . Criteria for the initiation, e existence, and decay of 


bs The second type of medium studied exhibits a constant modulus of mel 


cements. 
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ioe 2 Assoc. Prof. of Civ. aa Columbia University, New York - 
3 Cons. Engr., New 


— 
— 
— 
— — 
— — 
— — 
q 
as 
=: ts id i are of stress and strain under a pressure pulse is given and it is shownthat even- _ is 
tually the medium_will reach an unstressed state and remain at rest but vith 
|, 
= 


INTRODUCTION 
if 
: Tht paper is italiana with the behavior a of ma 1 Can 
7% compacted (or locked) under certain types of loadings. A compacted layer 
in such a medium is characterized by a response which can be described as | 
hea rigid, as compared to the rusponse of the uncompacted or loose por- a 7 
tion. In a first type of medium, the compaction is initiated abruptly as a criti- 
E 2 strain is reached, as shown by the stress-strain diagram of Fig. 1. Com- ee. 
oe is initiated gradually in a second type of material, which exhibits the a 
os str characteristics shown on Fig. 5 (to be presented subsequently). 


| p _ The the y for the dispersion « of stress waves in media with gene ral stress : / 


strain curves, concave upwards or downwards, has been explored on a mathe-_ 
matical basis,4 and applied to various problems.2}6 = = 
For the most part, the assumption ofa continuous and smooth stress-strain 
t ae makes it difficult to obtain results in closed form. However, the general _ a 
_ types of behavior discussed herein have counterparts in earlier work dealing - a } 
with such smooth stress- strain curves, M. White and L. Griffis used5 
rve similar to Type I locking media, Fig. i. _ They also considered =, 
case whichis similar to Type Il media, Fig. 5. The idealizations of the pres-_ 


ent paper permit closed solutions and 


a ante and the density. For stresses above greater than 09 the stress-strain — 
“1 1 curve is parallel to the stress axis and the medium, therefore, responds as a 


to stresses larger than Cop On unloading a permanent strain le re- 


e limiting cases of te stress- strain relation are of interest: 
(a) = the material elastic up to the strain €,. 
(b) €o = 0; the material is rigid-plastic up to the stress % 
Oo = 0; the material is plastic up to the strain 
4 “Propagation of Plastic Deformation in Solids,” by Th. Von Karman and P. E. Dawez, 
roc. VI, Int. Congr. for Applied Mechanics, Paris, 1946.0 - 
_ 5 The Propagation of Plasticity in Uniaxial Compression,” by M. P. White and 
Jour. Appl. Mech., Vol. 15, 1948, p. 256. 


“Permanent Strain in a Uniform Bar due to Longitudinal Impact,” . ae M. P. White 
L. Gri ffis, r. Appl. h,, Vol. ‘14, (1947, 
— i Jou ppl Mec ol P. Aj 


¥ 
| 
— material which, under a compressive load, has the stress-strain 
4: 
— 
— 
— 
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dium which exhibits a ‘definite maximum strain has been termed? ; a om 


"stocking material. ” In various practical conditions real materials approxi- 
mate this behavior, as mentioned? by W. Prager. Granular media, such as ” 
Certain types of soils, also behave like this and are of particular interest. Up- a "a 
on initial loading, the response of such materials is nearly elastic8 but further 
increase of the compression produces yielding and breakup of the elastic 
‘ grains, which are then rearranged in a denser packing. This results in an ir- - 
reversible compaction; a sudden increase in stress is manifested, beyond * 
which only negligible increase of strain and density will occur. These phe 
momen are idealized in the behavior of the locking medium postulated. = > 
OF A SEMI- INFINITE. LOCKING MEDIUM 


applied to the y = 0,a pressure wave of finite amplitude fewill 
propagate with a plane wave front. In what follows, stress, density, and particle — 
velocities are assumed to have finite discontinuities across the wave front, and — 
the analysis of the wave propagation is carried out for finite ata tetameeate Ja 


..... 
— response of the medium with the stress-strain characteristics, shown > 


Po in Fig. 1, may be e expected to to be influenced by the magnitude ¢ of the initial value = 
Po of the applied pressure. Three distinct types of response are anticipated: 


vave with acoustic ‘velocity will occur if: 


An elastic wave of constant stress by a subsonic compaction a 


front, behind v which the m material is is located, occur if: 


A supersonic compaction fro front, behing the is with 


out any elastic snes wi ill oce occur alll 


in which Per is a critical value of ‘the initial pressure Po which separates tes (b) z 
from (c). The existence of these three types of responses, and the value of 
ES Pcr is established in subsequent sections. mj will also be shown that undera 
oe decaying pressure pulse, an initially supersonic compaction front degenerates 


into a subsonic front with an elastic precursor, which, in turn, is transformed 


“On Ideal 175. Ww. Transactions 


— 
ag 
i Consider my=0tc 
4 
— 
— 
— 
| 
- 
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ek, 8 “Qn the Theory of Elastic Waves in Granular Substance,” by T. Takahashi and Y. | ae 
lletin, Earthquake Research Inst., Tokyo Univ., Vol. 27, pp. 37-43. 
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s WAVES AND SHOCKS 

into : a aiid elastic wave front. An initially subsonic wave front with an pai 


precursor will be transformed in time into elastic wave 


% ELASTIC WAVES IN A LOCKING MEDIUM 
elastic are by a constant or a decaying pulse, fol- 


De 4 (a) If the applied pressure does not excess 0, Poni edie does not lock, 
“4 but responds elastically. The classical one- dimensional wave equations9 govern — 
_ the propagation, and the pulse is transmitted indefinitely ' without distortion at 
a the acoustic velocity of Eq. 1. _ Thus, for a decaying or a constant pressure 


2 <0 


“the in the elastic range is identical mons that ‘of a Tocking material be- 
cause the medium behind the wave front is ata strain moves as 
rigid | body the constant particle velocity 


O< € <e. 


c) An elastic wave c a compa 


astic wave cannot exist behind a compaction because the 


: sure at the compaction front which is aregeal than ren locks the medium at the 
“strain and compacts it intoa rigidbody. = 
— @) If the applied pressure po > Oo, and the compaction front velocity isless _ 
‘than the acoustic velocity Co, then an elastic precursor with a constant a 


_ This can be seen ia noting that the stress inthe etastic seein cannot be &§ 


to dy are propagated at the acoustic velocity 


Dassematans WAVES IN A LOCKING MEDIUM 


Let the ‘location of the compaction wave front be given by z(t) and the dis- i 

placement of the free boundary by u(t), where z andu are measuredfromthe 

initial position | of the boundary. . The velocity + of | the front is is 2 and that of the 

boundary is u, where dots” indicate derivatives with respect to t. Since the 

medium is locked behind the front, the particle displacements and velocities © 

are alsou and respectively, (Fig. 2) 

_-- The state of the locked medium behind the front is characterized by ett 

- Pe, strain Eos — u(t), and stress o(y, t), where y is the distance from the | ai 


— 
— 
|| 
| | 
| 
— 
— 
— 
— 
generates the particle velocity = j= 
asingstresses 
| 
— 
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iIfz> be the front is supersonic. 7 Then the medium ahead of the front will 
a: be at rest, unstressed, and at initial density p.. If 2 <c, the front is subsonic | 

: ~ and there will be an elastic precursor wave. The material just ahead of the 
compaction front is then at density , Po» ‘Strain. stress 0, and velocity ip = 

— €0¢ Co. To permit the simultaneous treatment of both subsonic a and sil > 
cases it is convenient to define 6= such that: 


i 

‘strain: 


“velocity: 6 8 


Conservation of mass requires, that: 


or in te Nhe 


: "and, since the m medium behind the front is compressed from its original length 


_ The momentum of ‘the mass Pc (zp -u r) between the smeenuiens front and at 
é.. surface is Po (zy - - ur) i,. This momentum changes at a rate equal to the a 


— 
4 &§ 
.& 
im 
f 
— (8b) * 


Substituting ‘Eqs. 1, 8b, and into equation of of motion becomes 
‘ 


(fc £9) 


stress aty) can be obtained from an another form the equation 


(13) 


_is obtained by sub- 
“r 7 


‘aie Eqs. 8b, 9, and 11 into gaa 13: 4 


_ Differentiating the left-hand of Eq. 11, and setting x0) = 0, the we 
locity of the front is given by: 


The erttsenl separating the supersonic from the subsonic 
motion, is obtained by setting 2(0) = =Cy in Eq. : a, which yields for ae ee 


The significance of on the s diagram is shown in Fig. 3. 
SUPERSONIC SHOCK IN A LOCKING MEDIUM 
“When Po > Pers > 1, (6= 0), the of motion 


“O€ 
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Integrating twice with the init initial z2=2z=0 at t=0, the motion 
_ the shockfront is given by: we Or 
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iategrals are taken between | the limits: 0 unless otherwis indicated. 


— 
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rticle velocity i is same 


_ As the surface pressure i denen. , the front velocity 2 decreases, Let ttgb 

defined as the time at which: 


| dt at | 


t 


ak 


t 


particle and ac are given a at t= 


4 


s 
s “cr 


— 
«Using 
— 
— 
and the stress at the shock front 
Also at ts the acceleration — 
(28) 
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and, ‘therefore, the shock front does not degenerate into an elastic wave, The 


motion must continue a subsonic compaction preceded by an elastic 


_SUBSONIC FRONT IN A LOCKING MEDIUM 


‘becomes: 


“gard 


at r/ 


= 


1/2 


initial al velocity of the front from Eq. 15 is: 


a time t = initial for the subsonic are: 


a 
— ubsonic front, ivative can be 
Since both the s Z_ and 


‘This transformation introduces a new time t -ts+ and of 
motion becomes Eq. 30 with t, Zr, 2, replaced t, as Zz 
_and , may be integrated twice. “(The solu- = 


solution by setting tg = = 0. 
motion of the front ‘is given by: 4 


4 


2 fat at 


beg 


2 fat Sp(t) * 


‘The stress at the front, from 14, iss an 


and the particle velocity from Eq. 9 is: 


oT: 
gas 
is s obtained aby Eq. 37 


— 
— 
— 
— 
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— 
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— 
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; At the time t = tg, the stress o z at the front, the particle velocity, and the par- 


ticle acceleration are equal to “the values derived for the end of the s supersonic 


‘motion given by Eqs. 26, 28, and 29. ‘The | front acceleration however, eeaoenae a 


and the a change ir in slope at that time 
=. _ Let time t=te be the instant when the motion of the compaction front, ‘rela- 
tive to the medium ahead, ceases: 
2 


Eq. MRS 


, the wiatinn velocity also has the valu 


The stress at the tront dec ‘to he limiting stress at t = = te: 


‘The ‘equations ‘deriv in in the are also valid for the limit- 

“ing cases described under the heading “Physical Characteristics of Type “e 


im 
— 
his 
The entire locked mass, therefore, moves at the particle velocity of the pre- 
cursor. The time te is defined by Eq. 40, that is, by: 
3 } and the compaction ceases, while the elastic wave continues to propagate in-— im ig 
— 


= Supersonic compaction n exists if Po > Do | and the ‘equation. of motion with k2 


3 


4 ‘If po <a el elastic wave propagation will occur, and ‘therefore phe- 


5) 


Case (b): 


4 


strées Jo and the compaction front propagation is + gove rned by 


= 


the pressure is instantaneously propagated to without com-_ 


conditions imply = 0, therefore the is 


a defined in the dimensional form 


att 


ne 


4 
iim 
— 
— 
asticwaves 
— 

nal time T is given in terms of atime constantt,: 


dimensional forms: 


ib 


= 


where primes derivatives with respect to T. 4 shows the 
sults of a numerical example in of the non- -dimensional for 


PHYSICAL C CHARACTERISTICS OF TYPE I II MEDIA- 
second medium is defined by the stress-strain diagram of Fig. 5. 
. a stress-strain curve is a straight line OA on initial loading to A. The slope 


of this line defines an initial modulus, E Upon ‘unloading the stress-strain 7 
curve is another straight line AB which dete a second modulus, Ej. If the 
: _ material is reloaded, it follows line BA toA and then continues along the i in- 
2 


a _ ‘iia a plane compression wave propagating into a virgin material. The 


‘ distribution of pressure ata given instant is shown in Fig. 6. ; Point A repre- 


_ sents the leading edge of the wave propagating in the positive y direction. In 
the portion of the wave AB the pressure is continuously | and, 


mies, and stresses may be expressed in the non- 

Fig. 4 shows the results up to the time T= 2.0, the time at which compac- 
ceases. After this time the system behaves as a rigid mass (the 
* 
| 
7 

Throughout the remainder © paper only small rains are 

> — 
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portion is “indistinguishable the wave that wc would result in an n an elastic 
terial of modulus: Thus, this portion of the wave will tend to propagate 
distortion!0 atthe velocityco: 


‘The one- dimensional v wave equations» 


particle in the. — of wave ‘Propagation, and u is the 


"TYPE i ‘VIRGIN MATERIAL 


In the wn pee thet behavior isthe same as that of a linear elastic ma- ~ 
pet in which the stress, particle  eraaeais particle displacement, and strain 
3 ‘In the portion of the wave BC, ‘Fig. 6, the stress is less than at B. The 7 
strain values in this region will lie on an curve as AB 


ogee 
if 
where pt ion in the zone 3 
 &§ 


SPACE-TIME DIAGRAM FOR DISCONTINUOUS 


Since a se Bat is greater ‘than Cor the rear portion of the wave, BC, 
tends to overtake the forward portion, AB. There is a continuous reflection — 
“Process at the functure B as will be shown subsequently. 


‘semi- ‘is to » y= Otoy=o. Ifa pres— 
- p(t) is applied to the surface fe 0 at time t = 0, there will be a wave front vs 


at y= o t for allt>0. This implies »s that the material below the line OA in - 


‘line OA is the initial sound d velocity hoe 


— 
— 
— 
a a In the unloading phase the particle velocity and stress obey the wave equa- ig aa 
— ] tions 8,9 since the stress-strain curve is linear. The displacement and strain — ii 
not governed by the elastic wave equation because the unloading stress- 
m 

|. 


__The response for 0O<t<t ys a step wave of press pressure po propagating ‘ at 
Co. The conditions in 0 of Fig. 7: are = Po and particle 


velocity Ug: 


At time ty the om in surface pressure f from o to Pi occurs at the ie 


4 = 0 and propagates at bnateatie: cy which is the slope of the line separating 


The time tg at which the pressure pj reaches the mT wave front may 
be determined fr from the of Fig. 7: 


= At the time tg r= entire length of material from y=0 to y= Co tg is at 
stress = py and particle velocity It may be considered that this materi- 
3 is a body of elasticity E; striking the undisturbed mass of elasticity Ep. A 


wave is propagated inte the undisturbed material at velocity co anda reflected 


_ wave travels upwards? at t velocity | c1. The reflected wave front is enema 
by the foundary between zones 1 and 2, Fig. 
_ The stress gg and particle velocity tg throughont zone 2m 


— ressure history consider The value py — 
— The particular s p(t) rises abruptly from 
0, the surface pressure 
iy: 
iz 
j 
— 
— — 
— 
= 
a 4 a 
— 
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ae: The Salis wave reaches the surface y= 0 at time 2 tg which i is determined — 


q 
Conditions 3 at t time a are bala ‘adi at t time x The stress and ve- 
lip at time are e replaced by fra at Hence the e values of 


y 


may be continued indefinitely and leads to the following general expressions a 


may be analyzed in the same way as the reflection at time tz. The — 


— — 
— 
— 
ta: 
— 
4 
— 
— 
— 
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. The ‘limiting particle bs from both Eqs. 68 and 71 is: 


— ‘the solution for values of t much larger than ty to 
final pressure throughout with no memory of Po as 
as stress and velocity distributions are concerned. 
n Cy Co , (purely elastic material) one obtains -1l,a- -0 and 
The surface pressure reduction from po to pj never overtakes the original — 
= front and the boundary between Zones 0 and 1 is parallel to OA in Fig. 7 = 
When | cj approaches o while Co is constant, 8 approaches 1 from above and 
i approaches 1 from below. — _ The physical action approaches that of a medium — 
which locks _irrecoverably at the maximum strain reached at each point. As 
the limit cy = «© is approached, the zones in Fig. 7 collapse i in the t-direction. 
Solving Eq. 66 for an/2; 


ab 


oO 

Eq. 73 in 67 and 68, one obtains for n even 


— 
— 
ae 
— 
— it area exerted by the pressure onthe surface. The al 
_ is the impulse per unit ar ical with those which may be derived by adirect — a. 
results Eqs. 74 and 75 are identica wit he fact that the compacted mass moves 


WAVES AND : 


front may seen by reference to the zone pattern in Fig. 7, 


= 


‘The 
ies mee: For cahite not on the surface nor at the wave front, the value of on is an odd 
s or even zone value depending on the time. _ The fraction of the time when even | 
values are in effectis equal tothe distance .from the surface to the point divid- 
ed by the distance from the surface to the wave front. During the remainder — 
qi the time, odd zone values are extant. ‘ Thus, in the limit the stress fluctu- oon 
q = with infinite frequency at each intermediate point. By averaging the fluc 
at is seen that the actual value of the stress is: 
)=o,+ y [e +e 


| This result may y also be derived by secrete’ of the equation of ‘motion — 


fe) 
wave analysis this section can be written in the foltowing form: 

wr 


d 
F,( nists Plt) 


ats 


a 
. 
F ‘(31 


se stress and odditinn in each new zone of Fig. 7 to be made up of values a hel 


: a the previous zone plus a contribution f, due to the wave traveling upwards ox 
at velocity, The term Fo t + represents the sum of all 


4 


-— contributes the original wave 


ee pressure Pp and the reduction starting at a: Fig. 1. Bea 


_ RESPONSE TO CONTINUOUSLY VARIABLE SURFACE PRESSURE 
a decreasing surface pressure such as in Fig. pres- 
; sure p(t) rises from zero to pp att=0. Fort> 0, the pressure decays con- 7 
to zero, This surface pressure be regarded | asa Step wave of 


~ 
Po series of infinitesimal, negative dp = 
bal pressure change, dp, may therefore be treated in the manner of the singe ores- 


fects will yield the complete solution for the variable pressure Pp. 
i = _ When the surface pressure shown in Fig. 8 is applied to the boundary of a 
- semi-infinite type II locking medium, there will result fil erat aterm like Fy in in 


‘this series will be again a series like Eq. 83. 
of form Eq. 84, will 


86 represents the integration all of 


— 
a 
4 
— 
— 


Pglt) = Po - .. 


Bas 
os ‘The value of Pa 7 shown ¢ graphically in Fig. 8. 
In a similar manner, contributions of the form fn(t), Eq. 84, may be inte- P 


_ grated in the present case to yield ee 


(t) = a" p, 


to F(t), Eq. 83: 


The typical form of pp(tP is shown in ‘Fig. 
Lh The complete solution for the stress in the present case has a form simi- = 


(91) 


_ The ® velocity, in in the present case, is given by an equation of the form of small 


6 Eqs. - and 92 give the stress and velocity for any y and t in the range 0 a 
y<Cog The strains at any y andt and t may be derived from the stress at 


int in uestion. 
the point in question. 


TION, Pr 


tere pg(t) is a) 
— 
— 
a 
= 


face pressure v(t) decays to zero ina finite time. In such cases, the pressure ~<a 
_ reduction pg rises se Po in a finite time and the reflected pressure is — 


-a@ 


ero for long times. This: is to be expected a dissipative material. 


‘The limiting case a exhibits no 
on unloading. The formulas for stress and particle velocity, Eqs. 91 and eal 
reduce to simple forms in this case. 

By use of Eqs. 88, the particle Eq. 92, is in general 


where has been replaced by Lo 62 


cy 2, then a1 from below and (1 - 


inn, change Therefore, the sum may be by an 
in the limit. The summation es p in Eq. 96 may then be einnienean as i 
tine 
(a t) maf -a)a p(a@ t)dn. 


— 
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| 
j 
- 
iii 
— 
pow. | — 
+22 
| 
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— 
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Changing th the variable of Antegration n to where ald 


ta 
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L? 
"where is is the per | unit area exerted by the surface pressure during 
: The sum containing Po in Eq. 97 is evaluated by re- 
by in on — 99 in Eq. 97 yields the final result: 


fashion to ‘the particle ‘in the case cj +. tThe 


pars the sums ‘sums by nee and integrating in a ae manner as 


* 


above, yields 
w 


ry 
lim te = lim 


4 


Ae 


— -=G 


— 
Eq. 100 shows that the entire mass from the suriace to the wave front at 
he _--y = ¢, t moves as a rigid body, since the particle velocity is independent of y. $i 
— 
— 
hh 
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car Now as a— ~ the values of si and t' both approach t, so that ~~ 105 becomes 


result of Eq. 1 106 6 may be written in terms of the acceleration if 


‘The stress at the surface y = 


(0, t) = = 
and at at the wave front: _ 


The given by Eqs. “100, 107 and 109° "may y also directly 


by writing the — of motion of the compacted mass which moves as a 


Let the 


. (110) 
where Po is the initial peak intensity of the and 6 isa ti isa time constant, 


— 
The properties of the medium are as fol 


Velocity | at Unloading, 


Recovery) Eqs. 110 and 88 


‘where ais given by Eq. 62. 
Using the formulas of Eqs. 91 92 the stress 

, “Gcan be computed at any point and time after the two functions p(t), Eq. 110, _ 
ve Pr, Eq. 112, have been tabulated. © Fig. 10(a) shows p(t) which is also the 
_ gtress at surface and o (c, t, t), the stress at the front. Fig. 10(b) gives etl 

particle velocities at the surface u(0,t) and at the front, u(cg t, t) “al Ag oN 
b: Propagation Velocity at Unloading, oe (No Elastic 

‘By Eqs. 110, 108 and 109 the stress given 1 by Loe 


and 4 t) and ‘consequently o 1 t, indicating that the material with» a 

no recovery gives an average value of particle velocities and a lower bound ss 


for the peak stress. In both instances particle velocities and stresses | 


= = acoustic velocities; ‘ 


= ‘Th 


stress wave » front; 


= ordinates of the locked (Fig. 2); 


4 
these values are also plotted in 


—Gj (STRESS AT WAVE FRONT, MATERIAL WITH NO RECOVERY) | 
— 


— 


8 


& 


u(0,t) (VELOCITY AT SURFACE) ~~ 


=~ | 


ecaying pressure pulse; 
initial value of p(t); 
= critical initial 
= non-dimensional pressure; 
non-dimensional stress; 
= non-dimensional 


= non- Gimensional time; 


= elastic modulus during unloading and. reloading r Fig. 5); 

a 


= 
impulse of 


WAVES d of the wave front; — 
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COMMENTARY ON PLASTIC DESIGN IN STEEL: CONNECTIONS 
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the sixth in a series of reports ¢ on Plastic Design ‘emanating 


‘iia paper is 
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(which appeared i in the July, 1959 issue of the Journal of the Engineering — f 
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‘pears at the beginning of this Progress Report . Chapter 8, entitled “Connec- 
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— 
— Proceedings of the American Society of Civil Engineers _ 
— 
— 
= i 
* 
— 
4 


comMenrARY ON PLASTIC DESIGN IN STEEL 


2. BASIC PRINCIPLES 1 ‘Behavior of Material and 


3. 3 Mechanism Method of Analy sis 
Structural Ductility 
‘5 Yield Stress Level 
46° Plastic Moment 
1 loads 
a 8 Load Factors 
VERIFICATION OF PLASTIC 1 Basic Concepts 


_ ‘THEORY .2 Continuous Beams 
ADDITIONAL 6.1 Shear Force ia, 
. COMPRESSION MEMBERS 7.1 Introduction 
Reduction of the Plastic 
to Axial Thrust 
1A Rotation Capacit 
15 The Influence of Lateral q 
8a 1 Straight Corner ( Connection 
Haunches 
8.3 Tapered Haunches _ 
4 Curved Haunches 


8. 
8.5 Beam-to-Column Connections 
8.6 Details with Regard to ) Welding 
Details with Regard to Bolting 
9.1 

9. 


e- Deflections in the Elastic Range 


9. 4 Deflection: in the Plastic Range 


9.4 Deflection at Ultimate Load 

9. 5 Step-by-Step Calculations = © 

9.6 Approximate Deflection at = 


: 
ice 
— 
i 


7 
Connections play a key role in assuring that a structure can reach the co 
-_- ultimate load. Since connections frequently are ineshed at points of a - 


in design. The principal requirements for connections are: 


. Overall stiffness for or maintaining g the of all structural units: rela- 
_ tive to each other. 


on 


The various types of connections which will be discussed a: and whichare typi- 


a calof those that might be be encountered i in steel framed cae are designated 


Column pechoroge 
Miscellaneous 


8.1.— 


in Fig. 8.1. “These include corner connections haunched), 
to-column connections, beam-to-beam connections, splices, c column anchor- 
ages, and miscellaneous connections (purlins, girts, an and bracing). 
_ Primary attention is given herein to corner connections and beam- -to-— 

~ column connections. Methods of analysis are based on assumptions of stress 
: distribution at ultimate load which satisfy equilibrium but do not violate the — 

condition. Solutions: thus” constitute lower bounds to connection 


capacity. The same principles would be betta however, for ‘the analysis: 


ott other connection types. 


— 
— 
— 
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_ Straight corner connections are formed by Seadtie joining two rolled sec- 
tions. _ Straight corner connections in which the rolled sections are joined at i> 


right angles (as at 1 in Fig. 8.1) are sometimes called square corner connec- — 
tions. Studies of the theory, design, and behavior of square corner connections 
may be found in Refs. 8.1 to 8.6. The basic principles of the theory of con- 
ections will be illustrated by considering an unreinforced square ‘corner con- 
ection. It is more critical than a connectionin which the members do not join 
t right angles. Fig. 8.2(a) shows the moment diagram for a typical rectangu- 


lar frame loaded with a a uniformly distributed load. typical unreinforced- 


n 
n 
a 


Moment Diagram for. 
Portal Frame 


Forces in in Corner 


corner connection is sketched in in Fig. 8. 8.200). shear 
acting on the connection are depicted in Fig. - 8.2 -2(c). 
md In arriving at a simple analysis of the forces in a knee it is assumed that 
normal stresses caused by bending moment and thrust are all carried in the 
o flanges and that shear stresses are all carried by the web. In Fig. 8. 2(d), the 
action of the applied forces on the parts of the corner and of the parts on each > 


_ other are represented by arrows. The tensile force in the outer flange of the 


— 
— 
aa 
— 
— 


“4 beam. The resultant of these two forces is carried into the corner as a prone 
q along line DC. A similar pair of force components exerted on the ied 
stiffener | causes shear along line BC of the w web. Live 


corner connection, it is rather obvious that an unsatisfactory condition would | | 
exist: if there were insufficient material to carry the forces without Suckling 


general yielding. Assuming that the horizontal rolled ‘section | 
through the knee, its flanges AB and CD of Fig. 8.2(d) are sufficient to resist 


_ the flange of the column. - The vertical stiffener BC must have sufficient area _ 
; to carry the column compression flange force into the beam web. In all cases _ 
the welds must be sufficient to transmit the required shear or normal force. | 
Study of the square web panel reveals that the shear forces would tend to de- 
form the panel as shown in 
Consideration of the equilibrium of the horizontal forces on the portion of 
“the outer flange between A and B in — 8. 2(d) will give an expression for the ; 


Ty w 


depth o of column 


Equating the shear and flange the r requ 


on for wy is to in} Ref. 8. 1. It hes 
been shown that the use of ‘more exact analysis does alter the results of 


— PLASTIC DESIGN 1 
beam is carried into the web in shear along line AB. Ina like manner,the ten- 
sileforce inthe outer flange of the column is carried through the endplate into 
the web as a shear along line AD. In each case, the tensile stress in the flange 
— 1 .-_—s is: assumed to be linearly reduced from oy at the edge of the corner (BorD) | ie. 
to zero at the external corner i 
—— 7 The prolongation of the inner flange of the beam carriestwoexternal forces: _ oa. 
— 
— 
intheflangeisgivenby 
Bmust notexceed 
= ired web w: oo 
os — _ According to the Mises yield criterion the limiting shear yield stress Ty — os i 
OyA/3. Using this criterion, the required web thickness becomes 
— 


‘ which case inieoatenael is required. This may take the form ofa doubler 
plate which increases the total thickness of the web to the required amount. 

_ However, it is nearly always more practical to provide a symmetrical pair of 
diagonal stiffeners. — Diagonal stiffeners act somewhat like the diagonals of a 
truss panel in preventing shear deformation. Fig. 8.2(f) shows diagonal stiff- 
eners between corners A and c Of the web. The stiffeners are 


= dy/Ae) 


Noting that Tr must equal Ty + + Ts the area is found as 


The actual sete of straight ‘corner connections has been studied in os 
E several tests. In Fig. 8.3 are shown the results of tests on four straight con- — 
— A, K, L, and M madefrom 8B13 members (8.1), The moment-rotation A 
curves of the connections are compared with the moment-rotation curve for an 4 
4 8B13 ‘beam as shown by the heavy line. «Of major interest are strength, stiff- 
‘ness, and rotation capacity. Each connection reached a maximum moment 
greater than that of an 8B13 beam. The rotation of each of the knees was great 
enough to be considered adequate to allow a structure to form a mechanism. 
The initial stiffness of the knees was approximately the same as that of the 
7 “rolled beam, but larger rotations occurred in the knees at a lower moment be- a 
- cause of residual stresses. However the larger rotations were not severe > 


enough to impair the practical effectiveness of the cor ctions. — 


Results of five additional tests on straight corner connections are given in 
Fig. 8.4. These tests were performed on knees joining several sizes of wide- 
flange (8.4), ‘The shaded zone at Mp/My = Mp is so marked to depict 
the range of. spread in ‘shape factor of the five wide flange shapes that were 
tested. These tests confirm that these knees also satisfied the requirements 


q ody 
ig q 
y considering equilibrium of the forces on the top flange, the requiredarea 
of the diagonal sti Ss mav be obtained. Th 
— 
— 
=A 
= 

4 
4 
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Web Stiffening; 


_— the web of a connection is deficient in aimee. it may be rein- — 

forced by doubler plates to meet the requirements of Eq. 8.1 or by di- 
agonal stiffeners welded to the flanges and to the web. The area of a 
symmetrical pair of Gagonal stiffeners should 


/a 


8. 2 HAUNCHES 


Haunches of either the tapered or r the curved are ‘sometimes used to 


achieve a a pleasing appearance. Their use in elastic design makes it possible 
“to adapt the section of the haunch to follow more closely the shape of the mo- 
ment diagram, furnishing approximately the resisting moment required at a 
of given sections with resulting economies of material. . Similarly, in 
- plastic design the use of haunches also makes possible a reduction in size of 
the main member, which may mean the difference — using a built-up 
member and using a more economical rolled shape. - Although the haunch will 
"Permit a reduction in main member size, it will be costly to fabricate and thus 

_ may otfeet some ot the savings that will be realized by using the smaller main — 


Tests of naenaiod connections designed by methods intended primarily for 
_ elastic structures revealed that the connections exhibited good strength, but — 
that some lacked sufficient rotation capacity. The lack of rotation capacity 
: attributed to premature lateral buckling (8-1). Recently developed methods _ 
: of analysis confirmed by tests show that plastic hinges can function properly 
in the haunch if adequate provision is made to prevent such buckling (8. 7 a ia 
“= The effect of haunches on the analysis of a frame is to increase the number 
of sections at which plastic hinges may form. However, the methods of analy-_ 
sis are unchanged. Fig. 8.5 gives the results of a mechanism analysis of < 
is frame with haunches. The correct solution depends on the loading and 
geometry of the structure. For the given geometry and loading the correct so- 


‘lution would be Mechanism 4 (Fig. 8.5(e)) and the resulting moment diagram 


hinge | moment of the main members i smaller as a result of 

moment Mp, whichis greater than this From the final moment dia- 

gram (Fig. 8.5(g)) the moment, thrust, and shear at any cross section of he 


may be determined for purposes of design. 


— 
— STRAIGHT CORNER CO 
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— 
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: OA sketch of a typical tapered haunch is is given in 1 Fig. 8. 6. The aii to be 
considered in the design of the haunch are the moment, thrust, and shear “4 
i Section 1 where the rolled section is joined to the haunch. ‘These would be de- d oe, 


termined from an analysis of the complete structure. “ss 
Conn} 
Conn. 
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- kips 


in - ki 
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‘All Connections 
8 B13 Members 


“Average ‘unit re rotation in radians per | inch 
FIG. 8.3.- —MOMENT-ROTATION CURVES OF 


_ Three mainproblems must be considered in the design of tapered haunches. 


‘These are: — (a) resistance to bending in the tapered portion of the haunch, (2) 

“4 resistance to local and lateral buckling, and (3) shear stresses a the web and 
flange forces in and around the corner panel BDFE (Fig. 8.6), 

ee In analyzing the tapered portion of the haunch for resistance to bending, 

4 :: ‘ser sections perpendicular to the layout line of the structure may be con-. 


sidered. The web thickness and width of the are made 


68.3 TAPERED 
q 
— — — 
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oni to those of the adjoining rolled section. This insures that the webis able 
to carry at least as much shear and axial force as the webof the rolledsection. _ 
‘Changes in the resistance to bending moment of the haunch may be controlled, a 
therefore. , by the thickness of the flanges and by the depth of the haunch. | The — 
outer flange is usually parallel tothe layout line and the inner flange makes an 
angle o of taper 8 with the layout line, thus defining the depth of the haunchat any 
¥ os At any distance x along the haunch, the moment of resistance of the 
section perpendicular to the layout line must equal or exceed the — 
om paar at that section determined from the plastic analysis. = 
a ‘The he plastic m moment of ai catalial of any cross section i is oo wen by (8.1 
where Zx is the plastic 1 modulus at the particular section. It has been shown 
that if the plastic moment of resistance is adequate at the rolled section (sec- i | 


iG. ve 4. —MOMENT- ROTATION CURVES FOR SIZE- 


tion 1 Fig. 8. ae and at the re-entrant corner (section 2, Fig. 8.6), there will | 
re no necessity of checking at other sections along the haunch (8. 7), ‘oa 
nfs In proportioning tapered haunches, it is desirable that the flanges have the a 


game size as the nominal dimensions of the adjoining rolled section. In this 


case, the flange area presented on a cross section perpendicular to the layout 
; line would be different for a sloping flange at angle 8 than for a flange parallel 


to the layout line. _ The effective cross section would be unsymmetrical with i 
respect to its major axis with resulting complexities in the calculation of the 4 


— 
— 
> 
— rr «= | 
— 
— 


| 
i In Ref. 8. 1, a somewhat simplified analysis is suggested to reduce the com- 
.. plexities: of computations. As shown in Fig. 8.6, ‘the stress in the sloping com-— 
 yeeeaien flange is assumed to act in the direction of the flange. To calculate + 


(c) Mechanism 


Mechanism (e) “Mechanism 4 


(g) Diagram for 

FIG. 8.5.-ILLUSTRATIVE EXAMPLE—ANALYSIS OF PORTAL 


increasing the of the sloping flange, the 
in each flange m may be made equal, thus providing effective symmetry to the 
eross section. The necessary thickness the sloping flange is 


‘bg = thickness of (compression) flange plate 


— — 
(d) 
— 
— 
| 
a 


The modulus at any section will be 


ob = width of | 

dy = depth of haunch at section x 


web — 


cases where the is is less than 20°, the required increase in thick- 
“y ness suggested in Eq. 8.4 would not exceed 6%. Most practical design ‘pro-— 
cedures: would neglect a difference. For all practical purposes, “bot 


x >. 


‘ 


8.6.—FORCES ON TAPERED KNEE = 
flanges should probably be ‘made equal and still use the assum; assumption that 
_ cross section is symmetrical and that its centroid is at mid- a 


assured by adhering to the ‘recommendations of Chapter 6. pers 
Lateral buckling of tapered haunches depends onthe slenderness. of the com 
- pression flange and on the extent to which the flange is yielded and strain 
hardened. It is assumed that positive lateral support of the inner flange will 
be provided at points A, B,and C. (Fig. 8.6). Since the web provides restraint — 
against buckling of the flange about its weak axis, the flange must buckle in its i 
conservative approximate expression. for the critical buckling length 


has been obtained 8.7 by using the tangent modulus 


— 
Ki. 
> 
= 
— — 
- 
1 
7 


buckling concept (8.8, 8. 14) 


It is assumed that the compression flange is uni- 
_formly stressed to 0, , that the strain in the flange just reaches € st» and that os 
_ the flange buckles as a pin-ended column between points of lateral bracing. | 


These result in acritical slenderness of: — 


= 


Ty == radius of gyration of flange = 

. Est = strain hardening modulus of steel (900 ksi for A7 steel) 
yield stress s of steel assumed in the design (33 ket) 
By substituting the material properties, Est and Oy, the following critical Tength 


is obtained for fully strain- hardened compression flanges of A7 steel: 
Ler = (8. 8) 
‘This expression a any restraint to buckling offered by adjacent portions 
of the flange or by the web. 
Frequently it will be found that oo given by Eq. 8.8 will be less than the 
_ value obtained from the original haunch layout. In such a case, additional points Ls 
a of lateral support could be provided. Also the flange thickness or the depth of — - 
the haunch at Section 2 (Fig. 8.6) could be increased. Reference 8.7 suggests | 
| for accomplishing such modifications. = 
The shear stresses in the web and flange forces around corner ner panel 
BDFE (Fig. 8.6) are of the same character as those shown in Fig. 8.2(d) for ._ . q 
square corner, and would result in deformation similar to that in Fig. 8.2(e) 
2 the forces exceed the shear strength of the panel. The problem may be 7 


- examined in more general form with the aid of the sketch of a haunched knee ¥ 


of a gable bent shown in Fig. 8.7. In this knee, the forces on the panel BDFE a 
are to be determined. For the most severe loading, the cross sections BD and | 

' _BE will be fully yielded and the total area of the flanges at these sections will df 

be stressed to Oy. . At the re-entrant corner B, the inner flanges will have an 
_ unbalanced component of force in the direction of BF in the absence of a web | 


stiffener. . This unbalanced component of flange force would have to be carried» 


z od by the web. The e web could be fully yielded due to bending stresses at point B, 


ae soit is advisable to neglect the carrying capacity | of the web and furnish a sym- 
metrical pair of diagonal stiffeners along FB. = 
i At the outer corner F, the situation is slightly different because part of the a 
q , flange force could be carried by shear in the web along lines EF and DF thus 
reducing the components of flange force to be carried by the diagonal stiffeners 


atl F. For this reason, the selection n of stiffeners will be 


— 
7," 
i 
: 
| 
— 
— 
At 
4a 
Nog 
Ty 4 
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‘The maximum possible inthe and in the diagonal stiffener 
at corner B are shown in Fig. 8.7(a). The required area of the diagonal stiffen- 
er may be determined from an equation of equilibrium of horizontal il components, = 
ofthese forces asfollowss 


As cos Act 08 Oy Ace sin Bo = Bo =0.. 9) 


> 


cos(f,+ 7) - Aca sin Sin 


re 


— 


‘of inner flange o o sr haunch 
le of ti aft ate 
e of inclina ion of ‘raft er 
= angle of slope of diagonal stiffener 2 


angle of taper of 1 rafter 


Ati cos y os cos °° cos + cos 


a 


os2 


= area of outer flange of beam haunch 
= web thickness of 
deh ch = depth of column haunch at section EB (Fig. 8.7 


is that the ‘shear capacity of the web is adequate to transmit the outer flange 
oS force, and that diagonal stiffeners are needed only to transmit the unbalance of 
S the inner flange force. Since in any case the web will carry some force, Eq. 
8.10 based on equilibrium at the inner corner will control for es ileal 
= Transverse stiffeners at the junctions ot of the tapered and sainmaenadibies 
_ (points A and C) may be designed to > carry the unbalance of the flange force due 
to the sudden change in direction. Since the angles 81 and $2 are rcnchle 
— quite small, the size of the sitinete usually will be governed by” minimum — 
A series of tests of haunched corner connections has given experimental 
‘results agreeing well we the concepts" (8. 9), Moment-— 


q 
— 
- | = 
— 
_ where | ‘gal 
— 
— 
7 
= 
4 


a _ rotation curves s for these tests are shown in Fig. 8.8. Specimens 44 and 45 were ne 


if made with the critical lengthof compression flange as given by Eq. 8.8. Speci- 
= =€68©men 44 was not braced against twisting at the inner corner while specimen 45 
: was braced. Specimen 45 performed satisfactorily while specimen 44 buckled | 
_ before it reached the computed ultimate load although it | it exhibited considerable 
post strength. Specimens = and 47 had compression flange lengths: 


8. 1. DIAGONAL. |STIFFENER F FORCES IN IN A GAI GABLED HAUNCHED KNEE Gell: 


greater sreater than that that Eq. 8. 8, ‘but were m modified by increasing the haunch 
_ depth and by increasing the flange thicknesses respectively. These specimens — 
also performed satisfactorily in that they met the requirements stated earlier. . 
‘The results of these tests indicate that safe and adequate connections om 


_be made by meeting the requirements derived by the theory. In summary, the 
zi following suggestions are made with vages to proportioning of — — 


120° 
— 
— 
| 
— 
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Re quired Plastic Modulus: 


‘ ‘moment at the Ae end and the shallow end of each haunch. _ 


of Chapter | 6 should be followed to assure that 
ture local | buckling of flanges will not occur. 


The « critical length of ac flange which be 
strained to the point of strain-hardening throughout its entire teats es 
lateral buckling is conservatively estimated as: 


_ The critical unbraced length of a compression mn flange may be in- 
creased substantially by increasing the — of taper 8 or by increas- 


To resist an unbalance of inner flange ‘reinforce ‘the 


os rae be > provided having a total area: aia 


ve Diagonal stiffeners should be welded to the web and to both iia. : 


_- Transverse stiffeners at the junctions of the tapered haunch and 
prismatic section serve to carry the unbalance of the inner flange force — 

” due to the sudden change in direction. Minimum size requirements will 


usually governthe thickness. 


In ~~ 8.9 are shown the layout t and the applied forces for a typical ee 


ome in the haunched portion as a result of the moment gradient. ~~ ‘q 
a” with tapered haunches, the three main problems to be considered in the 
design of curved haunches are: (1) resistance to bending in the haunch, (2) re- - 

ae local and lateral buckling, and (3) shear stresses in the corner | 


| 
— 
— 
— 
— 
— 
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| 
— 
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compensate for ‘the reduction in effective section due to. the inclination 3 of atthe 
flange and would make the section effectively symmetrical about its half depth. — 
The theoretically required thickness is given by (8. 


te = = ttr/cos B 


angle between point of tangency and id given section of knee 
For all practical cases , the controlling angle B will be small (approximately 
12°) and symmetry can be assumed without increasing the compression flange — 


symmetry, the section at any section xis 


via d, = = depth of haunch at any section (x) 


R-=r radius 0 of cur of in inner flange 


_ for the thickness of the tension flange gives ie - 


x \b-w/ (b-w) 


At some station the by a controlling of the ang B, 
the values of d, and Mx will require a maximum flange thickness. For values 
_ of the radius of flange curvature and distance to the inflection point that oe 


be encountered ordinarily in practice, the controlling angle 8 has been found to | 
The resistance to local buckling of the flanges of curved haunches may be — 
assured by adhering to the recommendations of Chapter 6. eee << 
Lateral stability of curved haunches is dependent upon to 
_ buckling of the compression flange perpendicular to the plane of the web. Refer-_ 7 
= 8.7 contains an approximate solution obtained by analyzing the problem as * 
a ¢ buckling of a curved beam simply supported at the points of tangency and 
7 at the diagonal stiffener. _ To simplify the analysis the curved flange is con- 
- servatively assumed to have a uniform stress distribution over its whole area 
re. and length. The tangent modulus equation for the ‘buckling of fully strain- 


hardened plates gives the critical arc 
flange (8.7,5-8), 


uy ; _ The plastic moment of resistance of any cross section x of the ea 

 —pendicular to the layout line is (Fig. 8.9), 

, = 
q = 

ii 
. 

Ais equation Wwe suDstitutea into the expression 8.3 for — 

— | — any 
il 
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= central angle of curved flange oteonteth points of lateral support. It 
is equal to the angle (45° . - (y/2)) in Fig. 8.10) unless there ; are ad— 
ditional points of lateral s support between A and B. | 


Tx *=rad radius of gyration of compression flange in strong direction = b/W12. _ 


‘The e critical arc when expressed in of the width becomes 


be provided by the 1 segment and since the stress distribution is 
‘not actually uniform along the curved inner flange, a moderate increase could : 
be made in the critical arc length. . For a 90° connection, Eq. 8.8 reduces to: 

R = 6 b 


bing: 
not braced oat A 


| 


FIG, 8.8.-MOMENT CURVES URVES FOR HAUNCHEL ) CONN! NS 
% ED CONNECTIONS | 
are length according to Eq. 8.8. Other possible measures are to o change the 
radius of gyration by increasing the flange width b or by changing theshape of __ 
the flange cross section. - Also, ‘the strains in 1 the flange could be controlled by 


—— pas 49 
— 
— 
q 
ii tor some reason it 1S desired to increase the radius above the vaiue given 
_ (S by Eq. 8.15, additional points of support mav be added to decrease the critical eee 
im” 


Ae 


- An additional problem not pee with tapered haunches is the cross 
bending of the curved flange. * Cross bending is a phenomenon in which radial 
: ‘components ¢ of force inthe curved flanges tend to bend the flange edges toward 
‘the web. It has been shown (8.7) that cross bending will have only a negligible 
_ effect on haunch connection is so 


_ effect of shear force on the web will be discussed with the aid of Fig. " 
8.10 showing a gabled haunch. 7’ the two ‘segments of the haunch bounded by 
- eross sections perpendicular to the layout lines at A and D and at E and C, the 
problem may be dismissed with the observation that the webarea is equal to or 
greater than the web area of the adjacent prismatic beam which | carries the 


SURVED 


_ same shear force. In the corner panel BDFE, the forces are similar to those ~ 
in the tapered haunch described in Article 8. .y and therefore Eq. 8.12 would — 
indicate the necessary stiffener size at the outer corner F. ae ee 
ss sirragrwin reason for using diagonal stiffeners in the corner panel is to re- 


flange. . An ‘approximation of the required stiffener s size ‘may be determined + 
neglecting the radial compression resistance of the web and assuming that the 
diagonal stiffeners resist two components of inner flange force directed 
chords between AB and BC as shown in Fig. 8. som). eae the forces in a 

= 20, A, A, sin (22.5° 
inner 


8.18) 


— 
— 
— 
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PLASTIC DESIGN 
Eq. 8. 18 will tend to control the design of diagonal stiffeners —_— onan curved ™ 
-haunches rather than Eq. 8.12. In the case of steeper gables (y -large), it is 
that Eq. 8.12 may require larger stiffeners 
* Results of tests of curved connections (8.9) are given in Fig. 8.8. - Test 49° 
_ was made on a specimen which had the maximum unsupported flange length 


; _ recommended by Eq. 8.8 and it performed ve: very well. Test 48 was — i 


A 


m 


‘FIG, 8.10. —DIAGONAL STIFFENER FORCES IN A GABLED 
asp specimen with a greater | unsupported co compression ‘flange tos but with 
the flange thickness increased to reduce strains. This specimen also oe 
. _ formes satisfactorily, but had somewhat less rotation capacity. . Results of some | 


— 
— 
im 
— 
i — 
= 
|. 
oo Sa _ us of curvature on rotation capacity. Connection H with R - 5.5 b had somewhat 7 = 
capacity. Connection G w yn capacity, and 


April, 1960 


mi ight be considered i inadequate. Connection | Be of a slightly di different dailies, 


dR R = 3.4 b, and exhibited very good rotation capacity. ia ie 


ase The} plastic m modulus: furnished at any any point in 1 the ry must be ac ade- | 
quate to resist the applied moment at that point. For typical values & : _ 
| R andmoment gradient, the centralangle 8 to the section controlling an 
_ thickness of theflange is about 12°. If the web of the haunch is no thinner 
¢ than the web of the beam, the thickness of flange required to annie a wane 
adequate plastic modulus is given by 


Ae 


+ Lateral Bx | Buckling 


Thec critical unbraced length of a curved compression flange strained 
** to the point of strain- ~hardening without premature lateral buckling is 


Bo The critical unbraced dlength ofa a compressiontangem may b be increased ~ 


To reinforce the web against to in 
the curved inner flange, a of — stiffeners may 


For steep gable 


| be checked by means of Eq. 8.12. a. 


8.5 BEAM-TO-COLUMN CONNECTIONS © 


Fig. 8. 11 shows s three common types of beam-to-column connections used 
i: on planar structures. The function of the “Top” and “Interior” connections is _ 


to transmit moment from one beam to another, the column carrying any un- = 


balanced moment. The “Side” connection transmits beam moment to the upper 


— — 

— 
7 
— 
— 
— 
| 
effect of cross bending on connection behavior will be negli 
= 
— | | 
im 
j 
— 


ati 8.11.—TYPES OF BEAM-TO-COLUMN CONNECTIONS | 


‘ 8. 12. .—FORCES ON BEAM- TO- ~COLUMN CONNECTIONS 


Stiffener ‘Stiffener 


8.13. _STIFFENING OF BEAM TO- COLUMN ‘CONNECTIONS 


— 
— | Peme ty q 
— 
iii 
‘FIG. 
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to carry rryfull- moment beam s without Other 


 _— require stiffening of their webs or flanges to aid in carrying the con- 
= forces from the flanges of the connected beams. = 


ie Studies on the design and the behavior of full sane interior — 


treating problems relating to case of 
a ctions without stiffening will be examined as to the adequacy or inade-_ 


acy of the connections, and then some methods of sides —" necessary stif = 
} 


“a Fig. 8. 12(a) are shown eubemationty the moments and forces on a typi- 


1 


cal interior beam-to-column connection. In Fig. 8.12(b), the effect of a beam — 


the beam web forces being of secondary importance. Significant effects of the © 
beam flange forces can occur in two regions in the column. The first region | 
is the column web where ‘yielding’ may be accompanied by buckling due to the = 
beam compression flange force or by fracture due to the beam tension flange 7 
force. The second region is the column flange where bending may contribute 

q to the fracture of welds cecil the beam flange | to the column flange. . 


_ moment on a column is shown as ; a couple composed of the two flange forces, 


\Zone of Possible Initiation 


«FIG. 8.14. .—COLUMN } FLANGE 


tensile force in the beam flange tends ‘to outstanding 


flanges" as shown in Fig. 8. 14. At the toes of the column flanges, flexibility — 
the beam and column flanges to deform However, at the 


7 
— 
— 
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: in regions of web yielding and flange bending. Stresses in the column ae | 
es 4 by a concentrated beam flange force will spread out as they penetrate into the — 
~ column. ‘Because of this, the intensity of stress decreases with deeper pene- 4 @ 
ts - tration. If the spread of stresses is insufficient to reduce their intensity to - 
= yield level at the depth of the base of the columnflange fillet (k-depth),the web _ 
will not be able to provide sufficient reactive resistance to the beam am flange — 
fore orce. e. This effect is most serious in the region stressed by the beam c com= 
pression flange. Rational analysis of the spreadof stress in a wide-flange sec- a 
_ - _ tion is difficult and is usually replaced by a linear assumption based on test 
Reference 8. 11 gives | the results of ‘several tests and 
that the stress in the column be assumed tobe distributedon a 2.5:1 slope from 
7 the point of contact to the column “k-line.” ’ As shown in Fig. 8.12(b), this as- 
_ sumption implies that the force of a beam flange is resisted by a length of 
_ column web equal to (ty - + 5 Ke) at the column “k- line,” ¢ where | th is the beam > 
_ flange thickness and Kc is the column fillet depth. For equilibrium, the re- 
sistance of the effective area of the web must equal or exceed the applied con- 
-centrated force of the beam tension or compression flange. 


Column Web Resistance = °y We( (te +! 5 Ke). 
Beam Flange Force = “Oy y At 


= area of one e flange of beam 4 


- if the column web resistance is insufficient | tocarry the beam flange forces, : 


_ Static equilibrium in cases where the column flange thickness is small in- 

_ volves the consideration of bending in the column flange due to the beam tension © 7 
or compression flange forces. The analysis is further complicated by the fact — 
_ that the flange bending occurs in two directions, both longitudinal and trans- 
j verse to the axis of the member as can be visualized from Fig. 8.14. To solve 
the problem it is necessary t to make assumptions regarding the distribution of 
the beam flange force on the column —— the extent of the zone of bending 


Beam flange force = oy Af... . bes (8. 
‘a thick center portion of the column flange between the ends of the fillets a 
may be assumed to resist the beam flange force with an it 


‘Column flange direct resistance = = Oy 


— 
- = Oy We (to +5 Ke) 2 Oy (8.21) 
— 
 flanga inine tha ealumn woah Tt maw he acenmad that the yielded beam 
— 
— 
i 


m= = Distance between fillet extremities of flange of ioe 
| The remaining portion of the column flange resistance is due to bending of the 


7 projecting portion of the flange, and is affected by the projecting width of flange, _ 
= thickness of flange, length of line load, length of the zone affected by ie 
_ _ and strength of the material. In general, the action can be looked upon as a 

_ case of plate bending, with the flange thickness being the most important —_ + 
‘metric property. _ force of resistance to the beam flange » line load may be 4 


ey is a coe icient t depending on ti he width of column and beam flanges, ex- 
>) tent of two- -way bending, distance between column flange fillets, and boun- 


develop th the full yield of the be beam the column flange r 


Ag = = Oy tm + C1 te“ 


it the column resistance is found to be e deficient, ‘the column: is 
thin, and is effect is most serious in the region 


ing a plate bending mechanism which extends for a length of column flange 
equal: to 12 te: assumptions of the relative dimensions of beams and 


From this osedeniinis of assumptions, ‘the column mn flange thickness to. resist 


te = “es ee ee 
Results of tests showed that connections proportioned according to Eq. 8.26 7” 
carried the plastic moment of the beam satisfactorily. Fig. 8.15 shows the re- 2 
: "sults of tests on connections designed to meet the appropriate design criteria. 
‘The a tests are those without stiffeners which satisfied both Eqs. 8. 21 and 
If the column web tl satisfies Eq. 8.21, ‘column are 
“needed adjacent to the beam c compression flange. if the column flange thick- 
ness satisfies Eq. 8.26, column stiffeners are not needed adjacent to the econ " 
tension flange. If the column proportions: fail to satisfy either of these equations, © 
however, stiffeners should be provided. Two types of stiffeners are shown in 


8.13, the horizontal stiffener and the vertical stiffener. of 


= 


7 - Yield line theory is an upper bound plastic analysis of bending of plates in in which 
are ssoumed to form along lines in ina to form a mechanism. 


— 
— 
— 
— 
= 
4 
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» stiffener may be proportioned by considering the additional amount of resist- ; 


: ing force required to achieve equilibrium. The stiffeners should also be p pro- 


portioned to maintain stability under their full load. 
Horizontal stiffeners should preferably be furnished in symmetrical pairs 


eee both beam flanges. » They should be welded to the column flange and = 
ne, web by either groove or fillet welds. 7 Vertical stiffeners should also be fur- _ 
‘x nished in symmetrical pairs and should be deep enough toallow thebeam flange _ 


: force to be dispersedin the stiffener inthe same manner in which it is assumed 


‘to be dispersed in the web in Eq- 8.19. 


FIG. 8.15.—RESULTS OF Two- -WAY TESTS 
Results of tests on connections using both types of given in 
‘Ref. 8.11. Some of them are shown in Fig. 8.15. These tests and tests without — a4 
‘stiffening were used to help | establish ‘semi- empirical methods of deciding the 


| 

oslt]  \ BEAMS: 16 WE 36 | 

i 
% 
need ior stilieners and the proportioning eners to Tesist Of Deal 
flanges. If the moments applied by the beamsto the columns wereunsymmetri- 
cal, the same methods would be appropriate. 


The preceding analyses have shown that the use e of stiffeners may may depend 
either on compression in the column web (Eq. 8.21) or on tension normal * 


the column flange (Eq. 8.26) resulting from the end moments of the beams. 
For most of the rolled column sections normally only one type of failure need 
be checked, depending on the proportions of the column. _ The sections con- 
trolled Ay each type of failure may be tabulated, thus simplifying design. — 


4 


taking th/ ke = 


> 


8. 16 shows a - of the values rs of te/V keWe for all wide-flange columns 
in 1 the 8 in., 10 in., 12 in., and 14 in. - Series. It can be seen that in most cases zi 
the column need only be checked for compression. For values of te/ WkeWe KeWe be- 
Ss tween 0.91 and 0.96 the need for column stiffening will depend on the beam, and - 


q both and i compression sh should be ‘he 


When the the ot abt in the two beams at an interior connection differ by a 
large amount, they may cause large. shears in the column web, tending to de- i 
form the web in the same manner as in a corner connection (Fig. 8.2(e)). With 
such an unbalance of moments, the shear in the web should be checked, and if 
_ necessary, diagonal stiffeners or a doubler plate should be added (8. 12). In Fig. 
: 8.17(a) the shears and moments acting on a typical connection are ‘shown. Fig. 
8. 17(b) is a free body diagram of the forces acting on the column web just be- 


The forces are V, the horizontal shear present 


| 
— — 
‘From Eqs. need stiffeners in 
a y For all practical combinations of rolled beams and columns that might be —. 

tension region 4 
— 
— 
— 
| 
a taking tp /Kg ° t need stiffeners in 
\ 
— 
— 


é in the column above the connection ree tensile flange f forces T; and T.. 
‘The flange forces are obtained approximately by dividing the moment in the 
respective beams by the beam depths. _ The net result of these forces must ad -_ 
resisted by a shear stress T — on an area of column web equal to Wede. 


M2 M2 My (8.33) 
_ Assuming web to resist shear is a 


"stiffeners similar “to those in or web plates could be used 


_ to-column connection is obtained vd when M, equals zero. © 
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ates FIG. 8.16. —CRITICAL PARTS OF 
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MOMENTS ON BEAMS del 


FIG. 8. —FREE BODY DIAGRAM OF FORCES ACTING 
Py. abet, ON TOP STIFFENER AND COLUMN WEB 


= eee a four-way beam-to-column connection. Two beams may be framed 
into a column web either by direct welding or by ‘plate plus seat type of con- 
_— ‘nections. When connections such as these are used, the question has been 
_ raised as to whether the triaxial stresses in the column web might cause pre 
4 mature failure, or on the other hand, whether a beneficial effect might be ob- > 
3 tained through the partial stiffening action 1 of the beams fran framing into the column — 7 
web. 4 limited number of tests have been made on four- -way connections. 


— 
«| 


the specimens tested, there was no adverse of {triaxial stresses 


Sketches of the four '-way ‘connections tested are shown in Fig. 8.18. The 
direct-welded type, shown in Fig. 8.18(a) had a column web thickness meeting — a 

_ the requirement of Eq. 8.21 and hence the column web required no reinforce-_ 


ment to prevent crippling due to the forces a to the column flanges by oe 


Wel 

4 


Split-Tee Web Stiffener_ 
‘FIG. 8. 18.—TYPES OF FOUR-WAY INTERIOR 


bie web, appeared to stiffen the web. In the second type of connection (Fig. 8. 18(b)), oe 
5 the horizontal plates which served as top plate and seat plate for the beams . es 

framing into the column web also served as ‘stiffeners for the column. . Be- 
cause the beams framing into the web were not as deep as those welded to the . 
at column flanges, the bottom stiffener was 4 inches away from its ideal location 7 
se opposite the lower flange of the deeper beam. However, the connection per-_— ei 


formed connection in Fig. ‘8. used split- -tee 


| 


= 
— 
— 


e stems of the tees served to support both the | stiffeners ‘sand 


the column web. The flanges of the split tees served as suitable surfaces to 
which beams could be directly w welded to complete the connection. — 
tests of Ref. 8.11 were performed with considerable ‘axial load in the ro comer 
. ‘ts order that the results could reflect the influence of this practical loading 
; _ Situation. _ All of these tests were performed under a symmetrical moment cc con- Ga 
dition. Some special problems still remain to be solved when unequal moments 
For the design of four-way connections Ref. 8. 11 recommends that the con- 
| nection of the beams to the column flange be designed exactly as ina — 
as if the beams to the column web were not 


_ BEAM-TO-COLUMN CONNECTIONS (Beams Joined to Column a 
column stiffeners are not needed | adjacent anit the beam compression — 
Column stiffeners are ne t needed adjacent to the beam tensionflanges_ 


4 


Column Stiff 


‘proportioned to carry the excess of beam flange force over that 
| the column web and flange are able to — 7 (Stiffeners should also 

2 be ee so as not to buckle. 


When unbalanced moments or forces cause shear stresses 


the column. Such stiffening take e the form stiffeners or 


butt welds, the forces are carried or tension, and are 


“limited by the tensile or compressive resistance of the base metal or weld na 
- metal depending on which is least. Since weld metal is stronger than the base . . 
metal normally used in any joint, the strengthof butt welds should be calculat-— 
 edon the basis of oy, the yield strength of the base metal. rg 
_ A nominal shear stress of 22.4 ksi was selected for fillet welds. The basis 


13.6 ksi as presently (1960) permitted by the AWS code for building snl 
“tion Thus the calculated shear stress at load exceeded the 


; in excess of the column web capacity, stiffening should be provided for . A 
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ps se. stress. of 20 ksi. Use of these stresses at perio loadresults in weld 


ae for a simple beam whichare usually the same as those which have proven 
* ‘satisfactory for years in allowable stress design. 


The results of the tests in connections described Pon ere 6, and 8. 9 
verified that the welds designed using these assumptions we were re adequate in that 
“no o weld failure oc occurred. 


Welds m may | be assumed of on 
_ mum throat section the tensile yield stress of the base metal. seta. j 
Fillet Welds may be assumed to be capable of developing at | 
_ the shearing yield stress of the weld metal on their minimum throat: : 
; section. A safe design value may be obtained by multiplying the value : 
4 presently used in allowable stress design by the ratio o. /Ow, where ow m p 
_ is the the allowable tensile ile stress of the the material being used. yu anal ‘wi 
Bolted joints in plastic design will probably have their gi greatest 
= It could be argued that judicious selection of the locations of such joints cd 
_ place t them at points of low moment so that the joints will remain elastic. The lg 
Hover of these joints then could be handled by conventional elastic methods. i 
However, it sometimes becomes advantageous to make field connections ved 


‘Two possible approaches available for full capacity connections are (1) to 
-keepall parts of the joint “elastic” at maximum moment, and (2) to design parts 
of the joints so that fasteners and plates may yield at maximum moment along : 
; with the main members. Recent studies have emphasized the latter approach. 
A limited amount of research has been on bolted joints designed» 
a to develop the full plastic moment « of beams. References 8.15 and 8.16 give — 
the results of tests of connections using ASTM A-325 high tensile bolts. The | 
- connections in Ref. 8.15 included both beam splices and beam-to-column con- 7 
nections and for each case there were examples with the bolts in shear andin _ 
tannin. The connections in Ref. 8.16 were butt type connections with the bolts 
. i in tension. The results of the test programs show that it is possible to design - 
and fabricate high tensile bolted joints 3 which will develop the full bending strength - 


“factorily as reported 
___ Methods of analysis and design are not discussed in detail in Ref. 8.1 “ 
ae Methods are discussed in Ref. 8. 16, but the author suggests | that better pro- 
cedures based on the same fundamentals can be developed. cations paiil 
a _ The essential features of any analysis of a full-moment bolted connection — 
are the assumptions as to the load capacity of each individual bolt, the manner 4 
i in which the compressive components of force are re to be carried, the location Hy 
ae of the neutral axis, and the strength of any detail material (fittings) used in — 
forming the joint. It seems reasonable to conclude that precise solutions for 


4 


— 
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_ An example of the assumptions that might be made in designing a high strength . ai 
“4 bolted joint is given in Ref. 8.16 as exemplified by Fig. 8.19(a), 8.19(c), and — a 
o .19(f). The number and location of bolts in a joint were first assumed. For ey 
_ the bolts in tension, an “ultimate” load capacity was taken in this reference as a 
wie ~ 3 times the proof load of the bolts. 2 _ ‘The shear resistance provided by cach 
bolt at ultimate load due to friction was taken at one-fourth its initial tension. 
a structural steel parts of the joint were assumed to be stressed to yield 


esign, a compression area \ was computed sufficient 


Conventional Moment Splice fe) Seat ‘and Plate Type 


8. 1.19. ‘PES OF PLASTIC MOMENT CARRYING H H TENSILE 


tension side of the neutral axis (axis of rotation). The bending moment pro-- 
a F by the couple made up of the bolt tension forces and the — 72 
|: 


force, had to equal or exceed the applied bending moment. 
‘Reference 8.16 notes that plastically designed joints with high tensile ‘bolts 
in tension require a smaller number of bolts and less fitting material than con- =e 

ventional moment splices which use the bolts only inshear. 
Although accurate leading to the most economical safe 
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DESIGN 
be designed to. develop the plastic | moment 
_ This report is based on a series of preliminary reports weneat for he. = 
- Welding Research Council and the American Society of Civil Engineers by a 
i "research group at Fritz Engineering Laboratory, Department of Civil Engi- = 
as neering, Lehigh University. This group has been engaged through tk the Institute _ 
:, of Research in an n investigation « entitled “Welded Continuous Frames and Their — 
Components”; the staff members include L. S. Beodie, G.< + Detocell, Jr. ‘ 
a _ T. V. Galambos, R. L. Ketter, T. Kusuda, G. C. Lee, T. Lee, L. Ww. Lu, ~e 
and B. Thurlimann. The original draft was written by G.C. . Driscoll, 
h. Jr. William J. Eney is Director of Fritz Laboratory and Head of the Civil - 4 
‘The organizations | that supported the Research projects out of which the a 
Se report was prepared were the following: / American Institute of Steel 
Construction, American Iron and Steel Institute, the Department of the Navy | 
(Office of Naval Research, Bureau of Ships, Bureau of Yards and Docks), the — ; 
Welding Research Council, and the Column Research Council. 


; : Although much of the experimental and theoretical work reported \ was per - : 
4 _ formed at Lehigh University, the Welding Research Council and ASCE com- 
- -mittees have attempted to broaden this report by including the results of re- ot 


search at other institutions, both here and abroad. The inclusion ¢ of this 
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g CONNECTIONS FOR WELDED PORTAL FRAMES, ‘The 
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INVESTIGATION OF THE STRENGTH OF WELDED PORTAL 
; a _ ‘THE DESIGN OF KNEE JOINTS FOR RIGID STEEL FRAMES, BWRA iq 
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J. W. Fisher, G. C. Driscoll, Jr., F. W. Schutz, Jr. a 
BEHAVIOR OF WELDED CORNER CONNECTIONS, The Welding Jour - 
37(5), Research Supplement, 216-s, (May 1958) 
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COMMENTARY ON PLASTIC DESIGN IN STEEL: DEFLECTIONS. po 


Progress Report No. 7 of the Joint WRC- ASCE Committee o 


FOREWORD 
WORD 


"appears at the beginning of Seeman Report No.6 immediately preceding this 
Chapter 9, entitled “Deflections,” comprises tl the present paper. 


‘a 


Note.—Discussion open until September 1, 1960. Separate Discussions should be sub- is 

itted for the individual papers in this symposium. To extend the closing date one fe 
onth, a written request must be filed with the Executive Secretary, ASCE. This paper | 
- part of the copyrighted Journal of the ry Mechanics Division, ee 


— 
— 
— 
— 
 -‘This pape is the seventn in a series of reports on Fiastic Design emanat- a 
pee ‘—l ing from a Joint Committee of the Welding Research Council and the American he = 
Society of Civil Engineers. Reference should be made to Progress Report No. 1 
se _ _ (which appeared in the July 1959 issue of the Journal of the Engineering Me- a 
| = + chanics Division of ASCE) for such items as the Foreword, Synopsis, Deli a 
— ; 7 _ the strong recommendation that it be read carefully byall structuralengineers. [im 
— 

— 
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The deflection may y exceed established limits set up up by a 
building code or specification. __ Deflections may hinder the operation of mov- 
ing parts or the closing of doors, or they may cause the cracking of plaster 
in finished ceilings. Thus, it is desirable to have methods for calculating the 
_ deflections of structures in both the elastic and plastic ranges. ae es 
chapter will outline the basic problem of structural deflections. 
discuss briefly the calculation of deflections in the elastic range. 

- ealeulation of deflections at ultimate load will be coveredin more detail. Illus- 

_ trative examples will show that the deflections at working load of ‘structures 
designed plastically usually | are no greater than those of structures designed 

tom calculation of the deflections of structures in 1 the elastic rang range is ;cum- ay 
_ bersome even after the strength analysis is complete. For this reason, it is ; 
os procedure to dispense | with the calculation of deflections w when suf- 
experience has been gained with a particular type of structure to indi- 
cate that deflections will not be excessive. In plastic designthe same practice | Bt : 
‘Structural ‘deflections are 


riate conditions, of the differential equation: 


_ = curvature caused by sities ‘Thus it is a function of moment which 
= in turn i isa a function of x. In the elastic range M/E IL real 


There are of solving this in both the elastic and 
- plastic ranges. The methods differ principally in the means of expressing o 
and in the methods of handling In the sections, 


- _ Methods for calculating deflections inthe elasticr range are eavailable in 7 
- on structural theory so these methods will not be explained here. To makea ~~ 
BS precise calculation of deflections at working loadof structures proportioned ae 


4 

= plastic design it would be to use one of these methods. will be 
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PLASTIC DESIGN 


moment of inertia of the cross section of the member. 
feature of of the elastic methods is that eupressions giving 
ee the slope at all points of a deflected member are continuous functions, and that mi 


with other members may be fully rigid or completely free. 
rs 


Certain of the methods also are adaptable for the consideration of facto 
a8 ‘cuaileding to deflection which ordinarily are neglected in the study of rigid _ 


= frame structures. These factors : are the deflections caused by axial force and 


: Several methods for the in plastic 
= are available. Five of these methods are outlined, discussed, and com- a 


pared with ‘experimental results in Ref. 9.1. _ In these methods, the moment- ; 
_ curvature relationship derived from the elastic and plastic stress-strain re- ‘= 
oe lationship is substituted in Eq. 9.1 and, using appropriate boundary conditions, o 
Slopes and deflections are obtained by integration. The methods differ only 


simplify the computation. Table 9.1 outlines the features of the five 

methods. Fig. 9.1 shows the stress distribution associated with the assumptions — 
_ for each of the methods. _ _ The simplest approach, the “plastic hinge” method, 

is most suitable for structural engineering applications. The additional work a 


involved in the theoretically more “exact” methods is not by the accu- 
: 
4 


= a The “plastic hinge” method assumes that the: structure has only ait re- y 


ae gions and localized plastic hinges. This allows the use of the curvature equation - 
3 in the form given as Eq. 9. 2. The boundary conditions are then changed to al- 
low for freedom of rotation at plastic hinges. . As pointed out by Symonds and 
i= 
Neal, the controlling boundary condition is continuity at the last plastic — BS. 
fi to form (9.2).* Any of the well-known methods of calculation of deflection in sf 


the elastic range may then be used with suitable handling of the boundary ¢ — 


. ‘ditions. Methods which have shown n promise ar are: “slope deflection,” ¥ “moment A 


area,” “conjugate-beam,” and the “method of work.” Some advantages « of the 
plastic hinge” method are e that it obviates the necessity of tae 


* An exception occurs ohin hinge rotation has camenee at one or more sections where _ 
a the bending moment at ultimate load is less than the fully plastic value. Fora discus - _ 
= 


later, an upper bound to this deflection may be determined without 
sorting to the elastic analysis of the indeterminate structure.) 
One thing elastic methods have in 
* al 
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1— METHODS FOR CALCULATING (G DEFLECTIONS 
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» STRESS DISTRIBUTIONS ASSUMED FOR CA LCULATING DEFLECTIONS 
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PLASTIC DESIGN 


“use of the deflection formulas tabulated in structural handbooks 
_ The other four methods discussed in Ref. 9.1 are “simple plastic theory,” 
Area,” “mathematical integration,” and “numerical integration.” These 
methods utilize more accurate representations of the moment- curvature re- 
lationship and properties of the steel to obtain a small i increase in accuracy of 


expense of. a considerable. increase in n complexity. trial-and-« -error r process 


_ must be used not only to determine the boundary conditions (last plastic hinge) 


put also to determine the moment and curvature. For this reason, these methods 


7 Wh, equation is used in this form, E I must be constant in segment AB.) 


: have only been used in the calculation of deflections of very simple structures 
such as continuous beams and rectangular portal frames (9.1, 9.4). For more 
complicated structures, the use of a digital computer would be required. 

_ _Deflections of more complicated structures have been calculated by the > 

“plastic hinge” method with reasonable agreement between calculations and 

tests. For most practical cases, therefore, the “plastic hinge” method i is recom- 

mended for deflection calculations (9. 2,9. 


In this section, the approximate dette: at t ultimate load of the simple 
portal frame shown in Fig. 9.2(a) will be calculated. The assumptions used will 


be those given in ‘Art. 9.3 for the “plastic hinge” method. One method of cal- — b 
is to the slope-deflection -2,9. 3, 8. 7,9.9,9. 


See Fig. 9.3 Clockwise M and are positive. 
The mechanism and moment diagram at t ultimate load are given in Fig. 9. 
and 9. 2(c), respectively. In Fig. 9.2(d) are given free body diagrams of the» 
‘portions of the frame. ‘To solve for the deflections 6y and 5y, two boundary _ : 
_ conditions a are needed. One of these is that 69; equals 693 because of continuity — 


at joint 2. The second boundary condition depends on whether the last plastic 


- hinge forms at section 3 rales 4. _ A trial solution of the problem will be 4 


vA 


— 
— | 
— | 
— 
7 
. 
— 
4 
P 
q 


gives one for bya and applies for either r of 
"Trial at section 3, at section (032 = = 


Considering segment - - 


Considering s segment 3 - 


In the above the ‘deflection continuity assumed 


9.6) 


iim 
4) 
a 
i 
ft 
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_ These values for by: and 5H would be the correct values for the dattnitiiniiis Ss 


“Trial at section 4, = 


<0 


and and solving with Eq. 9.4 one 


~ 33 
The deflections ‘calculated assuming at s A are 
_ those obtained on the basis of continuity at section 3. Consideration of the ro- - 
possible “reverse “kink” would occur at ‘section 4. Because of this “reverse 
icin,” the assumption resulting in such a solution is discarded. _ Therefore 7 
oe thee deflections s assuming continuity at section 4 are the correct deflections at 
: bi ultimate load. It has been established elsewhere (9.2) that the correct set of 
: x calculated deflections will be the largest set from any series of calculations 
based on different assumptions of the last plastic hinge. The last plastic hinge 
to form will = the he hinge at at which hcontinuity i is sassumed in rs the e larg- 


as example ‘of ee calculation of deflection at ‘ultimate load will be 
_ given using the “dummy load” method, variations of which are called the “unit , 
_ load” or “virtual work” method. _ This momee is used in many text books on 
__ The operations of the method are as follows: " first, the moments M due to | 
‘the applied loads are } ENS, second, a ———s unit t load is ; applied to the 


— 
— 
= 
— — 
— 
— — 0- £72 * \Mp - = ™ 
— 
— 
— 
| 
pee: 


unloaded structure at the point of the in the 

_ the desired deflection. A set of moments m in equilibrium with the unit load 


deflections at ultimate load are: 


(1) All members are assumed to be fully elastic between plastic hinges. 
Continuity is assumed at the last plastic hinge. 


be (3) The dummy load is applied to an auxiliary ‘structure with frictionless. 


hinges assumed at all but the last plastic hinge. 
le If the location of the last plastic hinge is unkown, a solution must be > 


4 | 
— 
— 
7 


_ After applying these conditions the deflection at the unit load one to the origi 4 
] ae loads is determined from a virtual work equation of the form ae 


a Deflection in direction of unit load. a ee 


‘Moments in structure due to applied loads. 
= Distance along 
= Young’s modulus of ‘elasticity. — 
= Moment inertia of member. 


The equation is integrated over the whole structure, , using intervals where al 


‘FIG. 9.3. FOR SLOPE DEFLECTION EQUATION. 
centrated loads and hinges, the discontinuities in the curvature functions are 
_-‘The vertical deflection at ultimate load of the ridge of the gabled frame in “a 
= Fig. 9.4(a) will be calculated. - The locations of the plastic hinges and the “M” 
= moment diagram due to the applied loads are shown in n Fig. (9. -Afa). ‘The neces- a 
additional from the analysis a: are: 
= 16.84 Mp 


if 


| 3 
m with unit load. — 
— 
— 


‘rafter. . By making two deflection calculations, one assuming ; each hinge in turn 
j to be the last, the correct deflection will be determined. = » 
‘In Fig. 9.4(b), a unit load is applied at the ridge of. the strt structure witha 

i tionless hinge as assumed at the knee. . Ih other words, there is continuity in ‘en 


Applied 


Last tangs Assumed in Rafter | Last t Hinge Assumed at Knee: 


om. 


with Hinge i in 


‘ne 9.4,—LOADS AND MOMENTS FOR DUMMY LOAD SOLUTION OF DEFLECTION 
the windward rafter in which the last plastic hinge | is assumed to occur. The 
resulting reactions due to the unit load and the “ m” moments are also shown. 


‘Moments M and m are in the proper ds: 


> 


— 
1960 
q 
4 
(b) Unit Load and 
7 
= 
— 
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Instead of a at the knee, the strectere with the dummy 1 load in 
‘Fig. 9. 4(c) has a hinge in the windward rafter. The structure in this condition — 
2 behaves as a three- -hinged arch with the “m” moments and reactions ‘shown. 


Using these “ m” "moments with the “M” ‘moments of 9. a second ex- 


zs 


(3 


ntegrating and substituting f for results in 


im 
9 The conclusion that the largest calculated deflection based on a given ultimate | a 
i moment diagram is the correct deflection is explained on pg. 192 of Ref. 9.3. 
_ It is stated that “any incorrect solution can be found from the correct solution a 
by superposing the displacements and hinge rotations to a backwards: motion 
the collapse mechanism.” 


9. 5 STEP-BY-STEP CALCULATIONS OF DEFLECTIONS. 
_ Occasionally, and especially for research purposes, ‘it becomes desirable 


construct the load- deflection (or rotation) curve fora structure. >» 


was ated deflections is the correct one, 
| 
a 
— | 


 - done by means “ts a step-by- step ete ‘because the load-deflection be- 
ri havior is approximately linear between the formation of "successive plastic 
hinges (9. 12), _ The calculations to be discussed here will be based vad on the 

r _ As an example, consider the deflection under load of the built-in beam in 


¢ 
Fig. 9.5(a). In the elastic range, the load-deflection curve follows the line 0-a. = 


+The curve will cease to follow line 0-a when the moment at the left support — 
_ Mp. _ This point is designated by A on the solid curve. Thereafter the 
= at the left ‘Support ¢ cannot increase, but the moments can increase at 


every other point in the beam. _ The subsequent increment of deflection of the 
sa _ beam will then be the same as the propped — of Fig. 9. “ss The elastic 
3 
‘6 


) 


DEFLECTION UNDER 


EP 


po of load corresponding to AB is iisranlnst by considering the — 


=. under the load. A certain moment corresponding to point A exists under the = 


load in the original structure (a). The difference between this moment and Mp 
is made up by the moment under the load of auxiliary structure (b). The _ 
_ crementof load is the amount required to raise the moment of auxiliary struc- — 
ture (b) from zero to the amount of the original difference from Mp, and the 
increment of deflection of the auxiliary y structure (b) caused by this incren 


a this stage, ‘mei are plastic hinges both at the load and at the left sup- 


port. Consequently, no increase in moment can occur at either of these points. 
auxiliary structure of Fig. 9.5(c) will be used to determine the next incre-_ 


ment of load and deflection in the same manner as was done for the first two 


_ 
j 
3 
| 
1 Structure Ww iginal beam 1s de 
letlecUlon Or flection of the or 
line O-b. 3 
— 
— 
— 


. increments . The. load at the end of the third step will be. based o d on the ‘the > moment — 

the right support as summed up for the three cases. 
_ The final stage of the load- -deflection curve is the borinontal eta cD 


7 matically by s structure (d) which has three hinges. . Neglecting catenary ef effects, 
the load- deflection structure (d) is represented by baseline d. 


«FIG, 9.6.— 

In summary, the load- deflection curve a structure loaded | to 

may be obtained by superimposing on the elastic load- -deflection curve of the 
primary structure, portions of the elastic load-deflection curves of auxiliary - 

Structures. The number of auxiliary structures will be the hager as the number 

Shy _ The same procedures may be followed for rigid frame structures. Fig. 9.6 

_ shows the load-deflection curve of a gabled portal frame and the load-deflection 
4) _ curves of the necessary auxiliary structures (9.7). The deflection at point (4) 4 


ei in n Fig. 9.6 is the same as would be obtained by the methods of Art. 9. 94, 


— 
7 
oe a 4 


ae 9.6 APPROXIMATE DEFLECTION AT WORKING LOAD 


‘It has been suggested (9. 13) that an upper bound for the deflect 
ing load may be obtained by dividing | the calculated deflection at ultimate load — ci 
we by the load factor. Calculations were made to compare this type of estimate > 7 a 
_ with calculated and measured deflections at working load. Bar graphs showing —e 
these « comparisons for five structures are. given in Fig. 9.7 . Each ¢ — . ] 
_ is plottedas a percentage of the calculated deflection 5y at ultimate load — - 
the deflection estimated at working load is always 54% of 6, whenthe load fac- 
tor is 1. 85 . The deflection calculated from an elastic analysis. at the plastic 
design working load and an observed deflection at working load from a a testare a 
~ shown where available. For a further comparison, a bar showing 1/360 of the 3 
Span is to same scale as the deflections.’ 


Y= 
| 


th 
Determinate fonts 6 
‘FIG. 9. 1. |. COMPARISON OF DEFLECTIONS AT WORKING AND ULTIMATE LOADS 
_ WITH ESTIMATED WORKING LOAD DEFLECTION. 


The bar graphs show that the estimate of deflections made by dividing 


th the load factor is more > and | more in error as the number of statically — - 


approximate estimate is very 

a good for one redundant, the estimate is four times the actual deflection at work- 
a ing load for the structure with six redundants. This estimate would imply ex- “a 
cessive deflection whereas the actual deflection here was well below L/360. 


cal 


te. L/360 is not suggested as an appropriate deflection limitation; it is shown merely 


- because it is ities of the same order of magnitude as any aaa: deflection limi- 

6 


= 
— 
y 
ds on the order of 
ight be. Much depe c hinges 
error the origina ic hinges. In true deflection. 
pe tion of the will equal the 
simultane — 
— 


edure suggested is entiatactory whenever esti- = 


One of the most important questions about deflections in plastic 
"whether or not they will constitute an undue design 1: limitation. The 

= this q question | on the basis of available test results and calculations is: no. 


~ _ Admittedly, on the basis of the few tests and calculations for examples — 


fe Variations in = and do change the relationships. However, it is 


‘ 


G om 14 ve 
FIG. 8. —LOAD- TION RELATIONSHIP “THREE DIFFERENT 


: deel BEAM DESIGN FOR SUPPORTING THE SAME LOAD. 


_ Consider the case shown in Fig. 9.8 (9.10). The problem is todesign a beam 


to carry a given load on a given span. First, picture the beam being designed 


and the load-deflection curve would be as shown by curve I in Fig. 9.8. A hori- _ 


in plastic The rolled section required would be an 
zontal mark on the graph indicates the value of the working load, and a oe 
if 


eal mark indicates the deflection which is equivalent to 1/360 of the span. _ 
‘Next consider a fixed-ended beam designed to carry the same load on . the 
Same span according to elastic design. Because of the negative moment at the 


Z.. fixed ends, only a a 16WF36 would be required. The load-deflection curve would 
fe as given by curve II in Fig. 9.8. It is seen that the deflection at working load 


mated deflection calculate = 
— 

— 
— — 
— 
— 
q 
— — 
— 
— 
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April, 1960. 
fixed ends would support a much ultimate lend than 


that of the simple beam. This indicates that, although both beams are designed 2. aa 


car case Ta I and be 

| 

9. 8 ROTATION R REQUIREMENTS 


we In addition to the problem of in the at 
ms plastic hinges must also be considered. For example, in the beam of Fig. 9.5, hs 
* F it is the inelastic rotation at the left support which allows redistribution of © 
Be = - moments to the load point and right support of the member, thereby wll 
or possible the increase in load from points A to D. Similarly, in any ee 


except in the case where all hinges form simultaneously. The measure of 
required rotation is referred to as the hinge angle - It is of interest to deter- 
mine the hinge angle that would be required in structures in order to — 
ertain secondary design guides which depend upon this — = 2 a 
_ The hinge angle may be calculated as part of the deflection calculations de- = 
 seribed in Art. 9.4. The results of numerous hinge angle calculations are aa 
a < _ sented in Ref. 9.5. Some of these results have been used in the determination — 
ie, of the required lateral bracing for bending 1 members given in Art. 6. 4. ie 
pe In some cases, the calculated hinge angle required to form a ‘nechaniam 
may be quite large, implying a limitation on the design. However, practical = 
considerations suggest that almost as good load carrying capacity may exist — ee 
with a much lower hinge angle requirement. Fig. 9.9 is plotted the theo- 
"retical load versus hinge rotation H at the first plastic hinge of a two-span eg 
a ae ‘ame. Seven labelled points on the curve indicate the formationof seven plas- 
“tie hinges in reaching tl the mechanism state. * Point (7) is the hinge angle that 
would be calculated in a rotation capacity study. Note, however, that at the 4 
_ formation of the next-to-last plastic hinge, point (6), the load is 98.2%of the _ a : 


maximum load. In order to attain this load, the required rotation is only 0.54 — ae 


L/(E D, or about one-third of the amount calculated for a full mechanism 


a 52 Mp L/(EI)). Assurance of reaching 98% of a predicted load, along with uk 
a reasonable load factor or safety factor is generally accepted as good met 
} neering practice. . Thus, it is seen that the problem of large required ee + 
y oe highly redundant structures is not as serious as might first be suspected. = at 
cnet Experimental verification that a structure may behave quite satisfactorily a ” 
even though it it may not meet the full rotation aa is given in Ref. 9. 6. Ba ‘ 
a did not quite attain the computed ultimate onthe basis of a "40 ksi couponstress, 
hakss it reached 97.4% of this value, or well in excess of anultimate load based upon — 
a 33 ksi yield ‘point (Fig. 9.10) . Observations indicated that six of the seven - 
plastic hinges were formed before occurred due to the insta- 


—— 
— 
third case in Fig. 9.8 indicates an even more economical solution by 
design. Here only a 14WF30 is required. The load-deflection curve 
aA 7 shows that the load capacity is equal to that of the simple beam. Though the as Bae <b 
ay - ay 
Ly 4 
= 
xg 
4 
‘There is every reason to believe that with adequate lateral support the seventh 
a Le plastic hinge would have formed, and an even greater load would have been i Yr 


a , reached due to strain- n-hardening. For highly indeterminate ‘multispan frames, 

large hinge angles m may be required for the full development of a mechanism. a 
However, for these same frames the load at the formation of the next-to- reed e- 

3 hinge is almost as great as the maximum load, while the rotation at that in- a 

stant is substantially less, as in the case shown in Fig. 9. 9. _ Hence a lower | 


— 


IG. 9.9.-HINGE ANGLE FOR A TWO-SPAN PORTAL FRAME. 
a limit for rotation capacity requirement, suitable for general design purposes, a 
need not be to theoretical r rotations 1s associated | wen 


be. 
| 


— 


_ Welding Research Council and the American Society of Civil Engineers by a oe 
research group at Fritz Laboratory, Department of Civil 


: fi Components” 5 the staff members include L. S. Beedle, G. C. Driscoll, ad. hl 
_ V. Galambos, R. L. Ketter, T. Kusuda, G. C. Lee, T. Lee, L. W. Lu, A. Osta-— 
- penko and B. Thurlimann. The original draft was written by G. C. Driscoll, Jr. 


i= 


= 40ksi) 76.3" 


dinkips 


Loa 


tical L 


(Elastic Design) = 


Vertical D Deflection of Left Ridge in 

FIG. 9.10, EXPERIMENTAL LOAD-DEF LECTION E OF TWO- 


as 


‘This report is on a series of preliminary for 


William J. Eney is Director of Fritz ree and Head of the Civil Engi-— 


organizations that supported the projects of which the 


basic report was prepared were the following: American Institute of Steel Con- 
‘struction, American Ironand Steel Institute, the Department of the Navy (Office 


of Naval Research, Bureau of Ships, Bureau of Yards: and Docks), the Welding ; 


Research Council, and the Column Research Council. 
Although much of the experimental and theoretical work reported was per- | 
formed at Lehigh University, the Welding Research Council and ASCE com - 

mittees have attempted to broaden this report by including the : results | of re- 


mation and of certain unpublished data is ‘gratefully acknowledged. Semel 


— 
— 

| 

OS 

| 
= 

— 


DEFLECTIONS 
“ADDITIONAL NOMENC 


(See Progress Report No. 1 for basic nomenclature) ja 


= moments in structure in equilibrium with a unit load a. 
= distance along a member 
non-dimensional parameter defining t the distance to the plastic hinge 
inthe rafter ofastructure 


deflection of a joint due to translation 


OA a = rotation of joint A due to transverse loads assuming simply - sain 
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their work. _ They wish to apologize for overlooking the fact that Reference 9 "| 
- had already appeared in published form at the time this report was printed. 
; Further correlation between the authors’ interaction | curves (Figs. 6 and 7) 
a 4 and the “C. R.C. Interaction Equation,” has been recently furnished by nae 
Massonnet.3 He has shown, in Figs. 21 and 22 of his paper, that 
correlation exists between the writers’ “exact” interaction curves and the 
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BLE ELECTRICAL A ANALOG. FOR yon YIELD 


Discussions by Glen V. Berg and Spiro S. Thomaides, T. K. Caughey, 
G. W. Housner and D. E. Hudson, and O. A. -A.Glogau 
GLEN V. BERG,} M. ASCE, and SPIRO S. THOMAIDES. 2—The authors de- 
serve commendation for having explored one of the lesser - -known facets of the 
earthquake-engineering problem. In so doing they have contributed a little 7 
more to the fund of knowledge which ultimately leads to better building codes _ 
and more competent procedures for aseismic design. 
_ General Discussion.—The paper suggests that a structure which survives “e in 
zi strong-motion earthquake in which it is distorted to such an extent that the a ' 


_ masonry walls crack, but otherwise appears sound, has worked through a yield 
7 cycle. This may betrue. However, for a steel-framed building of ordinary _ 
proportions, the distortion at which the frame members would yield is of the Sw 


ae of ten times the distortion at v which masonry cracks would | appear. The 


yield-cycle concept may be useful even for such cases. Experiments 
‘have been reporteds»4 in which concrete filler panels in reinforced concrete 
_ frames have been loaded well into the inelastic region and then reverse loaded 7 a 
the negative inelastic region. The load-deflection curves show a strong 
hysteresis loop, even for an unreinforced panel. A response curve of the elasto- z= 
wall tests conducted at Stanford University,> as well as other experiments,* 3 
show response curves close tothe elasto-plastic type for one-directional load- 
ing and unloading. ‘It has been suggestedS that reverse loadings of complete 
structures might produce load-deflection curves that mirror the curves for 
_ the first loading and unloading cycle, perhaps with softened features and a de- ie 
cing ultimate strength. If this is so, the procedures used here may produce — 
—— results for a larger class of structures than the authors claim. — 
_ The writers feel that yield excursion is the most important of the response 
| + tee shown because it is a measure of the energy dissipated in plastic — 
7 defor mation. _ Any building code that permits plastic deformation must, in our > 


type would represent a rough approximation to this behavior. Shear- 


8 October, 1959, by G. N. Bycroft, M. J. Murphy, and K. J. Brown. ae 
Assoc. Prof. of Civ. Engrg., Univ. of Michigan, Ann Arbor, Mich. 
2 Graduate Student, Dept. of Civ. Engrg., Univ. of Michigan, Ann Arbor, “Mich. abe : 
3 ’ “Frictional Effects in Composite Structures Subjected to Earthquake Vibrations,” 
Lydik S. Jacobsen, Stanford Univ., Calif., March,1959,. 
oP 4 «Studies of Shearing Resistance “of Reinforced Concrete Plates,” by M. Tomii, Re- 
_ port of the Inst. of Industrial Science, Univ. of Tokyo, January, 1957, (in Japanese.) ei 
ge “The Behavior of One-Story Brick Shear Walls,” by Jack R. Benjamin and Paes “ie ee 
Williams, ASCE, No. ST 4, July, 1958. 
Discussion of “Structural Dynamics in Earthquake Resistant Design,” by 
Jacobsen, ‘ASCE, No. ST 1, January, 18 1959. 
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Ther process susedis a -degree simplification of usinga 
i "Kutta technique. To permit comparison of results, digital analyses have i 
of the same systems reported inthe paper, responding to the same earth-— 
quake component. The results are shownin Figs. Bl through B7, in which both © 
= discussers’ digital results and the authors’ analog results are shown. oo . 
= curves for A = 0.10 have been omitted for clarity. A curve-by-curve com- 
parison discloses several points worthy of comment. 
.. sa ‘The digital spectra are not as detailed as the analog spectra, so the m minor 
© fluctuations are not strictly comparable. . However, the severe local fluctuations a 
that appear in the analog curves of Fig. B1 in the period range of 3.3 to 5. sec 
_ are absent in the digital curves, and the points for the digital curves, which _ 
were calculated at 0.1-sec intervals of period in this region, give no hint that — 
local exist. Perhaps the authors can suggest a reason for 
The yield- displacement curves, ies. B1 through BA, | show a qualitative a 
: of results, but the digital curves are lower than the analog curves © +e) 
os for low values of yield level, and higher for high values of yield level. For — ro 
 G=0.175 and G = 0.24 the digital maxima are _ roughly twice the corresponding ol 
The yield-excursion curves, Figs. B5 through B7, also show a qualitative 
_ similarity, but the digital curves lie well above the analog curves. The digital a: 
- maxima are roughly twice to three times the analog maxima, and the discrep- | - 
ancy tends to increase as the yield level increases. The disagreement here is 
especially disturbing because yield excursion is a measure of energy 
‘pation, which is an especially important parameter in limit design. 
Possible Reasons for Differences.—There are four possible causes for 
<= discrepancies ; namely, differences in input data, errors in digital computation, © 2 
cae errors in analog computation, or an inherent ill-condition of theproblem itself. 
aed It is presumed that the input data for the analog was obtained by photo-. “a 
_ graphic enlargement of the accelerogram record, and it seems reasonable to 
assume that the error here is small. Additional information on this point would | 
“aed = produce the digital input data, the writers approximated the first 30 sec 
of the accelerogram closely by a piecewise linear function of time, and put the — 
_ time- acceleration coordinates of the intersection points of successive anda 


Traneating the digital -accelerogram at 30 not | 

results. | The ground acceleration never exceeds 0.05 g after 27 sec, as com- — 
pared with a peak of 0.33 g in the first few seconds. By 30 sec the oscillation — a 
‘a would be damped sufficiently (the undamped case was not considered) to pre- 


_ One might question whether the piecewise linear approximation used for the 7 ; 
“digital was accurate. But the punched- 
“Dynamic Analysis of Elasto- -Plastic Structures,” Glen V. Berg and Donald 


Da Deppo, to be published in Proceedings, ASCE. = = 
“Punched Card Accelerograms of Strong-Motion Earthquakes,” by GlenvV. and 
Univ. of ‘Michigan Research Inst., 4958. 
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Moreover, velocity and curves” obtained 
_ integrating the punched card dataS agree well with previously published 

sults. 9 Linear response spectra have also been calculated and agree well with 
ai 


i. published spectra. It, therefore, appears that the secitil data should not be a | 

errors in digital computation are small. In using elementary fod 

shapes for which exact solutions are known, the digitally computed re- 
— sults ‘were accurate to four significant decimal digits. _ During an — 
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N, INCHES 
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a 
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06 
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ee 


TOTAL YIELD EXCURSION G = 


solution, the computer evaluates the energy | input, , strain cape 


a kinetic energy, energy dissipated in viscous damping, and energy dissipated 
= plastic deformation, all as functions of time. Different computation schemes 
: are used for the different energy components, so that an error generated in one | 
a component would not induce a compensating error in another. At the end of the | 
* solution, the energy input less the sum of the energy components is printed, — 
to Serve as ‘as ar an indicator of the > accuracy of t of the solution. The The total skould, of aa 
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“Lateral of Barthquake and Joint Committee Report, Transactions 
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course, be zero. Inn most. cases, vit wel substantially less than 1% ¢ of a 
SS the energy input. The largest discrepancy found to date has been less than , =. 
3 


ited If the problem were inherently ill-conditioned, small differences in the in 
_ put could produce large differences in the response. Although the non- -linearities 


_ increases of 7% to 22% in total energy input, as compared with 1% for a linear 


excursion. In the second perturbation the odd- oumbered acceleration ordinates 
were increased 0. 01 the even- i-numbered ones were decreased 01 


77 ‘The effect of this upon the response was practically nil. In view of these re- 
‘sults , it seems that the problem is not ill-conditioned. = a 
- 7 _ 8 Error control for the analog was not explained in detail in the paper. If < 
aa - error were of the same order as that encountered in the square wave _ F 
_ check (Fig. 4 of the paper) the results would be adequate. It is seenin Fig.4 
thal that t the analog response is consistently below the theoretical response for ‘the 
‘square wave, which produces only one yield. How does this pr project to the case 
_ where yield occurs many times in both directions? What checks were made to 
insure that the yield displacement and yield excursion circuits were accurate - 
oly recording all yields in such cases? Further information on the control of. a 
vided new researcharea is being odd chance has made it possible 


to analyze the same problem by two entirely different methods. _ The ae 


have an encouraging aga similarity but, contrary to ‘expectations, they 
onal 


. CAUGHEY, 10 G. w. HOUSNER,11 1. ASCE, and D. E. HUDSON.1 On.12_ 
ue generally recognized that the usual earthquake response spectrum 
curves, which have been obtained on the assumption that structures remained 

_ linear, are far from complete descriptions of the actual state of affairs. 
large energy dissipations which must come into play to limit the responses to 
_ large earthquakes must be associated with non-linear plastic action in the struc- 
= — itself. The problem discussed by the authors is, thus, a most important 
one for earthquake-resistant design, and the extension of — -spectrum 

techniques into the yield range is a major step forward. 

_ Although some calculations have been made of particular yield situations — 

by digital computer techniques, the present investigation is the first complete 

application of analog methods to the problem. Some very preliminary = | 

“ments made at the California Institute of Technology, with analog circuits, 

_ indicated that the problems of obtaining force-displacement relationships with 

accurately defined sharp breaks was a difficult one. _ By using silicon diode — 
limiter circuits the desired characteristics could be obtained over a limited 


piv, of Engrg., California Inst. of Tech., Pasadena, C 
au Div. of Engrg., California Inst. of Tech., Pasadena, Calif. — 


* Div. of Engrg., California Inst. of Tech., Pasadena, Calif. 


pede’ » - intuitively, that the near-random nature of the input would preclude any ex- 
treme magnification. As a rough check on this, several points for the case 
seo. _ G= 0.175 were re-run using two perturbations of the accelerogram. Inthe _ a 
— 
— 
a | 
— 
— 
— 
a 
— 


range > of The p of the variation of analog 

—= frequency has | been minimized by the scheme described in the present — 
=a paper, since it is the exciting frequency that is changed rather than the analog 
_ frequency itself. It would be of interest to see the measured force- a 


uf poses, the exact shape of the knee of this force- -displacement function is per-— 
_ haps of only secondary importance, in view of the fact that such relationships —_ 
" are usually not known for actual structures. The details of this function would 
cae be quite different for various types of structures, and it will require - 
| experimental studies of full-scale structures to definitely establish these pat- hd 
terns of behavior. It is because of these uncertainties that the linear assumption 
-— has an important place in design calculations, although it must always be a 
= in mind that little or no information about ultimate 


: be involved. The qualitative significance of the results should not, how-- 
ever, be affected by these basic difficulties of analog technique. 
ie is evident from the conclusions reached that such yield phenomena will 
be of the greatest importance for earthquake resistant design problems. itis : 
_ to be hoped that the authors will be able to continue these studies along the al 


A. GLOGAU.**—The authors are to be on their excellent 


work of bringing laboratory methods a step nearer to to actual conditions. . We ; 
_ must advance further along this line before analog results can be used directly 


for design code suggestions such as those in this paper. As an encouragement " 


to further research the following comments are made: 


very accurately determined, since ‘small differences of large quantities 5 


_ _In the following, a rigid building is meant to be a modern reinforced con- 


crete shear wall structure and a flexible to be a 


Earthquake resistant design is undoubtedly a prob 
5 aus On the basis of some earthquake spectra, and a number of assumptions, it 
has been shownthat under certain foundation conditions rigid buildings develop 
larger base shears than flexible buildings. = 


As a resulta ‘suggestion is made that rigid buildings should be designed for 


- a Rigid buildings when designed for only 10% ¢ g have an excellent record in 
7 a ‘severe earthquakes and would, if this suggestion was put into codes, be placed 
; a — atasevere disadvantage co compared to the flexible ones, about whose performance 
q aren It is, therefore, well to remember that there are a number of questions, 
with great bearing on the onto, which have not been peesmmaneuatiiead answered 


— relationship and natural period. researchers 


a - are of the opinion that they are inversely proportional and, thus, cancelling the 
advantage of the long period building (7410). 


References in parentheses to the 1956 Barthquake Conference. 
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. The effect of various types of materials on earthquake s epectra, 
_ 3. Probability of resonance: Is the flexible (steel frame) building more — 
—_ likely to become “fully” resonant in an earthquake having a coherent train of F<) 
_ periodic motion, than a rigid building with its ts large a in ——— 
The percentage of the inherent strength of a shear wall building 
nine to the normally uncalculated effect on damping and energy absorption of 
wa partitions, stairs, , etc. In the logically designed and constructed flexible balla air) 
this effect is difficult and expensive to achieve. 
_ 5. Collapse resistance to ground distortion: Rigid buildings, in alee : 
the box-type construction, have much greater inherent collapse resistance to a 
ground motions, in particular tilt, than flexible buildings. (The Polytechnical — 
- Schoo in Mexico City, appears to have been a dynamic design (flexible) but col- 


lapsed in 1957, probably due to foundation failure). af” 

Rigid buildings have been observed to rock without the existence of appreci- a 

"able horizontal stress -3). atthe actual ase sear 


_ Summary.—Apart from the uncertainty as to. what the actual base shear 

- amounts to, the uncalculable effects of item (5) above have as great a bearing 
_ the choice of a suitable load factor (and hence equivalent acceleration wood 
in design) as the ability of a building to absorb energy. It is good engineering — 

" practice to use relatively high load factorsfor new and little proved construc- 
tions, and to reduce load for those with 
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“THE NEUTRAL AXIS IN PLASTIC BENDING OF BEAMS2 

Hite. Discussions by John C. Sprague and H. Norman orman Abramson a ae 

C. F. ASCE. —Concrete is, at an imperfectly elastic 
material. Because it does not follow the rule of conventional elasticity, it is 
te considered more correct to refer to this property in concrete as the rate of 
_ change of stress with respect to strain. Stress-strain curves are not identical 

_ for the tension and compression areas of loaded, plain concrete beams, as as-_ 

j aa _ Since the shape of the stress. diagram and the location of the neutral ax axis i“ 

when the elastic limit in tension is exceeded, the stress-strain re- 

_ lationships in tensionand compression fibers of a plain concrete beam are un- RC: 

= equal when failure loadis approached;and the compression fibersare stressed _ 


- to about 12% of panei at the instant of failure. This has been illustrat- — 


Preset with increase in depth, resulting in differential ratios of compressive 
_ (Eg) to tensile (Et) moduli of elasticity for beams of varying — The ratio 
exceeds unity, for conventionally loaded beams. 
Apparent “plasticity” in plain concrete subjected to affords relief 
: a elastic peak stresses in the extreme tension fibers, accompanied by local — 4 
= destruction of cohesion. It would be interesting to know whether or not the 7 
“ author plasticity as synonymous with inelasticity. 4 
. NORMAN ABRAMSON. '—Many of the most well-known and widely used ; 
a ‘eiecae of analysis of bending of beams in the plastic range involve the as- a 
sumption that the neutral axis remains fixed with respect to the beam cross J 
However, i in many practical applications the beam section is unsym-_ 
_ metrical with respect toa principal axis, or else the plane of the applied loads 
is not parallel to a principal plane, which, in either case, leads to significant _ 
5 _ twisting (rotation of the cross section) o of the beam, and resultant motionof the — 
a neutral axis, as the loading increases. In fact, this rotation of the cross sec- 
is often the primary factor in limiting the ultimate strength. 
i et When end moments, not parallel to a principal plane of bending, are applied 
>both — 
a... Fig. 8 of the paper “Effect of Depth of Beam Upon the Modulus of ba saath of ‘Plain no 
October, 1959, by Aris Phillips. 
Dir., Dept. of Applied Mec 


— 
— 
— 
‘| 
— 
— 
a 
_ highly stressed, their modulus of elasticity tendsto remain substantially recti- [i 
| 


“Apri, 1960 
. _ principal axes, and will also twist. . When the twist at the mid- span is very i 
° the applied couple. However, as curvature and twist increase, the resisting Bh 
my couple must rotate with respect to the cross section, which in turn requires - 
rotation of the neutral axis. Also, as bending progresses and the outer cor- 
ners of the section enter the plastic condition, it can be seen that for the re- 
. - sisting couple to remain in the same plane, the neutral axis must rotate. This a 
ee is because the “soft” corners cause the points of application of the resultant — 
compressive and tensile forces to move toward the section centroid, and, there- 
fore, if the plane of the couple is to remain fixed with respect to the cross sec- q 
tion, the neutral axis must rotate in order to provide for the readjustment of z 
tress distribution. o . As similar discussion applies also for beams of un- = 
ae ymmetrical cross section, | even if the plane of ™ applied couple is parallel 


ia In view of the practical importance, then, of aed motion of the neutral axis 
Y ea the bending of beams in the plastic range, the author’s present paper | deal- 
with analytical aspects of the problem, even though restricted to beams with 
sections having a plane of symmetry with respect to a axis, is 
ic ~ ‘Turning again, however, to the more complicated problems discussedabove, 
a it may be of more or less direct interest to mention, even if only very briefly, 
some experimental observations4;5 relation to motions of the neutral axis when 
.- loading is unsymmetrical and the beam cross sections are either symmetri- _— 
or unsymmetrical (rectangular, I and L-shape). 
Rectangular Beam.—Asone of a series of tests,a 3/8 in. by 1 in. aluminum 
alloy 24 S-T, rectangular beam with a length of 40 in., , was subjected to bend- 
ing with the plane of the applied couple at an angle of 14°- -56" to the major | 7 
axis.4 The change in location of the axis was obtained graphi- 
; _ The neutral a axis rotated bs 


: bending moment of slightly o over 2,000 in.-lb, one in a rotation of the neu- 
Br _ I-Section Beam. inven ec beams of 40 in. . total length with an 18 in. span | 
rotation, and horizontal ‘components of ‘deflection -span. 

_ For loadings with the plane of the applied couple at both 15° and 30°to the major 

‘a principal axis, it was found that the flanges rotated in opposite directions with _ 
_ - respect to the web, with angles up to 3-1 /2°, while the web rotated toangles up be 3 
a 8°, but the neutral axis showed only small rotations of adegree ortwo. | 
_L-Section Beam.—Studies were also carried out with 24 S-T aluminum nalloy — 
2 section beams® of similar size as the I- beams (1 in. by 2 in. overall Lal As, 

7 _ tangular cross section) in which the plane of loading was at angles of 0°, 30° oo 


60°, and 76°- 07" the major principal axis. rotation of the 


a a. * “Stability of Rectangular and I Section Beams in the Plastic Range,” by H. A. i 
a liams, Tech. Report No. 8, Contract AF W33- -038 AC-16687, Stanford Univ., June 1, 1949. 
5 “An Investigation of the Bending of Angle Beams in the Plastic Range,” by H. N. 
a Journal of the Aeronautic Sciences, Vel. 22, No. 12, pp. 818-828, December, 


— 

— 
— 

— 

~ 

4 

a 

4 

4 

2 = 

— 

— 

_ 

> A. — 

4 


ne 


tat 


uy 


iderable in 


EM 3° to 4°. Again, the n ction itself ro- 
1e order of 3° to 4°. the cross se dto- 
axis was of the the cross section al axis also 
— = rotate less with Feape ri inal position. Tas on or der 
— — war em to indica ften exhibit an § 
sults would seem to a whole, may often exhil te 
erimental re section as a who to bending in 


sion c. Meinais, Stephen Likuan Tsai, and James R. Sims 
BAYLISS Cc. STEPHEN LIKUAN TSAI,2 an and JAMES R. SIMS,° F. 
ASCE.—Thewriters have ‘been very interested in this informative paper on the 
experimental analysis of flat plates using the Ligtenberg method, since simi- _— 
: lar experiments have been in progress at Rice Institute. As the author points — 
out, two methods of experimental analysis using reflected light have been de-— 
veloped. The one reported in this paper and developed by F. Ligtenberg* 
-_ determinest the surface slopes by a double exposure technique. Curvatures are 
then: obtained graphically. An earlier technique developed by M. Dantu®,6 ob- 
_ tains the image strain which is proportional to the bending strain of the plate. : 
‘The method used by the Presan Corporation is essentially the sa same that 
Experiments have been conducted by the w writers 1 using -a circular | plate cut 
SS sheet. If the surface is carefully polished to remove the © ba ; 
slightest imperfections from the surface, good photographs of the reflections 2s 
from the ruled screen can be obtained. Although somewhat better photographs a 
_ have been obtained using plastic models vacuum coated with aluminum, metal A 
7 models are preferred, since a constant value of the modulus of elasticity can 
: be used for | metals. _ As explained ” Mr. . Bradley, the modulus of elasticity of 
Obtaining sufficient fringes c on the photographs of the model for a reasonable — . 


Sprain yar me the thickness c of the | plate s so that the usual bending theory of 
oe thin plates applies, seems to be a problem in the apgtication of the Ligtenberg 
method. The low modulus of elasticity of magnesium sheet, approximately 

A“ 6.5 x 106 psi, makes possible the use of models 1/8 in. in thickness, eat wien f 
desired be obtained with conveniently loads. 


October 1959, William A. 
Graduate Student, Rice Institute, Houston, Tex. Tex. 
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ing to compare the author’s ex-_ 
srimental results with an approximate solution 


which has been used by the writersfor a square, 
uniformly loaded, clamped plate. Values for 
‘Gea and deflections obtaintd | by the author _ 
with Ligtenberg method, and values obtained by 
writers using the Rayleigh-Ritz method are 
plotted in Fig. M2. . Fig. M2 is similar to nd F 


~Ritz method are plotted in place of the author’ 
- The problem of a square, uniformly loaded 


clamped plate was solved applying the Rayleigh- 


method to thin plate theory. FIG, Ml 
en ergy of the system is expressed | by8 


Taking only two parameters, Ay and Ao, their values were determined by mini- 
the potential energy and were found to be Hine! 7 


y 

t 


To check the accuracy of this solution, tl the m maximum 


deflection in the plate, 

013937 q a 
= — 


i answer | is 1. 1.0% greater than the value obtained wy! I. A 


_ 8 “Applied Elasticity, by Chi- Teh | Wang, Mc Graw- Hill Book | New York, 1953, 
- 9 “The Calculation of Maximum Deflection, ‘Moment, and Shear for Uniformly Load- 


Rectangular Plate With Clamped Edges,” by J. Wojtaszak, Journal of 
Mechanics, Vol. 1937, p. A-174. 
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i very close agreement between the experimental results and ee. energy 
method seems to indicate that, in some cases, the experimental method may 


better results than finite difference method. 


since this value was used by the author. 
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BEHAVIOR OF OF STRUCTURES DURING RES DURING RARTHQUAKES! 


"Discussions Clarence J. Derrick nn Glen V. Berg 


in his “Conclusions: A 
“Tt is more reasonable to take the point of view that the design should — *, 
be such that the structures will survive the more frequent, ‘moderate 
4 ground motions without damage but in the more rare eventof verystrong ; 
ground motion, would be tolerated so longas it was not 
Implementation r requires two definitions: 
Of the circumstances which constitute “hazard to life and | limb. i id 


It is implied that the El Centro > earthquake « of May 18, 1940, n may be used as s: 
standard for maximum “very strong ground motion.” The station at which this ; 
earthquake was observed, instrumentally, was locatedon a “foundationcolumn” 
"i of “ “sedimentary rock.” There is reason to believe, from Gutenberg’s studies 2 7 
that for motion at periods greater than 0.6 sec, the surface effect upon “deep 
alluvium?” might have been considerably more violent. 

a _ ‘The “averaged spectra” " were produced by composition of records from dif- _ 
ferent stations located upon different “foundation columns.” This procedure 
tends to obliterate peculiarities of ground motion which Gutenberg attributed — 

‘= foundation influence. At periods greater than 0.5 sec, , these variations may 


‘It is implied that the values of the “averaged spectra,” multiplied by the = 


i “prescribed factors will reproduce actual spectra. ~ comparison for the El : 
Centro: earthquake of May 18, 1940, at 20% of “critica al” damping shows that a a 
; nie discrepancy exists between the deduced values and the published 4" 
_ spectra for the N-S component. 3 If the standard maximum is assumed as the 
resultant of the N- -S and E-W components, the deduced values are about 2/3 of Rs: | 
the published values. It would seem more reasonable to establish the absolute 9 
maximum, and then make allowances for increased strain, rather than to 
in the maximum conditions and increases in ‘strain allowances. 
a 2 “Effects of Ground on Earthquake Motion,” by B. Gutenberg, Bulletin of the Seismol. 
« Analysis of the Taft Accelerogram of the Earthquake of 21 July 1952, ' ” by G. 
‘Housner, Earthquake Research Lab., California Inst. of Tech., September, 1953. Ace 


celeration Spectrum” for El Centro 5/18/40, N-S, Fig. 11, p. 15; “Velocity Spectrum” 
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of an earthquake which actually has occurred. 
bia The “energy” analysis proposed by Mr. Housner is a generalized solution 

of the problem of response by a “damped” resonator, subjected to forced stimu-_ 
lation. In an actual design situation, the problems are specific. The engineer 

is not designing an average structure for average conditions but a particular 

structure for maximum adverse | conditions. The assumption that the distr i- 
| bution of mass does not influence the amount of work performed by an earth-_ 


Quake, is open to 


probable ‘performance of “modern type,” multi- -story, buildings, 
“subjected to “very strong ground motion,” may be deduced. Blume’s vy oe ; 
son of computed and observed free period, for the Alexander Building,4 indi- 
cates that much of the rigidity of older type multi-story buildings derives om 
‘elements no longer included in modern construction. In these older type build- 


- positive evidence on which to base either an analyticalor an empirical solution 
involving “ plastic action” or r “damping” in “modern type,” multi-story buildings. "7 


_ For _— situations, ‘the. case for “reasonable prudence” has been stated: © 
is also the opinion of the writers that all points of earthquake 


resistant design that cannot be established on the basis of reliable em- 
pirical data or incontrovertible analysis be treated with caution 


and that the design should be on the safe side.” 


With respect to “damping,” Jacobsen’s opinion is is illuminating: “Frictional 
in composite structures present a complicated picture, since frictional 


limited ranges of distortion .. .,” ‘and, “It appears relatively certain that the 
&g only practical way of taking: frictional forces into ) consideration i is to use tl the 
concept of equivalent viscous damping forces, an approximation which is very _ 
satisfactory in linear and ‘almost linear’ cases, but in highly non-linear cases 
the concept gives only crude approximations. 


we e Until more is known about the performance of “modern-type,” multi-story | 
buildings, subjected to “very strong ground motion,” the conservative designer 4 
has no alternative but limitation of distortions to “linear” or “almost linear” 
levels. . Many designers are already doing this. However, in some designs, the 
limitation is predicated upon shears deduced from “seismic factors” prescribed 
by building codes. Since shears so deduced usually are only fractions—from 
one-fifth to one-third—of the actual maximum dynamic shears, such limitations ; 
The assurance of public safety involves more than of collapse in 
“the structural framework and fall of debris. Ina disaster, such as an — 
quake developing “very strong ground motion,” the hazards of fire and a 


S 4 «Period Determinations and Other Earthquake Studies of a Fifteen- -Story a 
aed J. A. Blume, Proc. World Conf. on Earthquake Engrg., Berkeley, Calif., —* Wate 
5 -° Discussions of Paper 2514, by R. R. Martel, G W. Housner and 3. A. Alford, “a 
‘- ~ 6 *Frictio “Frictional Effects in Composite Structures Subjected to Earthquake Vibrations,” 


by S. Jacobsen, (A Technical under California Standard 
No. with Stanford ity 
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See §=6frame-work. What the San Francisco earthquake of 1906 would have done to a “a = 
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are extreme. The s framework must be capable, not only maintain- 


_ ing itself but of preventing other safety devices from becoming inoperable. — 
- ‘Fire separation diaphragms must remain intact to prevent the spread, not only - “> 
of fire but also of smoke and other asphyxiating gases. To prevent panic, ele- — ix 
-vator machinery must remain in operation, elevator doors freefrom jamming, | 
By and stairs and other exits unblocked. To assume that interior partitions, stairs, 
‘ete. , will be wrecked, paint 7 framework to distort plastically, is to beg © 
the building is to insure overall it is probable 


that frame stresses will not be critical. 


_ GLEN V. BERG, ,7 F. ASCE. ont the the effect of damping, Mr. Housner : states ee 
that “Moderate amounts of damping will give appreciable reduction of stresses. 
_ It is also clear that large amounts of damping will give large reductions in the 7 
_ forces and stresses in a structure.” A simple characterization of the quantita-_ 
tive effect of damping upon the forces and stresses" can be obtained d through 
a - Some strong-motion accelerograms show a long sequence of acceleration 


4 pulses of more or less uniform magnitude, and others show a brief series of ‘a 


strong pulses followed by a long series of relatively weak pulses. ‘¢ One might 
idealize the former as a random function that feeds energy into a structure at 
wy constant rate, and the latter as a single ganeeaeh Any real earthquake should uf 
fall somewhere between these two extremes. 
a. _ Responding to an earthquake of the first type, a simple ee — 
up essentially a sinusoidal response tuned to its own natural frequency, and = 
the amplitude would build up until the energy dissipated per cycle balanced the ~ 
energy input. For a given rate of energy input, the maximum amplitude (and 
stress) attained by the structure would be proportional to. vi/v, where V is the 
a single impulse would cause a simple structure to respond in an exponen- 
tially decaying sinusoidal 1 motion. For a given impulse the maximum displace- 
‘ment (and (and stress)attained i by | the structure wo would t be i teamed e (3 


"which is, is, approximately, 1 - -1. 57v + +2. 2, 


‘ideal cases. Displacement spectra have been computed by the writer for the 
_ strong-motion earthquakes of El Centro 1934, El Centro 1940, Olympia 1949, 


and Taft 1952. ‘7 B1 shows, for typical systems selected arbitrarily from a 


spectrum calculations, plots of 


mum displacement of the system damping to ve The 2% of critical 
_ damping was chosen as the reference point because it represents about the ng 
= least amount of damping that can exist in a real structure. (Zero was avoided | : 
; because vi/v 1/v becomes infinite there.) Nearly all of the displacements fit be- y 

_ tween the ideal cases, and the results are scattered over almost the entire > 
ra range between them. 

considered as bracketing the quantitative effect of damping, and that trying a 


tie down the characterization within narrower bounds would be futile. iia 


versus where Sq(v) i is the maxi- 
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for an equivalent, elastic, ‘undamped structure.” Current earthquake research 
at of includes an — this — 


vel af of tt. damping 

=. FIG, B1.—EFFECT OF DAMPING ON M XIMUM STRESS 


for elasto- ~plastic systems. Results to date indicate that t the s 
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Pa im] for an isolated case, ld err considerably 


COM MENTARY ‘ON PLASTIC DESIGN IN STEEL: ADDITIONAL =e 


Discussion by M. N. Fialkow 


~ 


~ 


N. FIALKOW, 1 F. _ ASCE.- e expression in n paragraph 6. 2, , of the Com- 

_ mentary for the maximum ratio of section depth to web thickness, which i is ade- 

- quate to insure against premature plastic buckling of the web for the case of. : 

- combined bending and compression, warrants further consideration ar and me 
cussion. This expression for -thickness 


100 3 (2 < o2 oat) 
-is shown in the paper to notd for € max/€y = 4.0 and, i er furnishes a 
safe ratio for use when € max/e y $4.0. The recommendation for its adoption 
by the Committee, therefore, implies that under “ normal” circumstances, plas - x 
fic occur without € max/€y exceeding 4 in any of the plastic 
hinges. This recommendation is not concurred by the writer since, as will be “a 
shown below, this value, (€max €y= 4), is frequently exceeded in normal, static, 

= design problems, its use imposes undue limitation on plastic deflection ~ 
needed in cases of dynamic loading, and it is not consistent with the corre- 7 

sponding strain ratios permitted, in beam flanges during hinge rotation, by the — 
ratio of flange width to thickness recommended by the Committee. In view of ;: 
_ the interdependence between the maximum possible strain and the section depth- a { 
_ web thickness ratio, the relationship of such items as beam- -shape factor and ~ 
computational method to this maximum strain is also noted herein. 
; ¥ a. Several instances of plastic collapse under static loading, which are not 
; 7 5 considered “abnormal” but which require that €max/€y exceed 4 inone or more — 


of the plastic hinges, are cited below: = = | 
Example 1.—Fixed- End Beam with at the ‘Third Point 
; _ (Fig. Fl)—Beedle, in Reference 1.2, on pages 202, 203, has solved this prob- F 
a _ lem and finds that, at the left support, oa, the rotation per unit length in the 
oa y an hinge just prior to collapse, is 8 times as great as ¢p, the rotation per unit 
length up to the elastic limit. 9 = € y0/Yo; it follows that € max/€y = 8. 
q by this analysis, the hinge length has been lh as one eighth the distance from. ; 
.¢ point of of contraflexure to the fixed end, which holds for a wide elie ve { 
with average factor loaded by a concentrated load. 
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17. 4x 10 radians per in. 


“4 
%p =. = - x -5 radians per in 
To provide in design, for large dynamic short duration, utiliz 
sien of the large energy absorbing capacity of structural components Ei 


_ deformed into the plastic zone is often necessary to achieve practical struc- — 


aa (1) Keeping the value of €max/Ey at 4 or less in flexuralmembers would _ 
- impose a severe restriction on the amount of plastic deflection and, therefore, 
on the ability to absorb energy loads. Thus, in the case of either a simple or — 
fixed- -end | beam subjected toa concentrated dynamic load at the midpoint, 
- computations indicate that the maximum deflection would be limited by this a 4 ites 


criterion to 2.5 times the elastic deflection. This contrasts with test data, for 
such beams, which indicate the ultimate deflection at | collapse to be 26 times» 4 
oe (2) Even if the design were developed on the basis of the above limit on Jf 
- €max/e y» the approximate nature of the assumed design overpressure and shock 
loading and the possibility of occurrence of loadings more severe than those i. 
assumed for design make it highly desirable that the structure be capable | of 

deformations beyond those computed for the assumed loading. 

q c. The Committee recommends that for beams in bending the maximum — 
= af ratio of flange width to thickness, be limited to 17. This permits the develop- 2 

a _ ment in the flange of strains to the point of strain- -hardening which corresponds _ 

2” to a ratio of €max/€y of 12 for A7 type structural steel. Reference 6.27 re- | 
_# ports that tests on beams with these flange dimensions confirm the et 
Paes rotation capacity with the further advantage of indicating no sudden drop in load 

capacity for strains beyond the point of strain hardening. For consistent 
i “ey design the critical web dimensions should permit maximum ‘strains e equal | to 

those in the flange, that is, to the point of strain- hardening. 
i 4 as In view of the preceding, it appears to the writer that the euseetehen eet 

an above for maximum ratio of section depth to web thickness for beams in com 
bined bending and compression should be revised. a _ Adoption of an expression 
a7 on a maximum strain in the web, at its junction with the e flange equal to 9 
that of strain-hardening, is believed justified;this corresponds to €max/ey = = 
.. To insure proper application of the Committee’s recommendations for sec + 
x properties, determination of the maximum strain in the hinge is necessary. 


— §6GH = 0.0002. Based on the characteristics of thisframe at collapse,as shown [i 
7 -_ - in Fig. 1.8 of Reference 1.5, the location of the point of contraflexure in the a q ass 
— 
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in Reference 12. this n on. is as the av average in 
the hinge, that is, the total hinge angle divided by the hinge length. The hinge y 
length is taken to be the length of beam over which the moment is greater than 
the yield moment, and is, therefore, dependent on the beam-shape factor and 


the type of loading. For the case of concentrated loads, as in Examples 1 and i 


"a 68=—s- from the point of contraflexure to the fixed end, as the shape factor varies be- 
om = from th t of contrafl to the fixed end, as the shape fact be- 
: tween 1.10 and 1.18, the usual range in shape factors for wide flange sections. 


a above, the hinge length varies between 0.09 and 0.15 times the distance, © 


7 
Thus, in lieu of €m ax/€y y= 8,as computed in Example 1 for a beam with shape 


—- of 1.14 , this T Tatio would | equal 11 for a beam with shape factor of 1. 10. 


. For uniform loading, the hinge length would vary between 0.3 and 0.4 of the 
- beam length for simple ‘spans, and between 0.05 and 0. 08 times the beam length | 
for cantilever beams. A different determination of hinge length is made in 
; Reference 1.5 on page 103, where Hodge uses one half of the above definition — 
7 of hinge length in order to compute the “average” value of the strain in alll 


‘4 hinge. It is suggested that the Committee’ s further reports include recom-_ 


mended procedures for determination of hinge strains 
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"DIVISION ACTIVITIES 
ENGINE ERING MECHANICS DIVISION 


of the American Society of ( Civil Engineers 


PROGRAM FOR PURDUE CONFERENCE > 


May 5, Thursday Morning - Registration 
5, Thursday Afternoon - Members and Frames 


__ “Stability of Frames in the Presence of Primary wigan Moments”, 


“Analysis of ‘Frames with Linear Behavior Subjected 
Deformations” » A. Ang, U.of Ml, 
“Restrained Columns”, M. Ojalvo, ‘Lehigh 


2 


“Ultimate Strength of ‘Over -Reinforced B. Kriz and S. L. Lee, 


Bearing Capacity of Floating Ice Sheets” ,G.G. . Meyerhof, Nova — 


“Load Replacement in Plastic Analysis”, B. C. Ringo, Mich. State. 


“Electrical Energy Analogs of Vibrating Beams” L. Ryder, Republic 


Stress Waves and Shock in Locking Media”, M. G. Salvadori and R. 

Skalak, Columbia U. and P. Weidlinger, engineer, N. Y. 
“Energy Relations for Elasto-Plastic Systems Responding to Earthquakes,” 


G. V. Berg and S. S. -Thomaides, U. of 


aa _ Note.—No. 1960- 12 is part of the copyrighted Journal of the Engineering Mechanics _ 
Division, Proceedings of the American Society of Civil Engineers, Vol. 86, No. EM - - 
1960 by ‘the American Society of Civil nip 


1960-12 

— 
Fr. q 
— 
Oe 
— — 
— 
= 
— 
— 
— 
— 
a 
— 


May 5, 5, Thurs. Evening - Banquet 

_ P.I., Valparaiso U., Purdue U., and Ind. Section of ASCE will be held at ee 
Purdue Memorial Union Bldg. Prof. M. G. Salvadori will ‘Speak meni. 

 (D. H. Pletta, Chairman, EMD, ill a meeting of the Executive Com 
for May 7 at Purdue. Members with items for the should 


_ dress remarks remarks to E. Wenk E. Wenk, Jr. Sr. Secretary.) 


PROGRAM FOR BROWN SYMPOSIUM 
‘An advance program has been distributed for the Second Symposium on 


reau of Ships. 
and Impact Loading— David Taylor Model Basin. 


y ‘Atomic ‘Theory of Plasticity and Fracture of Crystalline Solids” = 
ay 


One - Current Research Contributions—R. S. Davis, Arthur D. Little Inc.; and 
Current Research Contributions—D. C. Drucker, Brown U.;R.M. 
Haythornthwaite, U. of Mich.; W. + Hu, Pa. State and A. Phillips, 


Pay “Theory and of Work Materials” 


Current Research Contributions — —B. H. J. J. Greenber 
April 6, Wednesday Afternoo 
‘Current Research Contributions — -M. E. Lunchick, David taylor Model 
2 
April 6, “Wednesday Evening - 
4 April 1, T Thursday Morning 


W. S. Dorn, IMB Corp. 
‘gigi of Boundary Value Problems, Plates and Shells” 
Basia; J. H. Weiner and H. G. Landau, Columbia U. 
“Dynamic Loading and Plastic Waves” 
rveys—N. Cristescu, Rumanian Academy, N. M. Newmark, A. cial 
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Veletsos and Ss. Suteliffe, U. of ‘Tll.; T. H. i. H. Pian, MIT; ‘ELA. 


and J. Plass, U. of Tex. 


‘April % “Thursday day Afternoon 


s—J. Heyman, Cambridge U.; and L. Lehigh. 


Contributions on Current —W. Bureau of Weapons, 


SECOND CONFERENCE ON ELECTRONIC COMPUTATION 
_ As a follow-up to the initial announcement, the Committee on Electronic — 
i Computation of the Structural D Division announces an en enthusiastic response to _ 
the invitation for p papers to be p presented at the conference on September 8 and 
9, 1960, at the Pittsburgh Hilton Hotel. Approximately 75 abstracts have been 
received from which 25 to 30 papers will be selected for the conference. = 
Final selections will be made after receipt of the complete manuscripts, _ - 


which are due April 1: | The committee is anticipating an attendance in ex- : 
‘SIXTH CONGRESS, IABSE 


_ June 27 through July 1, 1960, in Stockholm, Sweden. ~~ | 
Morning— -Registration and meeting of permanent committees 


Basis of Structural Design, F. Stiissi 


Morning—New Developments in Butlding, H. Louis presiding. 


Sue June 30, , Thursday 


Afternoon—The Design of Pile Foundations. | ; 
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H INTERNATIONAL comanans ¢ OF )F APPLIED 
‘The Congress is scheduled for August 31 through September 1960, 

at Strsa, Italy. Apart from a number of invited general lectures, the techni- 
sessions will be held in two sections: 


Section 1—Fluid Dynamics (hydrodynamics and aerodynamics) 


Section 2—Mechanics of solids (rigid-body dynamics, vibrations, elasticity, : 
and theory of structures). 


of is by. the Italian ¢ Committee 
(President: Prof. G. Colonnetti. Secretary: Dr. F. Rolla). For information 
ee accommodations, write to: Dr. F. Rolla, Consiglio: Nazionale delle Ricer- 
che, Uffizio relazioni internationali, Piazza della Scienza Rome, Italy. 
ollowing excerpt from a letter dated January 5, 1960, to Dan H. wee 
Dletta, Chairman, EMD, from Harold Larsen, Manager, , Technical 


“Ay review of the increased publications activity in the 


— obviously has no trouble providing us with a suitable number — 
papers. As a result, we would like on a bi- 


Hansen of MIT rather than John Ss. 


* of New York City rather than C. W. Britzius. 
4 ae Fadum of North Carolina State is the ASCE representative to the 


Ina M. P. the 1961 Committee 
(Herrmann, Suer, Cermak, Goldberg, Throop, Hall, Haltiwanger, and a 
was asked to organize technical programs for the year beginning at the end of ae 


the 1960 Annual Convention and closing at the end of 1961 Annual: Co 
a which will be in New York City. Sel 
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aa -___ at full length and abstracts of papers accepted for presentation at one of the ; 4 
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Engineering Mecha 
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; fe T bers were asked to consider 


1. The Annual Convention which allows EMD to hold ten to twelve half-day 


EMD participation inthreetypesof 


sessions alone or jointly with another division. 


ie’ Regional ASCE Conventions such as the one to be held in Phoenix, _ 

_ 3. Purely Engineering Mechanics Conferences such as the annual ASCE- 

A ASME West Coast Mechanics C 

be held May 


cation in a Newsletter for circulation with the __ 


_ Donald L. Dean, Newsletter Editor, EMD 
‘Department of Civil Engineering 
of Kansas 
Lawrence, Kansas 
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